Constructive Methods of Invariant Manifolds for Kinetic Problems

We present the Constructive Methods of Invariant Manifolds for model reduction in physical and chemical kinetics, developed during last two decades. The physical problem of reduced description is studied in a most general form as a problem of constructing the slow invariant manifold. The invariance conditions are formulated as the differential equation for a manifold immersed in the phase space (the invariance equation). The equation of motion for immersed manifolds is obtained (the film extension of the dynamics). Invariant manifolds are fixed points for this equation, and slow invariant manifolds are Lyapunov stable fixed points, thus slowness is presented as stability. A collection of methods for construction of slow invariant manifolds is presented, in particular, the Newton method subject to incomplete linearization is the analogue of KAM methods for dissipative systems. The systematic use of thermodynamics structures and of the quasi--chemical representation allow to construct approximations which are in concordance with physical restrictions. We systematically consider a discrete analogue of the slow (stable) positively invariant manifolds for dissipative systems, invariant grids. Dynamic and static postprocessing procedures give us the opportunity to estimate the accuracy of obtained approximations, and to improve this accuracy significantly. The following examples of applications are presented: Nonperturbative deviation of physically consistent hydrodynamics from the Boltzmann equation and from the reversible dynamics, for Knudsen numbers Kn~1; construction of the moment equations for nonequilibrium media and their dynamical correction (instead of extension of list of variables) to gain more accuracy in description of highly nonequilibrium flows; determination of molecules dimension (as diameters of equivalent hard spheres) from experimental viscosity data; invariant grids for a two-dimensional catalytic reaction and a four-dimensional oxidation reaction (six species, two balances); universal continuous media description of dilute polymeric solution; the limits of macroscopic description for polymer molecules, etc

Spectral Methods for Numerical Relativity

Version published online by Living Reviews in Relativity.International audienceEquations arising in General Relativity are usually too complicated to be solved analytically and one has to rely on numerical methods to solve sets of coupled partial differential equations. Among the possible choices, this paper focuses on a class called spectral methods where, typically, the various functions are expanded onto sets of orthogonal polynomials or functions. A theoretical introduction on spectral expansion is first given and a particular emphasis is put on the fast convergence of the spectral approximation. We present then different approaches to solve partial differential equations, first limiting ourselves to the one-dimensional case, with one or several domains. Generalization to more dimensions is then discussed. In particular, the case of time evolutions is carefully studied and the stability of such evolutions investigated. One then turns to results obtained by various groups in the field of General Relativity by means of spectral methods. First, works which do not involve explicit time-evolutions are discussed, going from rapidly rotating strange stars to the computation of binary black holes initial data. Finally, the evolutions of various systems of astrophysical interest are presented, from supernovae core collapse to binary black hole mergers

A Surrogate Model of Gravitational Waveforms from Numerical Relativity Simulations of Precessing Binary Black Hole Mergers

We present the first surrogate model for gravitational waveforms from the coalescence of precessing binary black holes. We call this surrogate model NRSur4d2s. Our methodology significantly extends recently introduced reduced-order and surrogate modeling techniques, and is capable of directly modeling numerical relativity waveforms without introducing phenomenological assumptions or approximations to general relativity. Motivated by GW150914, LIGO's first detection of gravitational waves from merging black holes, the model is built from a set of $276$ numerical relativity (NR) simulations with mass ratios $q \leq 2$, dimensionless spin magnitudes up to $0.8$, and the restriction that the initial spin of the smaller black hole lies along the axis of orbital angular momentum. It produces waveforms which begin $\sim 30$ gravitational wave cycles before merger and continue through ringdown, and which contain the effects of precession as well as all $\ell \in \{2, 3\}$ spin-weighted spherical-harmonic modes. We perform cross-validation studies to compare the model to NR waveforms \emph{not} used to build the model, and find a better agreement within the parameter range of the model than other, state-of-the-art precessing waveform models, with typical mismatches of $10^{-3}$. We also construct a frequency domain surrogate model (called NRSur4d2s_FDROM) which can be evaluated in $50\, \mathrm{ms}$ and is suitable for performing parameter estimation studies on gravitational wave detections similar to GW150914.Comment: 34 pages, 26 figure

About regularity properties in variational analysis and applications in optimization

Regularity properties lie at the core of variational analysis because of their importance for stability analysis of optimization and variational problems, constraint qualications, qualication conditions in coderivative and subdierential calculus and convergence analysis of numerical algorithms. The thesis is devoted to investigation of several research questions related to regularity properties in variational analysis and their applications in convergence analysis and optimization. Following the works by Kruger, we examine several useful regularity properties of collections of sets in both linear and Holder-type settings and establish their characterizations and relationships to regularity properties of set-valued mappings. Following the recent publications by Lewis, Luke, Malick (2009), Drusvyatskiy, Ioe, Lewis (2014) and some others, we study application of the uniform regularity and related properties of collections of sets to alternating projections for solving nonconvex feasibility problems and compare existing results on this topic. Motivated by Ioe (2000) and his subsequent publications, we use the classical iteration scheme going back to Banach, Schauder, Lyusternik and Graves to establish criteria for regularity properties of set-valued mappings and compare this approach with the one based on the Ekeland variational principle. Finally, following the recent works by Khanh et al. on stability analysis for optimization related problems, we investigate calmness of set-valued solution mappings of variational problems.Doctor of Philosoph

Compact binary merger simulations in numerical relativity

The era of Gravitational Waves Astronomy was launched after the success of the first observation run of the LIGO Scientific Collaboration and the VIRGO Collaboration. The large laser interferometers incredible achievement prompted the need of extensive studies in the field of compact astrophysical objects, such as Black Holes and Neutron Stars. Today, seven years after this event, the field of study underwent a notable expansion, corroborated by the detection of a signal coming from a Binary Neutron Star merger, together with its electro-magnetic counterpart, and, more recently, by the first detections of signals coming from mixed compact binaries, i.e. Black Hole- Neutron Star binaries. In this thesis work we span our attention across different aspects of compact objects mergers, including the inclusion of new physics into the already performing numerical relativity code BAM and the study of specific systems of compact objects. We first explore the treatment of neutrinos in case of Binary Neutron Star mergers and a tool to identify and further analyze regions containing trapped neutrinos, in the hot remnant of such mergers. Neutrinos, play in fact a key role into the rapid processes that characterize the formation of elements in the dynamical ejecta expelled during these catastrophic events. In the following we explore a variety of configurations of mixed compact binary systems. After the development of the new ID code Elliptica, and the steps taken to verify its accuracy, we make use of its capability to evolve sets of physical system with various properties. Exploring the space of parameters we study different spin configurations and magnitudes of single objects and their effects on the merger dynamics