199 research outputs found
A dynamical adaptive tensor method for the Vlasov-Poisson system
A numerical method is proposed to solve the full-Eulerian time-dependent
Vlasov-Poisson system in high dimension. The algorithm relies on the
construction of a tensor decomposition of the solution whose rank is adapted at
each time step. This decomposition is obtained through the use of an efficient
modified Progressive Generalized Decomposition (PGD) method, whose convergence
is proved. We suggest in addition a symplectic time-discretization splitting
scheme that preserves the Hamiltonian properties of the system. This scheme is
naturally obtained by considering the tensor structure of the approximation.
The efficiency of our approach is illustrated through time-dependent 2D-2D
numerical examples
Approximated Lax Pairs for the Reduced Order Integration of Nonlinear Evolution Equations
A reduced-order model algorithm, called ALP, is proposed to solve nonlinear
evolution partial differential equations. It is based on approximations of
generalized Lax pairs. Contrary to other reduced-order methods, like Proper
Orthogonal Decomposition, the basis on which the solution is searched for
evolves in time according to a dynamics specific to the problem. It is
therefore well-suited to solving problems with progressive front or wave
propagation. Another difference with other reduced-order methods is that it is
not based on an off-line / on-line strategy. Numerical examples are shown for
the linear advection, KdV and FKPP equations, in one and two dimensions
Reduced-Order Modeling based on Approximated Lax Pairs
A reduced-order model algorithm, based on approximations of Lax pairs, is
proposed to solve nonlinear evolution partial differential equations. Contrary
to other reduced-order methods, like Proper Orthogonal Decomposition, the space
where the solution is searched for evolves according to a dynamics specific to
the problem. It is therefore well-suited to solving problems with progressive
waves or front propagation. Numerical examples are shown for the KdV and FKPP
(nonlinear reaction diffusion) equations, in one and two dimensions
Fast reconstruction of 3D blood flows from Doppler ultrasound images and reduced models
This paper deals with the problem of building fast and reliable 3D
reconstruction methods for blood flows for which partial information is given
by Doppler ultrasound measurements. This task is of interest in medicine since
it could enrich the available information used in the diagnosis of certain
diseases which is currently based essentially on the measurements coming from
ultrasound devices. The fast reconstruction of the full flow can be performed
with state estimation methods that have been introduced in recent years and
that involve reduced order models. One simple and efficient strategy is the
so-called Parametrized Background Data-Weak approach (PBDW). It is a linear
mapping that consists in a least squares fit between the measurement data and a
linear reduced model to which a certain correction term is added. However, in
the original approach, the reduced model is built a priori and independently of
the reconstruction task (typically with a proper orthogonal decomposition or a
greedy algorithm). In this paper, we investigate the construction of other
reduced spaces which are built to be better adapted to the reconstruction task
and which result in mappings that are sometimes nonlinear. We compare the
performance of the different algorithms on numerical experiments involving
synthetic Doppler measurements. The results illustrate the superiority of the
proposed alternatives to the classical linear PBDW approach
Eulerian models and algorithms for unbalanced optimal transport
International audienceBenamou and Brenier formulation of Monge transportation problem (Numer. Math. 84:375-393, 2000) has proven to be of great interest in image processing to compute warpings and distances between pair of images (SIAM J. Math. Analysis, 35:61-97, 2003). One requirement for the algorithm to work is to interpolate densities of same mass. In most applications to image interpolation, this is a serious limitation. Existing approaches to overcome this caveat are reviewed, and discussed. Due to the mix between transport and interpolation, these models can produce instantaneous motion at finite range. In this paper we propose new methods, parameter-free, for interpolating unbalanced densities. One of our motivations is the application to interpolation of growing tumor images
State estimation in nonlinear parametric time dependent systems using Tensor Train
International audienceIn the present work we propose a reduced-order method to solve the state estimation problem when nonlinear parametric time-dependent systems are at hand. The method is based on the approximation of the set of system solutions by means of a Tensor Train format. The particular structure of Tensor Train makes it possible to set up both a variational and a sequential method. Several numerical experiments are proposed to assess the behaviour of the method
State Estimation with Model Reduction and Shape Variability. Application to biomedical problems
We develop a mathematical and numerical framework to solve state estimation
problems for applications that present variations in the shape of the spatial
domain. This situation arises typically in a biomedical context where inverse
problems are posed on certain organs or portions of the body which inevitably
involve morphological variations. If one wants to provide fast reconstruction
methods, the algorithms must take into account the geometric variability. We
develop and analyze a method which allows to take this variability into account
without needing any a priori knowledge on a parametrization of the geometrical
variations. For this, we rely on morphometric techniques involving
Multidimensional Scaling, and couple them with reconstruction algorithms that
make use of reduced model spaces pre-computed on a database of geometries. We
prove the potential of the method on a synthetic test problem inspired from the
reconstruction of blood flows and quantities of medical interest with Doppler
ultrasound imaging
Reduced order modelling for direct and inverse problems in haemodynamics
International audienceIn this chapter we propose a review of reduced-order models (ROMs) to speed-up direct and inverse problems in the context of haemodynamics. In particular, we highlight three different ways of building them and comment upon the numerous contributions and methodological advancements in this field
An Information-Theoretic Framework for Optimal Design: Analysis of Protocols for Estimating Soft Tissue Parameters in Biaxial Experiments
A new framework for optimal design based on the information-theoretic measures of mutual information, conditional mutual information and their combination is proposed. The framework is tested on the analysis of protocolsâa combination of angles along which strain measurements can be acquiredâin a biaxial experiment of soft tissues for the estimation of hyperelastic constitutive model parameters. The proposed framework considers the information gain about the parameters from the experiment as the key criterion to be maximised, which can be directly used for optimal design. Information gain is computed through k-nearest neighbour algorithms applied to the joint samples of the parameters and measurements produced by the forward and observation models. For biaxial experiments, the results show that low angles have a relatively low information content compared to high angles. The results also show that a smaller number of angles with suitably chosen combinations can result in higher information gains when compared to a larger number of angles which are poorly combined. Finally, it is shown that the proposed framework is consistent with classical approaches, particularly D-optimal design
- âŠ