Asymptotic approaches in inverse problems for depolymerization estimation

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Doumic, Marie, Moireau, Philippe
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866910626964570112
author Doumic, Marie
Moireau, Philippe
author_facet Doumic, Marie
Moireau, Philippe
contents Depolymerization reactions constitute frequent experiments, for instance in biochemistry for the study of amyloid fibrils. The quantities experimentally observed are related to the time dynamics of a quantity averaged over all polymer sizes, such as the total polymerised mass or the mean size of particles. The question analysed here is to link this measurement to the initial size distribution. To do so, we first derive, from the initial reaction system $$\mathcal{C}_i \stackrel{b}{\longrightarrow} \mathcal{C}_{i-1} + \mathcal{C}_1,\qquad i\geq i_0\geq 2,$$ two asymptotic models: at first order, a backward transport equation, and at second order, an advection-diffusion/Fokker-Planck equation complemented with a mixed boundary condition at $x = 0$. We estimate their distance to the original system solution. We then turn to the inverse problem, i.e., how to estimate the initial size distribution from the time measurement of an average quantity, given by a moment of the solution. This question has been already studied for the first order asymptotic model, and we analyse here the second order asymptotic. Thanks to Carleman inequalities and to log-convexity estimates, we prove observability results and error estimates for a Tikhonov regularization. We then develop a Kalman-based observer approach, and implement it on simulated observations. Despite its severely ill-posed character, the second-order approach appears numerically more accurate than the first-order one.
format Preprint
id arxiv_https___arxiv_org_abs_2410_00460
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Asymptotic approaches in inverse problems for depolymerization estimation
Doumic, Marie
Moireau, Philippe
Analysis of PDEs
Depolymerization reactions constitute frequent experiments, for instance in biochemistry for the study of amyloid fibrils. The quantities experimentally observed are related to the time dynamics of a quantity averaged over all polymer sizes, such as the total polymerised mass or the mean size of particles. The question analysed here is to link this measurement to the initial size distribution. To do so, we first derive, from the initial reaction system $$\mathcal{C}_i \stackrel{b}{\longrightarrow} \mathcal{C}_{i-1} + \mathcal{C}_1,\qquad i\geq i_0\geq 2,$$ two asymptotic models: at first order, a backward transport equation, and at second order, an advection-diffusion/Fokker-Planck equation complemented with a mixed boundary condition at $x = 0$. We estimate their distance to the original system solution. We then turn to the inverse problem, i.e., how to estimate the initial size distribution from the time measurement of an average quantity, given by a moment of the solution. This question has been already studied for the first order asymptotic model, and we analyse here the second order asymptotic. Thanks to Carleman inequalities and to log-convexity estimates, we prove observability results and error estimates for a Tikhonov regularization. We then develop a Kalman-based observer approach, and implement it on simulated observations. Despite its severely ill-posed character, the second-order approach appears numerically more accurate than the first-order one.
title Asymptotic approaches in inverse problems for depolymerization estimation
topic Analysis of PDEs
url https://arxiv.org/abs/2410.00460