A BME Matematika Intézet Analízis és Differenciálegyenletek Tanszékének közös Alkalmazott Analízis Szemináriuma 2016. őszén indult Faragó István (Differenciálegyenletek Tanszék) kezdeményezésére az MTA-ELTE Numerikus Analízis és Nagy Hálózatok Kutatócsoporttal együttműködésben. A szeminárium célja, hogy elősegítse egy alkalmazott analízissel (funkcionálanalízis, differenciálegyenletek, numerikus módszerek) foglalkozó kutatói kör kialakítását az intézeten belül. A szemináriummal fórumot szeretnénk biztosítani az alkalmazott analízissel foglalkozó matematikusok és az analízist alkalmazó kutatók számára az együtt gondolkodásra. További cél az érdeklődő hallgatók (MSc, PhD) bevonása a kutatói munkába.

Szemináriumunk 2017-től felvette a Farkas Miklós Alkalmazott Analízis Szeminárium nevet. Ezzel szeretnénk emléket állítani egyetemünk egykori tanszékvezető matematikaprofesszorának, aki elindította egyetemünkön a matematikus-mérnök képzést, és a stabilitáselmélet valamint a biomatematika terén elért jelentős tudományos eredményeivel ill. könyveivel nagyban hozzájárult az alkalmazott matematika erősödéséhez. (English version of the introduction.)

** From the autumn semester of 2017 the talks will be in English on a regular basis. / 2017. őszi félévétől az előadásokat angol nyelven tartjuk. **

Organizers / Szervezők: Faragó István^{1,2,3}, Karátson János^{1,2,3} ,Horváth Róbert^{1,3} ,Mincsovics Miklós^{1,3} (^{1}BME, ^{2}ELTE, ^{3}MTA-ELTE NUMNET)

Prospective speakers:

- 28 February: László Simon (ELTE)
- 7 March: Orsolya Sáfár (BME)
- 14 March: Yiannis Hadjimichael (NUMNET MTA-ELTE reserach group)
- 21 March: - (Spring holiday at BME)
- 28 March: Tamás Insperger (BME)

**28 February (Thursday) 10:15, H306**

László Simon (ELTE)

**Multiple solutions of nonlinear partial functional differential equations**

We shall consider weak solutions of nonlinear elliptic boundary value problems and initial-boundary value problems for semilinear and nonlinear parabolic differential equations with certain nonlocal terms. We shall prove theorems on the number of solutions and find multiple solutions. These statements are based on arguments for fixed points of some real functions and operators, respectively, and existence-uniqueness theorems on partial differential equations (without functional terms).

Previous seminars:

21 February 2018

Gustaf Söderlind (Lund University, Sweden; guest professor at NUMNET MTA-ELTE reserach group)

**Logarithmic norms and Quadratic forms** with applications to Calculus of Variations

The logarithmic norm was introduced in 1958 for matrices, and for the purpose of estimating growth rates in initial value problems. Since then, the concept has been extended to nonlinear maps, differential operators and function spaces. There are applications in operator equations in general, including evolution equations as well as boundary value problems. The logarithmic norm is the extremal value of a quadratic form.

In this talk we outline how logarithmic norms of differential operators can be computed, and how they are related to variational calculus and ellipticity. Thus, while one typically seeks the minimizing function in a variational problem, the logarithmic norm is the corresponding extremal value of the functional associated with a particular symmetrized differential operator. There are also connections to eigenvalue problems for selfadjoint and non-selfadjoint operators. This will also be illustrated with an application to classical singular problems, such as the Bessel equation, as well as to the biharmonic operator.

14 February 2018

**Róbert Horváth (BME)**

**On some qualitative properties of parabolic problems**

In this talk we investigate some special qualitative properties of parabolic problems. At the beginning of the talk we review the remarkable qualitative properties of these problems. Then we will turn to two special properties: the first property says that the number of the so-called $L$-level points, or specially the number of the zeros, of the solutions must be non-increasing in time. The second property requires a similar property for the number of the local maximizers and minimizers. We show that linear equations and some special nonlinear equations fulfill the above properties in the continuous case. We use the maximum-minimum principles in the proof. Then we generate the numerical solution with the implicit Euler finite difference method and show that the obtained numerical solution satisfies the discrete versions of the above properties without any requirements on the mesh parameters. We show also some numerical test results.

