scholarly journals Applied Koopman Theory for Partial Differential Equations and Data-Driven Modeling of Spatio-Temporal Systems

Complexity ◽  
2018 ◽  
Vol 2018 ◽  
pp. 1-16 ◽  
Author(s):  
J. Nathan Kutz ◽  
J. L. Proctor ◽  
S. L. Brunton

We consider the application of Koopman theory to nonlinear partial differential equations and data-driven spatio-temporal systems. We demonstrate that the observables chosen for constructing the Koopman operator are critical for enabling an accurate approximation to the nonlinear dynamics. If such observables can be found, then the dynamic mode decomposition (DMD) algorithm can be enacted to compute a finite-dimensional approximation of the Koopman operator, including its eigenfunctions, eigenvalues, and Koopman modes. We demonstrate simple rules of thumb for selecting a parsimonious set of observables that can greatly improve the approximation of the Koopman operator. Further, we show that the clear goal in selecting observables is to place the DMD eigenvalues on the imaginary axis, thus giving an objective function for observable selection. Judiciously chosen observables lead to physically interpretable spatio-temporal features of the complex system under consideration and provide a connection to manifold learning methods. Our method provides a valuable intermediate, yet interpretable, approximation to the Koopman operator that lies between the DMD method and the computationally intensive extended DMD (EDMD). We demonstrate the impact of observable selection, including kernel methods, and construction of the Koopman operator on several canonical nonlinear PDEs: Burgers’ equation, the nonlinear Schrödinger equation, the cubic-quintic Ginzburg-Landau equation, and a reaction-diffusion system. These examples serve to highlight the most pressing and critical challenge of Koopman theory: a principled way to select appropriate observables.

Author(s):  
Nkosingiphile Mnguni ◽  
Sameerah Jamal

Abstract This paper considers two categories of fractional-order population growth models, where a time component is defined by Riemann–Liouville derivatives. These models are studied under the Lie symmetry approach, and we reduce the fractional partial differential equations to nonlinear ordinary differential equations. Subsequently, solutions of the latter are determined numerically or with the aid of Laplace transforms. Graphical representations for integral and trigonometric solutions are presented. A key feature of these models is the connection between spatial patterning of organisms versus competitive coexistence.


2017 ◽  
Vol 59 (2) ◽  
pp. 167-182 ◽  
Author(s):  
H. Y. ALFIFI

Semi-analytical solutions are derived for the Brusselator system in one- and two-dimensional domains. The Galerkin method is processed to approximate the governing partial differential equations via a system of ordinary differential equations. Both steady-state concentrations and transient solutions are obtained. Semi-analytical results for the stability of the model are presented for the identified critical parameter value at which a Hopf bifurcation occurs. The impact of the diffusion coefficients on the system is also considered. The results show that diffusion acts to stabilize the systems better than the equivalent nondiffusive systems with the increasing critical value of the Hopf bifurcation. Comparison between the semi-analytical and numerical solutions shows an excellent agreement with the steady-state transient solutions and the parameter values at which the Hopf bifurcations occur. Examples of stable and unstable limit cycles are given, and Hopf bifurcation points are shown to confirm the results previously calculated in the Hopf bifurcation map. The usefulness and accuracy of the semi-analytical results are confirmed by comparison with the numerical solutions of partial differential equations.


Author(s):  
A. KARAMI ◽  
H. R. KARIMI ◽  
B. MOSHIRI ◽  
P. JABEDAR MARALANI

Theoretical modeling of dynamic processes in chemical engineering often implies the numeric solution of one or more partial differential equations. The complexity of such problems is increased when the solutions exhibit sharp moving fronts. An efficient adaptive multiresolution numerical method is described for solving systems of partial differential equations. This method is based on multiresolution analysis and interpolating wavelets, that dynamically adapts the collocation grid so that higher resolution is automatically attributed to domain regions where sharp features are present. Space derivatives were computed in an irregular grid by cubic splines method. The effectiveness of the method is demonstrated with some relevant examples in a chemical engineering context.


2011 ◽  
Vol 2011 ◽  
pp. 1-37
Author(s):  
Jan Harm van der Walt

The method of order completion provides a general and type-independent theory for the existence and basic regularity of the solutions of large classes of systems of nonlinear partial differential equations (PDEs). Recently, the application of convergence spaces to this theory resulted in a significant improvement upon the regularity of the solutions and provided new insight into the structure of solutions. In this paper, we show how this method may be adapted so as to allow for the infinite differentiability of generalized functions. Moreover, it is shown that a large class of smooth nonlinear PDEs admit generalized solutions in the space constructed here. As an indication of how the general theory can be applied to particular nonlinear equations, we construct generalized solutions of the parametrically driven, damped nonlinear Schrödinger equation in one spatial dimension.


Sign in / Sign up

Export Citation Format

Share Document