8,005 research outputs found
Local Ensemble Transform Kalman Filter: a non-stationary control law for complex adaptive optics systems on ELTs
We propose a new algorithm for an adaptive optics system control law which
allows to reduce the computational burden in the case of an Extremely Large
Telescope (ELT) and to deal with non-stationary behaviors of the turbulence.
This approach, using Ensemble Transform Kalman Filter and localizations by
domain decomposition is called the local ETKF: the pupil of the telescope is
split up into various local domains and calculations for the update estimate of
the turbulent phase on each domain are performed independently. This data
assimilation scheme enables parallel computation of markedly less data during
this update step. This adapts the Kalman Filter to large scale systems with a
non-stationary turbulence model when the explicit storage and manipulation of
extremely large covariance matrices are impossible. First simulation results
are given in order to assess the theoretical analysis and to demonstrate the
potentiality of this new control law for complex adaptive optics systems on
ELTs.Comment: Proceedings of the AO4ELT3 conference; 8 pages, 3 figure
Parallel framework for dynamic domain decomposition of data assimilation problems: a case study on Kalman Filter algorithm
We focus on Partial Differential Equation (PDE)âbased Data Assimilation problems (DA) solved by means of variational approaches and Kalman filter algorithm. Recently, we presented a Domain Decomposition framework (we call it DDâDA, for short) performing a decomposition of the whole physical domain along space and time directions, and joining the idea of Schwarz's methods and parallel in time approaches. For effective parallelization of DDâDA algorithms, the computational load assigned to subdomains must be equally distributed. Usually computational cost is proportional to the amount of data entities assigned to partitions. Good quality partitioning also requires the volume of communication during calculation to be kept at its minimum. In order to deal with DDâDA problems where the observations are nonuniformly distributed and general sparse, in the present work we employ a parallel load balancing algorithm based on adaptive and dynamic defining of boundaries of DDâwhich is aimed to balance workload according to data location. We call it DyDD. As the numerical model underlying DA problems arising from the soâcalled discretizeâthenâoptimize approach is the constrained least square model (CLS), we will use CLS as a reference state estimation problem and we validate DyDD on different scenario
Local ensemble transform Kalman filter, a fast non-stationary control law for adaptive optics on ELTs: theoretical aspects and first simulation results
We propose a new algorithm for an adaptive optics system control law, based
on the Linear Quadratic Gaussian approach and a Kalman Filter adaptation with
localizations. It allows to handle non-stationary behaviors, to obtain
performance close to the optimality defined with the residual phase variance
minimization criterion, and to reduce the computational burden with an
intrinsically parallel implementation on the Extremely Large Telescopes (ELTs).Comment: This paper was published in Optics Express and is made available as
an electronic reprint with the permission of OSA. The paper can be found at
the following URL on the OSA website: http://www.opticsinfobase.org/oe/ .
Systematic or multiple reproduction or distribution to multiple locations via
electronic or other means is prohibited and is subject to penalties under la
Reduced order modeling of fluid flows: Machine learning, Kolmogorov barrier, closure modeling, and partitioning
In this paper, we put forth a long short-term memory (LSTM) nudging framework
for the enhancement of reduced order models (ROMs) of fluid flows utilizing
noisy measurements. We build on the fact that in a realistic application, there
are uncertainties in initial conditions, boundary conditions, model parameters,
and/or field measurements. Moreover, conventional nonlinear ROMs based on
Galerkin projection (GROMs) suffer from imperfection and solution instabilities
due to the modal truncation, especially for advection-dominated flows with slow
decay in the Kolmogorov width. In the presented LSTM-Nudge approach, we fuse
forecasts from a combination of imperfect GROM and uncertain state estimates,
with sparse Eulerian sensor measurements to provide more reliable predictions
in a dynamical data assimilation framework. We illustrate the idea with the
viscous Burgers problem, as a benchmark test bed with quadratic nonlinearity
and Laplacian dissipation. We investigate the effects of measurements noise and
state estimate uncertainty on the performance of the LSTM-Nudge behavior. We
also demonstrate that it can sufficiently handle different levels of temporal
and spatial measurement sparsity. This first step in our assessment of the
proposed model shows that the LSTM nudging could represent a viable realtime
predictive tool in emerging digital twin systems
4D Seismic History Matching Incorporating Unsupervised Learning
The work discussed and presented in this paper focuses on the history
matching of reservoirs by integrating 4D seismic data into the inversion
process using machine learning techniques. A new integrated scheme for the
reconstruction of petrophysical properties with a modified Ensemble Smoother
with Multiple Data Assimilation (ES-MDA) in a synthetic reservoir is proposed.
The permeability field inside the reservoir is parametrised with an
unsupervised learning approach, namely K-means with Singular Value
Decomposition (K-SVD). This is combined with the Orthogonal Matching Pursuit
(OMP) technique which is very typical for sparsity promoting regularisation
schemes. Moreover, seismic attributes, in particular, acoustic impedance, are
parametrised with the Discrete Cosine Transform (DCT). This novel combination
of techniques from machine learning, sparsity regularisation, seismic imaging
and history matching aims to address the ill-posedness of the inversion of
historical production data efficiently using ES-MDA. In the numerical
experiments provided, I demonstrate that these sparse representations of the
petrophysical properties and the seismic attributes enables to obtain better
production data matches to the true production data and to quantify the
propagating waterfront better compared to more traditional methods that do not
use comparable parametrisation techniques
Identification of weakly coupled multiphysics problems. Application to the inverse problem of electrocardiography
This work addresses the inverse problem of electrocardiography from a new
perspective, by combining electrical and mechanical measurements. Our strategy
relies on the defini-tion of a model of the electromechanical contraction which
is registered on ECG data but also on measured mechanical displacements of the
heart tissue typically extracted from medical images. In this respect, we
establish in this work the convergence of a sequential estimator which combines
for such coupled problems various state of the art sequential data assimilation
methods in a unified consistent and efficient framework. Indeed we ag-gregate a
Luenberger observer for the mechanical state and a Reduced Order Unscented
Kalman Filter applied on the parameters to be identified and a POD projection
of the electrical state. Then using synthetic data we show the benefits of our
approach for the estimation of the electrical state of the ventricles along the
heart beat compared with more classical strategies which only consider an
electrophysiological model with ECG measurements. Our numerical results
actually show that the mechanical measurements improve the identifiability of
the electrical problem allowing to reconstruct the electrical state of the
coupled system more precisely. Therefore, this work is intended to be a first
proof of concept, with theoretical justifications and numerical investigations,
of the ad-vantage of using available multi-modal observations for the
estimation and identification of an electromechanical model of the heart
- âŠ