13 December 2018

**Research reports of PhD students**

Teshome Bayleyegn: The method of multiple Richardson extrapolation

Benjámin Borsos: Variable preconditioning for strongly nonlinear elliptic problems

Szilvia György: The stability of systems of linear delay differential equations

Rahele Mosley: Reliable numerical modelling of malaria propagation

Bálint Takács: Epidemic models with spatial dependence

29 November 2018

Gábor Stépán (BME, Department of Applied Mechanics)

**Dynamics of hardware-in-the-loop experiments**

A widely used tool of engineering R&D is the hardware-in-the-loop (HIL) experiment. Instead of building the full prototype of a developed machine, only its most critical parts are constructed physically, while the rest of the machine is emulated by means of actuators, sensors, and digital control in between. If the mathematical model of the rest of the machine is available, the control unit can provide a realistic environment for the physically constructed test part. The control, however, introduces digital effects into this system which is originally continuous. From dynamical view-point, the most relevant digital effects are the appearances of delay and zero-order-hold (ZOH). The nonlinear dynamics of the real system and the one emulated by means of the HIL experiment are compared from stability and nonlinear vibrations view-point in case of a brake system where stick-slip phenomenon occurs. The limitations of HIL experiments are identified by means of Hopf bifurcation calculations, numerical simulations and dynamic measurements carried out on the corresponding experimental test rig.

The introductory part of the lecture will present the concept of the ERC advanced grant project and the role of HIL in the development of new milling tools for metal cutting, and the concluding part of the lecture will present the application of the results for the development of HIL experiments in case of high-speed-milling (HSM).

- Presentation and paper

22 November 2018

János Karátson (ELTE, BME)

**Extended finite element methods: a brief introduction**

The finite element method (FEM) is a fundamental tool of the numerical solution of real-life problems based on partial differential equations. In the recent decades, various generalizations of the standard FEM have been developed. A lot of such extensions have been motivated by difficulties, arising in physical or engineering problems, that may be cumbersome to overcome with standard FEM techniques. Such situations are the presence of boundary layers, singularities or discontinuities in the solution, complex and/or evolving geometry of the physical domains etc. The tools of extension of the FEM may be enriching the polynomial approximation space with non-polynomial shape functions, allowing general polygonal/polyhedral cells, or use a boundary-unfitted mesh and restricted shape functions (either to a bulk domain or to a surface). This survey type talk gives a brief introduction to the main ideas of some generalized FEMs that use the above ideas: XFEM, VEM, CutFEM and TraceFEM.

15 November 2018

Petra Csomós (ELTE)

**Operator semigroups, innovative integrators, and flood prevention**

The talk summarizes our results obtained when studying the possible application of innovative integrators (operator splitting procedures and exponential integrators) to optimal control problems. We introduce the innovative integrators considered, and present the advantages and drawbacks of their use. After introducing the abstract optimal control problem having a linear quadratic cost function, we prove the convergence of innovative integrators applied to it. To do so, we briefly summarize the main results in operator semigroup theory, used in the proof. As an example, we treat the linear quadratic regulator problem of the one-dimensional and two-dimensional shallow water equations, and illustrate our results by numerical experiments.

The results are based on the joint works with Johannes Winckler (Tübingen) and Hermann Mena (Innsbruck).

- Related papers:

8 November 2018

Barna Garay (Pázmány Péter Chatolic University)

**Moving average network examples for asymptotically stable periodic orbits of monotone maps**

For a certain type of discrete--time nonlinear consensus dynamics, asymptotically stable periodic orbits are constructed. Based on a simple ordinal pattern assumption, the Frucht graph, two Petersen septets, hypercubes, a technical class of circulant graphs (containing Paley graphs of prime order), and complete graphs are considered -- they are all carrying moving average monotone dynamics admitting asymptotically stable periodic orbits with period 2. Carried by a directed graph with 594 (multiple and multiple loop) edges on 3 vertices, also the existence of asymptotically stable r-periodic orbits, r=3,4,... is shown.

25 October 2018

Mihály Pituk (Departme**nt of Mathematics, University of Pannonia, Veszprém, Hungary)**

**Ergodicity in nonautonomous linear ordinary differential equations**

