Search

A flexible model for non-stationary Gaussian random fields on hypersurfaces is introduced. The class of random fields on curves and surfaces is characterized by a power spectral density of a second order elliptic differential operator. Sampling is done by a Galerkin--Chebyshev approximation based on the surface finite element method and Chebyshev polynomials. Strong error bounds are shown with convergence rates depending on the smoothness of the approximated random field., Comment: 32 pages, 4 figures

Incremental stability is assessed for the Unidirectional Flow Model, which describes the flow of a conserved quantity in a cascade of compartments. The model is governed by a cooperative system of nonlinear ordinary differential equations with a right-hand side that depends on an input function. Suitable conditions are given on the input function such that the solutions to the model converge to each other at an exponential rate. Milder conditions, compared with the latter, are given such that the solutions to the model converge to each other asymptotically. These conditions are related to inflow-connected and outflow-connected compartmental systems and are easy to check. Computational means to estimate the exponential rate are given and tested against numerical solutions of the model. A natural application of incremental stability is state estimation, demonstrated in the numerical experiments, put into a traffic flow context.

Control of mixed-autonomy traffic where Human-driven Vehicles (HVs) and Autonomous Vehicles (AVs) coexist on the road has gained increasing attention over the recent decades. This paper addresses the boundary stabilization problem for mixed traffic on freeways. The traffic dynamics are described by uncertain coupled hyperbolic partial differential equations (PDEs) with Markov jumping parameters, which aim to address the distinctive driving strategies between AVs and HVs. Considering that the spacing policies of AVs vary in mixed traffic, the stochastic impact area of AVs is governed by a continuous Markov chain. The interactions between HVs and AVs such as overtaking or lane changing are mainly induced by impact areas. Using backstepping design, we develop a full-state feedback boundary control law to stabilize the deterministic system (nominal system). Applying Lyapunov analysis, we demonstrate that the nominal backstepping control law is able to stabilize the traffic system with Markov jumping parameters, provided the nominal parameters are sufficiently close to the stochastic ones on average. The mean-square exponential stability conditions are derived, and the results are validated by numerical simulations.

Large or very large spatial (and spatio-temporal) datasets have become common place in many environmental and climate studies. These data are often collected in non-Euclidean spaces (such as the planet Earth) and they often present non-stationary anisotropies. This paper proposes a generic approach to model Gaussian Random Fields (GRFs) on compact Riemannian manifolds that bridges the gap between existing works on non-stationary GRFs and random fields on manifolds. This approach can be applied to any smooth compact manifolds, and in particular to any compact surface. By defining a Riemannian metric that accounts for the preferential directions of correlation, our approach yields an interpretation of the ''local anisotropies'' as resulting from ''local'' deformations of the domain. We provide scalable algorithms for the estimation of the parameters and for optimal prediction by kriging and simulation able to tackle very large grids. Stationary and non-stationary illustrations are provided.

In this work, we propose an algorithm performing short-term predictions of the flux of vehicles on a stretch of road, using past measurements of the flux. This algorithm is based on a physics-aware recurrent neural network. A discretization of a macroscopic traffic flow model (using the so-called Traffic Reaction Model) is embedded in the architecture of the network and yields flux predictions based on estimated and predicted space-time dependent traffic parameters. These parameters are themselves obtained using a succession of LSTM ans simple recurrent neural networks. Besides, on top of the predictions, the algorithm yields a smoothing of its inputs which is also physically-constrained by the macroscopic traffic flow model. The algorithm is tested on raw flux measurements obtained from loop detectors., Comment: 17 pages, 11 figures, 2 tables

A new numerical approximation method for a class of Gaussian random fields on compact connected oriented Riemannian manifolds is introduced. This class of random fields is characterized by the Laplace--Beltrami operator on the manifold. A Galerkin approximation is combined with a polynomial approximation using Chebyshev series. This so-called Galerkin--Chebyshev approximation scheme yields efficient and generic sampling algorithms for Gaussian random fields on manifolds. Strong and weak orders of convergence for the Galerkin approximation and strong convergence orders for the Galerkin--Chebyshev approximation are shown and confirmed through numerical experiments., Comment: Version submitted to a peer-reviewed journal. Changes: fixed residual typos, new outline. 33 pages, 5 figures

The main motivation of this work is to assess the validity of a LWR traffic flow model to model measurements obtained from trajectory data, and propose extensions of this model to improve it. A formulation for a discrete dynamical system is proposed aiming at reproducing the evolution in time of the density of vehicles along a road, as observed in the measurements. This system is formulated as a chemical reaction network where road cells are interpreted as compartments, the transfer of vehicles from one cell to the other is seen as a chemical reaction between adjacent compartment and the density of vehicles is seen as a concentration of reactant. Several degrees of flexibility on the parameters of this system, which basically consist of the reaction rates between the compartments, can be considered: a constant value or a function depending on time and/or space. Density measurements coming from trajectory data are then interpreted as observations of the states of this system at consecutive times. Optimal reaction rates for the system are then obtained by minimizing the discrepancy between the output of the system and the state measurements. This approach was tested both on simulated and real data, proved successful in recreating the complexity of traffic flows despite the assumptions on the flux-density relation., Comment: 31 pages, 18 figures, 4 tables

In this paper, we present a novel approach to geostatistical filtering which tackles two challenges encountered when applying this method to complex spatial datasets: modeling the non-stationarity of the data while still being able to work with large datasets. The approach is based on a finite element approximation of Gaussian random fields expressed as an expansion of the eigenfunctions of a Laplace--Beltrami operator defined to account for local anisotropies. The numerical approximation of the resulting random fields using a finite element approach is then leveraged to solve the scalability issue through a matrix-free approach. Finally, two cases of application of this approach, on simulated and real seismic data are presented., Comment: 25 pages, 8 figures

In this paper, we present finite element approximations of a class of Generalized random fields defined over a bounded domain of R d or a smooth d-dimensional Riemannian manifold (d $\ge$ 1). An explicit expression for the covariance matrix of the weights of the finite element representation of these fields is provided and an analysis of the approximation error is carried out. Finally, a method to generate simulations of these weights while limiting computational and storage costs is presented.

This paper presents an algorithm to simulate Gaussian random vectors whose precision matrix can be expressed as a polynomial of a sparse matrix. This situation arises in particular when simulating Gaussian Markov random fields obtained by the discretization by finite elements of the solutions of some stochastic partial derivative equations. The proposed algorithm uses a Chebyshev polynomial approximation to compute simulated vectors with a linear complexity. This method is asymptotically exact as the approximation order grows. Criteria based on tests of the statistical properties of the produced vectors are derived to determine minimal orders of approximation., Comment: 20 pages, 5 figures

An experimental study was conducted at the Ballistic Laboratory of NASA Glenn Research Center to study the impact of ice particles on a stationary flat surface target set at 45 degrees with respect to the direction of motion of the impinging particle (Figure 1). The experiment is part of NASA efforts to study the physics involved in engine power-loss events due to ice-crystal ingestion and ice accretion formation inside engines. These events can occur when aircraft encounter high-altitude convective weather.

Seven cast and one wrought gamma-alloys were ballistically impacted and tested under high cycle fatigue. The fatigue strength of each alloy was characterized as a function of initial flaw size and modeled using a threshold-based fracture mechanics approach.