Source author record

Dirk Lebiedz

Dirk Lebiedz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2026arXiv

Approximating normally attracting invariant manifolds using trajectory-based optimization

The numerical simulation of realistic reactive flows is a major challenge due to the stiffness and high dimension of the corresponding kinetic differential equations. Manifold-based model reduction techniques address this problem by projecting the full phase space onto manifolds of slow motion, which capture the system's long-term behavior. In this article we study the trajectory-based optimization approach by Lebiedz (2004), which determines these manifolds as minimizers of an appropriate entropy functional. Similar to other methods in this field, this approach is based on physical and geometric intuition and was tested on several models. This article provides a rigorous explanation for its effectiveness, showing how it approximates nonuniformly normally attracting orbits. It also outlines how the method can be utilized to approximate nonuniformly normally attracting invariant manifolds of higher dimension. Throughout the article we use a coordinate-free formulation on a Riemannian manifold. This is especially useful for systems subject to nonlinear constraints, e.g., adiabatic constraints.

preprint2022arXiv

On Differential Geometric Formulations of Slow Invariant Manifold Computation: Geodesic Stretching and Flow Curvature

The theory of slow invariant manifolds (SIMs) is the foundation of various model-order reduction techniques for dissipative dynamical systems with multiple time-scales, e.g. in chemical kinetic models. The construction of SIMs and many approximation methods exploit the restrictive requirement of an explicit time-scale separation parameter. Most of those methods are also not formulated covariantly, i.e. in terms of tensorial constructions. We propose an intrinsically coordinate-free differential geometric approximation criterion approximating normally attracting invariant manifolds (NAIMs). We translate some ideas behind existing approximation approaches, the stretching based diagnostics (SBD) and the flow curvature method (FCM) to tensors of Riemannian geometry, specifically to spacetime curvature in extended phase space. For that purpose we derive from flow-generating smooth vector fields a metric tensor such that the original dynamical system is a geodesic flow on a Riemannian manifold. We apply the resulting method to test models.

preprint2021arXiv

Generalized Liénard systems, singularly perturbed systems, Flow Curvature Method

In his famous book entitled \textit{Theory of Oscillations}, Nicolas Minorsky wrote: "\textit{each time the system absorbs energy the curvature of its trajectory decreases} and \textit{vice versa}". According to the \textit{Flow Curvature Method}, the location of the points where the \textit{curvature of trajectory curve}, integral of such planar \textit{singularly dynamical systems}, vanishes directly provides a first order approximation in $\varepsilon$ of its \textit{slow invariant manifold} equation. By using this method, we prove that, in the $\varepsilon$-vicinity of the \textit{slow invariant manifold} of generalized Liénard systems, the \textit{curvature of trajectory curve} increases while the \textit{energy} of such systems decreases. Hence, we prove Minorsky's statement for the generalized Liénard systems. Then, we establish a relationship between \textit{curvature} and \textit{energy} for such systems. These results are then exemplified with the classical Van der Pol and generalized Liénard \textit{singularly perturbed systems}.

preprint2014arXiv

On fundamental unifying concepts for trajectory-based slow invariant attracting manifold computation in multiscale models of chemical kinetics

Chemical kinetic models in terms of ordinary differential equations correspond to finite dimensional dissipative dynamical systems involving a multiple time scale structure. Most dimension reduction approaches aimed at a slow mode-description of the full system compute approximations of low-dimensional attracting slow invariant manifolds and parameterize these manifolds in terms of a subset of chosen chemical species, the reaction progress variables. The invariance property suggests a slow invariant manifold to be constructed as (a bundle of) solution trajectories of suitable ordinary differential equation initial or boundary value problems. The focus of this work is on a discussion of fundamental and unifying geometric and analytical issues of various approaches to trajectory-based numerical approximation techniques of slow invariant manifolds that are in practical use for model reduction in chemical kinetics. Two basic concepts are pointed out reducing various model reduction approaches to a common denominator. In particular, we discuss our recent trajectory optimization approach in the light of these two concepts. We relate both of them in a variational boundary value viewpoint, propose a Hamiltonian formulation and conjecture its relation to conservation laws, (partial) integrability and symmetry issues as underlying fundamental principles and potentially unifying elements of diverse dimension reduction approaches.

preprint2013arXiv

A continuation method for the efficient solution of parametric optimization problems in kinetic model reduction

Model reduction methods often aim at an identification of slow invariant manifolds in the state space of dynamical systems modeled by ordinary differential equations. We present a predictor corrector method for a fast solution of an optimization problem the solution of which is supposed to approximate points on slow invariant manifolds. The corrector method is either an interior point method or a generalized Gauss--Newton method. The predictor is an Euler prediction based on the parameter sensitivities of the optimization problem. The benefit of a step size strategy in the predictor corrector scheme is shown for an example.

