@misc{AscherChinReich1994, author = {Ascher, Uri M. and Chin, Hongsheng and Reich, Sebastian}, title = {Stabilization of DAEs and invariant manifolds}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-15625}, year = {1994}, abstract = {Many methods have been proposed for the stabilization of higher index differential-algebraic equations (DAEs). Such methods often involve constraint differentiation and problem stabilization, thus obtaining a stabilized index reduction. A popular method is Baumgarte stabilization, but the choice of parameters to make it robust is unclear in practice. Here we explain why the Baumgarte method may run into trouble. We then show how to improve it. We further develop a unifying theory for stabilization methods which includes many of the various techniques proposed in the literature. Our approach is to (i) consider stabilization of ODEs with invariants, (ii) discretize the stabilizing term in a simple way, generally different from the ODE discretization, and (iii) use orthogonal projections whenever possible. The best methods thus obtained are related to methods of coordinate projection. We discuss them and make concrete algorithmic suggestions.}, language = {en} } @misc{Reich1995, author = {Reich, Sebastian}, title = {Smoothed dynamics of highly oscillatory Hamiltonian systems}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-15639}, year = {1995}, abstract = {We consider the numerical treatment of Hamiltonian systems that contain a potential which grows large when the system deviates from the equilibrium value of the potential. Such systems arise, e.g., in molecular dynamics simulations and the spatial discretization of Hamiltonian partial differential equations. Since the presence of highly oscillatory terms in the solutions forces any explicit integrator to use very small step size, the numerical integration of such systems provides a challenging task. It has been suggested before to replace the strong potential by a holonomic constraint that forces the solutions to stay at the equilibrium value of the potential. This approach has, e.g., been successfully applied to the bond stretching in molecular dynamics simulations. In other cases, such as the bond-angle bending, this methods fails due to the introduced rigidity. Here we give a careful analysis of the analytical problem by means of a smoothing operator. This will lead us to the notion of the smoothed dynamics of a highly oscillatory Hamiltonian system. Based on our analysis, we suggest a new constrained formulation that maintains the flexibility of the system while at the same time suppressing the high-frequency components in the solutions and thus allowing for larger time steps. The new constrained formulation is Hamiltonian and can be discretized by the well-known SHAKE method.}, language = {en} } @misc{LeimkuhlerReich1994, author = {Leimkuhler, Benedict and Reich, Sebastian}, title = {Symplectic integration of constrained Hamiltonian systems}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-15653}, year = {1994}, abstract = {A Hamiltonian system in potential form (formula in the original abstract) subject to smooth constraints on q can be viewed as a Hamiltonian system on a manifold, but numerical computations must be performed in Rn. In this paper methods which reduce "Hamiltonian differential algebraic equations" to ODEs in Euclidean space are examined. The authors study the construction of canonical parameterizations or local charts as well as methods based on the construction of ODE systems in the space in which the constraint manifold is embedded which preserve the constraint manifold as an invariant manifold. In each case, a Hamiltonian system of ordinary differential equations is produced. The stability of the constraint invariants and the behavior of the original Hamiltonian along solutions are investigated both numerically and analytically.}, language = {en} } @misc{AscherChinPetzoldetal.1994, author = {Ascher, Uri M. and Chin, Hongsheng and Petzold, Linda R. and Reich, Sebastian}, title = {Stabilization of constrained mechanical systems with DAEs and invariant manifolds}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-15698}, year = {1994}, abstract = {Many methods have been proposed for the simulation of constrained mechanical systems. The most obvious of these have mild instabilities and drift problems. Consequently, stabilization techniques have been proposed A popular stabilization method is Baumgarte's technique, but the choice of parameters to make it robust has been unclear in practice. Some of the simulation methods that have been proposed and used in computations are reviewed here, from a stability point of view. This involves concepts of differential-algebraic equation (DAE) and ordinary differential equation (ODE) invariants. An explanation of the difficulties that may be encountered using Baumgarte's method is given, and a discussion of why a further quest for better parameter values for this method will always remain frustrating is presented. It is then shown how Baumgarte's method can be improved. An efficient stabilization technique is proposed, which may employ explicit ODE solvers in case of nonstiff or highly oscillatory problems and which relates to coordinate projection methods. Examples of a two-link planar robotic arm and a squeezing mechanism illustrate the effectiveness of this new stabilization method.}, language = {en} } @misc{Reich1980, author = {Reich, Sebastian}, title = {Algebrodifferentialgleichungen und Vektorfelder auf Mannigfaltigkeiten}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-47290}, year = {1980}, abstract = {In diesem Beitrag wird der Zusammenhang zwischen Algebrodifferentialgleichungen (ADGL) und Vektorfeldern auf Mannigfaltigkeiten untersucht. Dazu wird zun{\"a}chst der Begriff der regul{\"a}ren ADGL eingef{\"u}hrt, wobei unter eirter regul{\"a}ren ADGL eine ADGL verstanden wird, deren L{\"o}sungsmenge identisch mit der L{\"o}sungsmenge eines Vektorfeldes ist. Ausgehend von bekannten Aussagen {\"u}ber die L{\"o}sungsmenge eines Vektorfeldes werden analoge Aussagen f{\"u}r die L{\"o}sungsmenge einer regul{\"a}ren ADGL abgeleitet. Es wird eine Reduktionsmethode angegeben, die auf ein Kriterium f{\"u}r die Begularit{\"a}t einer ADGL und auf die Definition des Index einer nichtlinearen ADGL f{\"u}hrt. Außerdem wird gezeigt, daß beliebige Vektorfelder durch regul{\"a}re ADGL so realisiert werden k{\"o}nnen, daß die L{\"o}sungsmenge des Vektorfeldes mit der der realisierenden ADGL identisch ist. Abschließend werden die f{\"u}r autonome ADGL gewonnenen Aussagen auf den Fall der nichtautonomen ADGL {\"u}bertragen.}, language = {de} } @misc{Reich1992, author = {Reich, Sebastian}, title = {Differential-algebraic equations and applications in circuit theory}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-46646}, year = {1992}, abstract = {Technical and physical systems, especially electronic circuits, are frequently modeled as a system of differential and nonlinear implicit equations. In the literature such systems of equations are called differentialalgebraic equations (DAEs). It turns out that the numerical and analytical properties of a DAE depend on an integer called the index of the problem. For example, the well-known BDF method of Gear can be applied, in general, to a DAE only if the index does not exceed one. In this paper we give a geometric interpretation of higherindex DAEs and indicate problems arising in connection with such DAEs by means of several examples.}, language = {en} } @misc{AsgarimehrWickertReich2019, author = {Asgarimehr, Milad and Wickert, Jens and Reich, Sebastian}, title = {Evaluating impact of rain attenuation on space-borne GNSS Reflectometry wind speeds}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1132}, issn = {1866-8372}, doi = {10.25932/publishup-47344}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-473441}, pages = {20}, year = {2019}, abstract = {The novel space-borne Global Navigation Satellite System Reflectometry (GNSS-R) technique has recently shown promise in monitoring the ocean state and surface wind speed with high spatial coverage and unprecedented sampling rate. The L-band signals of GNSS are structurally able to provide a higher quality of observations from areas covered by dense clouds and under intense precipitation, compared to those signals at higher frequencies from conventional ocean scatterometers. As a result, studying the inner core of cyclones and improvement of severe weather forecasting and cyclone tracking have turned into the main objectives of GNSS-R satellite missions such as Cyclone Global Navigation Satellite System (CYGNSS). Nevertheless, the rain attenuation impact on GNSS-R wind speed products is not yet well documented. Evaluating the rain attenuation effects on this technique is significant since a small change in the GNSS-R can potentially cause a considerable bias in the resultant wind products at intense wind speeds. Based on both empirical evidence and theory, wind speed is inversely proportional to derived bistatic radar cross section with a natural logarithmic relation, which introduces high condition numbers (similar to ill-posed conditions) at the inversions to high wind speeds. This paper presents an evaluation of the rain signal attenuation impact on the bistatic radar cross section and the derived wind speed. This study is conducted simulating GNSS-R delay-Doppler maps at different rain rates and reflection geometries, considering that an empirical data analysis at extreme wind intensities and rain rates is impossible due to the insufficient number of observations from these severe conditions. Finally, the study demonstrates that at a wind speed of 30 m/s and incidence angle of 30 degrees, rain at rates of 10, 15, and 20 mm/h might cause overestimation as large as approximate to 0.65 m/s (2\%), 1.00 m/s (3\%), and 1.3 m/s (4\%), respectively, which are still smaller than the CYGNSS required uncertainty threshold. The simulations are conducted in a pessimistic condition (severe continuous rainfall below the freezing height and over the entire glistening zone) and the bias is expected to be smaller in size in real environments.}, language = {en} } @misc{NueskenReichRozdeba2019, author = {N{\"u}sken, Nikolas and Reich, Sebastian and Rozdeba, Paul J.}, title = {State and parameter estimation from observed signal increments}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {916}, issn = {1866-8372}, doi = {10.25932/publishup-44260}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-442609}, pages = {25}, year = {2019}, abstract = {The success of the ensemble Kalman filter has triggered a strong interest in expanding its scope beyond classical state estimation problems. In this paper, we focus on continuous-time data assimilation where the model and measurement errors are correlated and both states and parameters need to be identified. Such scenarios arise from noisy and partial observations of Lagrangian particles which move under a stochastic velocity field involving unknown parameters. We take an appropriate class of McKean-Vlasov equations as the starting point to derive ensemble Kalman-Bucy filter algorithms for combined state and parameter estimation. We demonstrate their performance through a series of increasingly complex multi-scale model systems.}, language = {en} } @misc{AcevedoReichCubasch2015, author = {Acevedo, Walter and Reich, Sebastian and Cubasch, Ulrich}, title = {Towards the assimilation of tree-ring-width records using ensemble Kalman filtering techniques}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, volume = {46}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, number = {892}, issn = {1866-8372}, doi = {10.25932/publishup-43636}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-436363}, pages = {1909 -- 1920}, year = {2015}, abstract = {This paper investigates the applicability of the Vaganov-Shashkin-Lite (VSL) forward model for tree-ring-width chronologies as observation operator within a proxy data assimilation (DA) setting. Based on the principle of limiting factors, VSL combines temperature and moisture time series in a nonlinear fashion to obtain simulated TRW chronologies. When used as observation operator, this modelling approach implies three compounding, challenging features: (1) time averaging, (2) "switching recording" of 2 variables and (3) bounded response windows leading to "thresholded response". We generate pseudo-TRW observations from a chaotic 2-scale dynamical system, used as a cartoon of the atmosphere-land system, and attempt to assimilate them via ensemble Kalman filtering techniques. Results within our simplified setting reveal that VSL's nonlinearities may lead to considerable loss of assimilation skill, as compared to the utilization of a time-averaged (TA) linear observation operator. In order to understand this undesired effect, we embed VSL's formulation into the framework of fuzzy logic (FL) theory, which thereby exposes multiple representations of the principle of limiting factors. DA experiments employing three alternative growth rate functions disclose a strong link between the lack of smoothness of the growth rate function and the loss of optimality in the estimate of the TA state. Accordingly, VSL's performance as observation operator can be enhanced by resorting to smoother FL representations of the principle of limiting factors. This finding fosters new interpretations of tree-ring-growth limitation processes.}, language = {en} } @misc{AcevedoFallahReichetal.2017, author = {Acevedo, Walter and Fallah, Bijan and Reich, Sebastian and Cubasch, Ulrich}, title = {Assimilation of pseudo-tree-ring-width observations into an atmospheric general circulation model}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, number = {627}, issn = {1866-8372}, doi = {10.25932/publishup-41874}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-418743}, pages = {545 -- 557}, year = {2017}, abstract = {Paleoclimate data assimilation (DA) is a promising technique to systematically combine the information from climate model simulations and proxy records. Here, we investigate the assimilation of tree-ring-width (TRW) chronologies into an atmospheric global climate model using ensemble Kalman filter (EnKF) techniques and a process-based tree-growth forward model as an observation operator. Our results, within a perfect-model experiment setting, indicate that the "online DA" approach did not outperform the "off-line" one, despite its considerable additional implementation complexity. On the other hand, it was observed that the nonlinear response of tree growth to surface temperature and soil moisture does deteriorate the operation of the time-averaged EnKF methodology. Moreover, for the first time we show that this skill loss appears significantly sensitive to the structure of the growth rate function, used to represent the principle of limiting factors (PLF) within the forward model. In general, our experiments showed that the error reduction achieved by assimilating pseudo-TRW chronologies is modulated by the magnitude of the yearly internal variability in themodel. This result might help the dendrochronology community to optimize their sampling efforts.}, language = {en} } @misc{Reich1990, author = {Reich, Sebastian}, title = {On a geometrical interpretation of differential-algebraic equations}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-46683}, year = {1990}, abstract = {The subject of this paper is the relation of differential-algebraic equations (DAEs) to vector fields on manifolds. For that reason, we introduce the notion of a regular DAE as a DAE to which a vector field uniquely corresponds. Furthermore, a technique is described which yields a family of manifolds for a given DAE. This socalled family of constraint manifolds allows in turn the formulation of sufficient conditions for the regularity of a DAE. and the definition of the index of a regular DAE. We also state a method for the reduction of higher-index DAEs to lowsr-index ones that can be solved without introducing additional constants of integration. Finally, the notion of realizability of a given vector field by a regular DAE is introduced, and it is shown that any vector field can be realized by a regular DAE. Throughout this paper the problem of path-tracing is discussed as an illustration of the mathematical phenomena.}, language = {en} } @misc{Reich1994, author = {Reich, Sebastian}, title = {Momentum conserving symplectic integrators}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-16824}, year = {1994}, abstract = {In this paper, we show that symplectic partitioned Runge-Kutta methods conserve momentum maps corresponding to linear symmetry groups acting on the phase space of Hamiltonian differential equations by extended point transformation. We also generalize this result to constrained systems and show how this conservation property relates to the symplectic integration of Lie-Poisson systems on certain submanifolds of the general matrix group GL(n).}, language = {en} } @misc{Reich1995, author = {Reich, Sebastian}, title = {On the local qualitative behavior of differential-algebraic equations}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-46739}, year = {1995}, abstract = {A theoretical famework for the investigation of the qualitative behavior of differential-algebraic equations (DAEs) near an equilibrium point is established. The key notion of our approach is the notion of regularity. A DAE is called regular locally around an equilibrium point if there is a unique vector field such that the solutions of the DAE and the vector field are in one-to-one correspondence in a neighborhood of this equili Drium point. Sufficient conditions for the regularity of an equilibrium point are stated. This in turn allows us to translate several local results, as formulated for vector fields, to DAEs that are regular locally around a g: ven equilibrium point (e.g. Local Stable and Unstable Manifold Theorem, Hopf theorem). It is important that ihese theorems are stated in terms of the given problem and not in terms of the corresponding vector field.}, language = {en} } @misc{Reich1991, author = {Reich, Sebastian}, title = {On an existence and uniqueness theory for nonlinear differential-algebraic equations}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-46706}, year = {1991}, abstract = {An existence and uniqueness theory is developed for general nonlinear and nonautonomous differential-algebraic equations (DAEs) by exploiting their underlying differential-geometric structure. A DAE is called regular if there is a unique nonautonomous vector field such that the solutions of the DAE and the solutions of the vector field are in one-to-one correspondence. Sufficient conditions for regularity of a DAE are derived in terms of constrained manifolds. Based on this differential-geometric characterization, existence and uniqueness results are stated for regular DAEs. Furthermore, our not ons are compared with techniques frequently used in the literature such as index and solvability. The results are illustrated in detail by means of a simple circuit example.}, language = {en} }