The weak and strong ergodic properties of the positive solutions of nonautonomous linear ordinary differential equations will be considered. We will show that if the coefficient matrix function is bounded, essentially nonnegative and uniformly irreducible, then the normalized positive solutions are asymptotically equivalent to the Perron vectors of the strongly positive transition matrix at infinity (weak ergodicity). If, in addition, the coefficient matrix function is uniformly continuous, then the convergence of the normalized positive solutions to the same strongly positive limiting vector (strong ergodicity) is equivalent to the convergence of the Perron vectors of the coefficient matrices.

This is a joint work with Professor Christian Pötzsche (Alpen-Adria University Klagenfurt, Austria).

18 October 2018

Stefan Filipov (Department of Computer Science, Faculty of Chemical System Engineering, University of Chemical Technology and Metallurgy, Sofia, Bulgaria)

**Successive application of the linear shooting method for numerical solution of nonlinear two-point boundary value problems**

The work to be presented concerns numerical solution of nonlinear two-point boundary value problems. We show that the quasi-linearization method (Newton method on operator level) can be used as a basis to derive (i) the FDM with Newton method and (ii) the shooting by Newton method. The same relation holds for the Picard and the constant-slope methods. Based on these results, we propose (i) a replacement of the FDMs for nonlinear problems (the relaxation methods) by respective successive application of the linear shooting method and (ii) a shooting by Picard method (shooting-projection method). We discuss the advantages of the proposed approaches and present examples.

11 October 2018

Róbert Kovács (BME, Department of Energy Engineering)R

**Anal****ytical and numerical solutions of dissipative systems**

Our world is not ideal, in reality the processes are dissipative. The framework of non-equilibrium thermodynamics offers lot of possibilities to derive models, constitutive equations that describe the behavior of a dissipative system. The level of modeling is arbitrary, depends on our choice which is reflected by these models. However, they should be remain as simple as possible to be applicable for practical problems in question.

First, the non-equilibrium thermodynamical background and the structure of equations are discussed. Here mostly the parabolic - hyperbolic properties of the resulted partial differential equations (PDE) are emphasized. The role of boundary conditions and their effect on solutions are also presented through different examples. Such example is related to a particular experimental arrangement called heat pulse (or laser flash) experiment that used to detect different dissipative wave propagation phenomena.

In this presentation the way from generating the equations to their solutions for experiments is presented. It covers analytical solution of a PDE for time dependent boundary condition and a particular numerical method that allows us to eliminate certain boundary conditions and related to the specific structure resulted by non-equilibrium thermodynamics.

- Presentation
- Homepage of the Montavid Thermodynamical Research Group
- Related papers:
Numerical solutions: https://www.sciencedirect.com/science/article/pii/S0017931018322002 (arXiv preprint arXiv:1805.01108)

Analytic solutions: https://www.sciencedirect.com/science/article/pii/S001793101831843X (arXiv preprint arXiv:1804.05225)

Theoretical background: https://www.sciencedirect.com/science/article/pii/S0017931014011533

4 October 2018

Miklós Rontó (University of Miskolc)

**New successive approximation method of general non-local boundary value p****roblems**

We propose a new successive approximation technique for the solvability analysis and approximate solution of general non-local boundary value problems for non-linear systems of differential equations with locally Lipschitzian non-linearities. It will be studied the non-linear boundary value problem

$\frac{dx(t)}{dt}=f(t,x(t)),\,t\in \left[ a,b\right] , \Phi (x)=d, $

where $\Phi :C(\left[ a,b\right] ,\mathbb{R}^{n})\rightarrow \mathbb{R}^{n}$ is a vector functional (possibly non-linear), $~f:\left[ a,b\right]

\times \mathbb{R}^{n}\rightarrow \mathbb{R}^{n}~$is a function satisfying the Caratheodory conditions in a certain bounded set $D$, which will be concretized later, $d$ is a given vector and $f\in Lip(K,D),$ i.e. $f$ locally Lipschitzian

$\left\vert f(t,u)-f(t,v)\right\vert \leq K\left\vert u-v\right\vert ,~\text{for all}\ \left\{ u,v\right\} \subset D~\text{ }a.e.t\in \left[ a,b\right] . $