preprint2013arXiv

A Numerical Slow Manifold Approach to Model Reduction for Optimal Control of Multiple Time Scale ODE

Time scale separation is a natural property of many control systems that can be ex- ploited, theoretically and numerically. We present a numerical scheme to solve optimal control problems with considerable time scale separation that is based on a model reduction approach that does not need the system to be explicitly stated in singularly perturbed form. We present examples that highlight the advantages and disadvantages of the method.

preprint2013arXiv

An optimization approach to kinetic model reduction for combustion chemistry

Model reduction methods are relevant when the computation time of a full convection-diffusion-reaction simulation based on detailed chemical reaction mechanisms is too large. In this article, we review a model reduction approach based on optimization of trajectories and show its applicability to realistic combustion models. As most model reduction methods, it identifies points on a slow invariant manifold based on time scale separation in the dynamics of the reaction system. The numerical approximation of points on the manifold is achieved by solving a semi-infinite optimization problem, where the dynamics enter the problem as constraints. The proof of existence of a solution for an arbitrarily chosen dimension of the reduced model (slow manifold) is extended to the case of realistic combustion models including thermochemistry by considering the properties of proper maps. The model reduction approach is finally applied to three models based on realistic reaction mechanisms: 1. ozone decomposition as a small test case; 2. simplified hydrogen combustion for comparison with another model reduction method; 3. syngas combustion as a test case including all features of a detailed combustion mechanism.

preprint2012arXiv

A robust optimization approach to experimental design for model discrimination of dynamical systems

A high-ranking goal of interdisciplinary modeling approaches in the natural sciences are quantitative prediction of system dynamics and model based optimization. For this purpose, mathematical modeling, numerical simulation and scientific computing techniques are indispensable. Quantitative modeling closely combined with experimental investigations is required if the model is supposed to be used for sound mechanistic analysis and model predictions. Typically, before an appropriate model of a experimental system is found different hypothetical models might be reasonable and consistent with previous knowledge and available data. The parameters of the model up to an estimated confidence region are generally not known a priori. Therefore one has to incorporate possible parameter configurations of different models into a model discrimination algorithm. In this article we present a numerical algorithm which calculates a design of experiments which allows an optimal discrimination of different hypothetic candidate models of a given dynamic system for the most inappropriate parameter configurations within a parameter range via a worst case estimate. The design criterion comprises optimal measurement time points. The used criterion is derived from the Kullback-Leibler divergence. The underlying optimization problem can be classified as a semi infinite optimization problem which we solve in an outer approximation approach stabilized by a homotopy strategy. We present the theoretical framework as well as the numerical algorithmic realization.

preprint2009arXiv

A variational principle for computing slow invariant manifolds in dissipative dynamical systems

A key issue in dimension reduction of dissipative dynamical systems with spectral gaps is the identification of slow invariant manifolds. We present theoretical and numerical results for a variational approach to the problem of computing such manifolds for kinetic models using trajectory optimization. The corresponding objective functional reflects a variational principle that characterizes trajectories on, respectively near, slow invariant manifolds. For a two-dimensional linear system and a common nonlinear test problem we show analytically that the variational approach asymptotically identifies the exact slow invariant manifold in the limit of both an infinite time horizon of the variational problem with fixed spectral gap and infinite spectral gap with a fixed finite time horizon. Numerical results for the linear and nonlinear model problems as well as a more realistic higher-dimensional chemical reaction mechanism are presented.

preprint2009arXiv

Minimal Curvature Trajectories: Riemannian Geometry Concepts for Model Reduction in Chemical Kinetics

In dissipative ordinary differential equation systems different time scales cause anisotropic phase volume contraction along solution trajectories. Model reduction methods exploit this for simplifying chemical kinetics via a time scale separation into fast and slow modes. The aim is to approximate the system dynamics with a dimension-reduced model after eliminating the fast modes by enslaving them to the slow ones via computation of a slow attracting manifold. We present a novel method for computing approximations of such manifolds using trajectory-based optimization. We discuss Riemannian geometry concepts as a basis for suitable optimization criteria characterizing trajectories near slow attracting manifolds and thus provide insight into fundamental geometric properties of multiple time scale chemical kinetics. The optimization criteria correspond to a suitable mathematical formulation of "minimal relaxation" of chemical forces along reaction trajectories under given constraints. We present various geometrically motivated criteria and the results of their application to three test case reaction mechanisms serving as examples. We demonstrate that accurate numerical approximations of slow invariant manifolds can be obtained.