By a solution of the problem, one means an absolutely continuous function satisfying the differential system almost everywhere on $\left[ a,b\right] .$ The analysis is constructive in the sense that it allows one to both study the solvability of the problem and approximately construct its solutions by operating with objects that are determined explicitly in finitely many steps of computation. The practical application of the technique is explained on a numerical example.

27 September 2018

Ilona Nagy (BME, Department of Analysis)

**Two Limit Cycles in a Two-Species Reaction**

Kinetic differential equations, being nonlinear, are capable of producing many kinds of exotic phenomena. However, the existence of multistationarity, oscillation or chaos is usually proved by numerical methods. Here we investigate a relatively simple reaction among two species consisting of five reaction steps, one of the third order. About this reaction we show the following facts (using symbolic methods): the kinetic differential equation of the reaction has two limit cycles surrounding the stationary point of focus type in the positive quadrant. The outer limit cycle is always stable and the inner one is always unstable. We also performed the search for partial integrals of the system and have found one such integral. Application of the methods needs computer help (Wolfram language and Singular) because the symbolic calculations to carry out are too complicated to do by hand. Even to make characteristic drawings is very far from trivial: the inner limit cycle is very small. Some of the methods we use are relatively new, and all of them have never been used in reaction kinetics, although they can be used to have similar exact results on other models.

20 September 2018

Professors Gustaf Söderlind and Carmen Arévalo will present three talks during their visit to Budapest. All three talks deal with adaptive methods of ODEs. Professor Söderlind’s first talk will be held at the Seminar on Applied Analysis of the Department of Applied Analysis at ELTE university (details below), while his second talk together with Professor Arévalo’s talk will be presented at the Miklós Farkas Seminar. Although professor Söderlind’s talks form a complete presentation together, yet they can be followed separately.

Gustaf Söderlind (Lund University, Sweden)

**Adaptive Time-Stepping. Part II. Time transformations applied to reversible Hamiltonian dynamics and weakly dissipative systems**

In the second talk on time step adaptivity, we focus on the special needs of conservative dynamical systems. This includes Hamiltonian problems, and weakly dissipative systems. In integrable Hamiltonian problems, the mathematical solution is time reversible, which precludes the use of classical controllers, which adapt the step size to manage the error observed in previous steps. Instead, a time reversible tracking algorithm is developed, which allows full reversibility of the adaptive computational process. This is shown to preserve first integrals over long times, and even improves the accuracy over constant step size symplectic integrators. We demonstrate the procedure in two examples from celestial mechanics, and then proceed to demonstrate how a similar approach can be combined with splitting methods in weakly dissipative systems. The latter approach has been put to effective use in rolling bearing dynamic simulation.

The first part of the talk:

Adaptive Time-Stepping. Part I. How control theory and digital filters enhance the performance of IVP solversSeminar on Applied Analysis, 17 September 2018, 14:00, Department of Applied Analysis and Computational Mathematics, ELTE, Pázmány P. 1/C, room 3.719

In two talks we will describe a systematic approach to time-stepping adaptivity in the numerical integration of initial value problems. The first part will deal with the proper construction of discrete controllers that keep the local error equidistributed along the solution trajectory. The controllers will be based both on classical constructions such as proportional-integral controllers, and on modern digital low-pass filters, which are designed to maintain a smooth step size sequence. We will demonstrate that a regular step size sequence is of benefit as it increases the computational stability — small changes in the error tolerance will only produce a small change in the accuracy of the computed solution. While computational speedup is often only modest, the quality of the computed solution as well as the reliability of the software can sometimes be remarkably improved. Presentation

Carmen Arévalo (Lund University, Sweden)

**An intrinsically adaptive formulation of multistep methods**

Multistep methods are important tools for solving ordinary differential equations with initial conditions. In order for these methods to be efficient they must be adaptive, that is, they must allow the choice of an appropriate step-size for each integration step. We present a comprehensive way of formulating multistep methods that is adaptive by construction and show how this methodology can be applied to particular situations. We also show the application to strong stability preserving methods, used to solve ODEs arising from the semi-discretization of time-dependent partial differential equations (PDEs), especially hyperbolic PDEs with shocks. We finally demonstrate how we apply the formulation to differential algebraic equations (DAEs), an ODE system coupled with algebraic constraints.

Presentations in 2017/18, 2016/17