首页 > 最新文献

Inverse Problems in Science and Engineering最新文献

英文 中文
Analysis for two-dimensional inverse quasilinear parabolic problem by Fourier method 二维拟线性逆抛物问题的傅里叶分析
IF 1.3 4区 工程技术 Q3 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2021-02-25 DOI: 10.1080/17415977.2021.1890068
F. Kanca, İ. Bağlan
In this work, two-dimensional inverse quasi-linear parabolic problem with periodic boundary and integral overdetermination conditions is investigated. The formal solution is obtained by the Fourier approximation. Under some natural regularity and consistency conditions on the input data,the existence, uniqueness and continuously dependence upon the data of the solution are proved by iteration method. The inverse problem is first examined by linearization and then used implicit finite difference scheme for the numerical solution. Also predictor corrector method is considered in the numerical approach. Some results on the numerical solution with two examples are presented with figures and tables. The sensitivity of the scheme with respect to noisy overdetermination data is illustrated.
本文研究了具有周期边界和积分超定条件的二维拟线性抛物型反问题。形式解是通过傅立叶近似得到的。在输入数据上的一些自然规律性和一致性条件下,用迭代方法证明了解的存在性、唯一性和对数据的连续依赖性。逆问题首先通过线性化进行检验,然后使用隐式有限差分格式进行数值求解。数值方法中还考虑了预测-校正方法。通过两个实例给出了数值解的一些结果。说明了该方案对噪声过定数据的敏感性。
{"title":"Analysis for two-dimensional inverse quasilinear parabolic problem by Fourier method","authors":"F. Kanca, İ. Bağlan","doi":"10.1080/17415977.2021.1890068","DOIUrl":"https://doi.org/10.1080/17415977.2021.1890068","url":null,"abstract":"In this work, two-dimensional inverse quasi-linear parabolic problem with periodic boundary and integral overdetermination conditions is investigated. The formal solution is obtained by the Fourier approximation. Under some natural regularity and consistency conditions on the input data,the existence, uniqueness and continuously dependence upon the data of the solution are proved by iteration method. The inverse problem is first examined by linearization and then used implicit finite difference scheme for the numerical solution. Also predictor corrector method is considered in the numerical approach. Some results on the numerical solution with two examples are presented with figures and tables. The sensitivity of the scheme with respect to noisy overdetermination data is illustrated.","PeriodicalId":54926,"journal":{"name":"Inverse Problems in Science and Engineering","volume":"29 1","pages":"1912 - 1945"},"PeriodicalIF":1.3,"publicationDate":"2021-02-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1080/17415977.2021.1890068","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43649244","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Bayesian estimation and uncertainty quantification in models of urea hydrolysis by E. coli biofilms 大肠杆菌生物膜水解尿素模型的贝叶斯估计和不确定度量化
IF 1.3 4区 工程技术 Q3 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2021-02-24 DOI: 10.1080/17415977.2021.1887172
B. Jackson, J. Connolly, R. Gerlach, I. Klapper, A. Parker
ABSTRACT Urea-hydrolysing biofilms are crucial to applications in medicine, engineering, and science. Quantitative information about ureolysis rates in biofilms is required to model these applications. We formulate a novel model of urea consumption in a biofilm that allows different kinetics, for example either first order or Michaelis–Menten. The model is fit to synthetic data to validate and compare two approaches, Bayesian and nonlinear least squares (NLS), commonly used by biofilm practitioners. The shortcomings of NLS motivate the Bayesian approach where a simple Markov Chain Monte Carlo (MCMC) sampler is applied. The model is then fit to real data of influent and effluent urea concentrations from experiments with biofilms of Escherichia coli. Results from synthetic data aid in interpreting results from real data, where first-order and Michaelis–Menten kinetic models are compared. The method shows potential for general applications requiring biofilm kinetic information.
尿素水解生物膜在医学、工程和科学领域的应用至关重要。需要关于生物膜中尿素分解速率的定量信息来对这些应用进行建模。我们建立了一个新的生物膜中尿素消耗模型,该模型允许不同的动力学,例如一阶或Michaelis–Menten。该模型适用于合成数据,以验证和比较生物膜从业者常用的贝叶斯和非线性最小二乘(NLS)两种方法。NLS的缺点激发了贝叶斯方法,其中应用了简单的马尔可夫链蒙特卡罗(MCMC)采样器。然后将该模型与大肠杆菌生物膜实验的进水和出水尿素浓度的真实数据进行拟合。合成数据的结果有助于解释真实数据的结果,其中对一阶和Michaelis–Menten动力学模型进行了比较。该方法显示出在需要生物膜动力学信息的一般应用中的潜力。
{"title":"Bayesian estimation and uncertainty quantification in models of urea hydrolysis by E. coli biofilms","authors":"B. Jackson, J. Connolly, R. Gerlach, I. Klapper, A. Parker","doi":"10.1080/17415977.2021.1887172","DOIUrl":"https://doi.org/10.1080/17415977.2021.1887172","url":null,"abstract":"ABSTRACT Urea-hydrolysing biofilms are crucial to applications in medicine, engineering, and science. Quantitative information about ureolysis rates in biofilms is required to model these applications. We formulate a novel model of urea consumption in a biofilm that allows different kinetics, for example either first order or Michaelis–Menten. The model is fit to synthetic data to validate and compare two approaches, Bayesian and nonlinear least squares (NLS), commonly used by biofilm practitioners. The shortcomings of NLS motivate the Bayesian approach where a simple Markov Chain Monte Carlo (MCMC) sampler is applied. The model is then fit to real data of influent and effluent urea concentrations from experiments with biofilms of Escherichia coli. Results from synthetic data aid in interpreting results from real data, where first-order and Michaelis–Menten kinetic models are compared. The method shows potential for general applications requiring biofilm kinetic information.","PeriodicalId":54926,"journal":{"name":"Inverse Problems in Science and Engineering","volume":"29 1","pages":"1629 - 1652"},"PeriodicalIF":1.3,"publicationDate":"2021-02-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1080/17415977.2021.1887172","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41889341","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 3
Approximation error model (AEM) approach with hybrid methods in the forward-inverse analysis of the transesterification reaction in 3D-microreactors 三维微反应器中酯交换反应正逆分析的混合方法逼近误差模型(AEM)方法
IF 1.3 4区 工程技术 Q3 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2021-02-22 DOI: 10.1080/17415977.2020.1870973
P. C. Pontes, J. M. C. Junior, C. Naveira-Cotta, Mrigank Tiwari
This work advances the approximation error model approach for the inverse analysis of the biodiesel synthesis using soybean oil and methanol in 3D-microreactors. Two hybrid numerical-analytical approaches of reduced computational cost are considered to offer an approximate forward problem solution for a three-dimensional nonlinear coupled diffusive-convective-reactive model. First, the Generalized Integral Transform Technique (GITT) is applied using approximate non-converged solutions of the 3D model, by adopting low truncation orders in the eigenfunction expansions. Second, the Coupled Integral Equations Approach (CIEA) provides a reduced mathematical model for the average concentrations, which leads to inherently approximate solutions. The AEM approach through the Bayesian framework is illustrated in the simultaneous estimation of kinetic and diffusion coefficients of the transesterification reaction. For this purpose, the fully converged GITT results with higher truncation orders for the 3D partial differential model are employed as reference results to define the approximations errors. The results highlight that either the non-converged solutions via GITT or the reduced model solution obtained via CIEA, when taking into account the model error, are robust and cost-effective alternatives for the inverse analysis of nonlinear convection–diffusion-reaction problems.
本文提出了一种近似误差模型方法,用于3d微反应器中大豆油和甲醇合成生物柴油的反分析。考虑了两种计算成本较低的混合数值解析方法,为三维非线性耦合扩散-对流-反应模型提供了近似的正解。首先,通过在特征函数展开中采用低截断阶,利用三维模型的近似非收敛解,应用广义积分变换技术(GITT);其次,耦合积分方程方法(CIEA)为平均浓度提供了一个简化的数学模型,从而导致固有的近似解。通过贝叶斯框架说明了AEM方法在同时估计酯交换反应的动力学和扩散系数。为此,采用三维偏微分模型具有较高截断阶的完全收敛的GITT结果作为参考结果来定义近似误差。结果表明,在考虑模型误差的情况下,通过git得到的非收敛解或通过CIEA得到的简化模型解都是非线性对流-扩散-反应问题逆分析的鲁棒且经济的替代方案。
{"title":"Approximation error model (AEM) approach with hybrid methods in the forward-inverse analysis of the transesterification reaction in 3D-microreactors","authors":"P. C. Pontes, J. M. C. Junior, C. Naveira-Cotta, Mrigank Tiwari","doi":"10.1080/17415977.2020.1870973","DOIUrl":"https://doi.org/10.1080/17415977.2020.1870973","url":null,"abstract":"This work advances the approximation error model approach for the inverse analysis of the biodiesel synthesis using soybean oil and methanol in 3D-microreactors. Two hybrid numerical-analytical approaches of reduced computational cost are considered to offer an approximate forward problem solution for a three-dimensional nonlinear coupled diffusive-convective-reactive model. First, the Generalized Integral Transform Technique (GITT) is applied using approximate non-converged solutions of the 3D model, by adopting low truncation orders in the eigenfunction expansions. Second, the Coupled Integral Equations Approach (CIEA) provides a reduced mathematical model for the average concentrations, which leads to inherently approximate solutions. The AEM approach through the Bayesian framework is illustrated in the simultaneous estimation of kinetic and diffusion coefficients of the transesterification reaction. For this purpose, the fully converged GITT results with higher truncation orders for the 3D partial differential model are employed as reference results to define the approximations errors. The results highlight that either the non-converged solutions via GITT or the reduced model solution obtained via CIEA, when taking into account the model error, are robust and cost-effective alternatives for the inverse analysis of nonlinear convection–diffusion-reaction problems.","PeriodicalId":54926,"journal":{"name":"Inverse Problems in Science and Engineering","volume":"29 1","pages":"1586 - 1612"},"PeriodicalIF":1.3,"publicationDate":"2021-02-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1080/17415977.2020.1870973","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44754425","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Experimental static data based damage localization of beam-like structures considering axial load 考虑轴向载荷的基于实验静力数据的类梁结构损伤定位
IF 1.3 4区 工程技术 Q3 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2021-02-17 DOI: 10.1080/17415977.2021.1872565
M. Hashemi, R. Izadifard, O. Yazdanpanah
In this research work, a crack diagnosis method for beam-column structures is proposed considering axial load effects through experimental data. The proposed method is employed for the detection of damage locations including single and multiple damage scenarios considering four cases of simply supported beam-column. The results show that the locations of single and multiple damage scenarios can be well recognized in low axial loads considering random noise effect. However, in the vicinity of the critical load, healthy and damaged static data are uncommon. On the other hand, increasing the axial load, especially when it reaches to the critical load, has a negative effect on the static responses that is precisely considered by the proposed damaged index.
在本研究工作中,通过实验数据,提出了一种考虑轴向载荷影响的梁柱结构裂纹诊断方法。将所提出的方法用于检测损伤位置,包括考虑四种简支梁柱情况的单个和多个损伤场景。结果表明,在考虑随机噪声影响的低轴向载荷下,可以很好地识别单个和多个损伤场景的位置。然而,在临界载荷附近,健康和损坏的静态数据并不常见。另一方面,增加轴向载荷,特别是当轴向载荷达到临界载荷时,会对静态响应产生负面影响,而静态响应正是所提出的损伤指数所考虑的。
{"title":"Experimental static data based damage localization of beam-like structures considering axial load","authors":"M. Hashemi, R. Izadifard, O. Yazdanpanah","doi":"10.1080/17415977.2021.1872565","DOIUrl":"https://doi.org/10.1080/17415977.2021.1872565","url":null,"abstract":"In this research work, a crack diagnosis method for beam-column structures is proposed considering axial load effects through experimental data. The proposed method is employed for the detection of damage locations including single and multiple damage scenarios considering four cases of simply supported beam-column. The results show that the locations of single and multiple damage scenarios can be well recognized in low axial loads considering random noise effect. However, in the vicinity of the critical load, healthy and damaged static data are uncommon. On the other hand, increasing the axial load, especially when it reaches to the critical load, has a negative effect on the static responses that is precisely considered by the proposed damaged index.","PeriodicalId":54926,"journal":{"name":"Inverse Problems in Science and Engineering","volume":"29 1","pages":"1729 - 1745"},"PeriodicalIF":1.3,"publicationDate":"2021-02-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1080/17415977.2021.1872565","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47822838","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 2
A simple method of reconstructing a point-like scatterer according to time-dependent acoustic wave propagation 根据时间相关声波传播重构点状散射体的一种简单方法
IF 1.3 4区 工程技术 Q3 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2021-02-09 DOI: 10.1080/17415977.2021.1886290
Bo Chen, Yao Sun
This paper investigates the approximate solutions to the time-dependent acoustic scattering problem with a point-like scatterer under some basic assumptions and provides a simple method to reconstruct the location of the scatterer. The approximations of the solution to the forward scattering problem are analysed utilizing Green's function and the retarded single-layer potential. Then, based on the approximate solutions, a sampling method is proposed to solve the inverse scattering problem for the location of the scatterer. The proposed method is easy to implement since no equations or matrices have to be computed for the reconstruction. Numerical experiments are provided to show the effectiveness of the method.
本文在一些基本假设下研究了类点散射体的含时声散射问题的近似解,并提供了一种重建散射体位置的简单方法。利用格林函数和延迟单层势分析了前向散射问题的近似解。然后,基于近似解,提出了一种采样方法来解决散射体位置的逆散射问题。所提出的方法易于实现,因为重建不需要计算方程或矩阵。数值实验验证了该方法的有效性。
{"title":"A simple method of reconstructing a point-like scatterer according to time-dependent acoustic wave propagation","authors":"Bo Chen, Yao Sun","doi":"10.1080/17415977.2021.1886290","DOIUrl":"https://doi.org/10.1080/17415977.2021.1886290","url":null,"abstract":"This paper investigates the approximate solutions to the time-dependent acoustic scattering problem with a point-like scatterer under some basic assumptions and provides a simple method to reconstruct the location of the scatterer. The approximations of the solution to the forward scattering problem are analysed utilizing Green's function and the retarded single-layer potential. Then, based on the approximate solutions, a sampling method is proposed to solve the inverse scattering problem for the location of the scatterer. The proposed method is easy to implement since no equations or matrices have to be computed for the reconstruction. Numerical experiments are provided to show the effectiveness of the method.","PeriodicalId":54926,"journal":{"name":"Inverse Problems in Science and Engineering","volume":"29 1","pages":"1895 - 1911"},"PeriodicalIF":1.3,"publicationDate":"2021-02-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1080/17415977.2021.1886290","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47555066","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
Inverse scattering problem for detecting a defect in a magnetoelastic layer 检测磁弹性层缺陷的逆散射问题
IF 1.3 4区 工程技术 Q3 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2021-02-09 DOI: 10.1080/17415977.2021.1884246
Khaled M. Elmorabie, Rania R. Yahya
ABSTRACT This work is devoted to studying a direct and inverse scattering problem for a magnetoelastic layer having a defect, in the frame of the electromagnetic theory. In terms of the displacement field over the defect's contour, a coupled system of boundary integral equations is formulated, for magnetically permeable and impermeable defects. To identify the position and size of the defect, an efficient numerical algorithm is developed by using the quasi-Newton iterative method. In order to check the influence of the magnetic field upon the scattering waves from the layer, a series of numerical examples is presented with different noise levels. The results showed that the magnetic field has a sensitive effect on the identification process when the external magnetic field increases, especially for the materials having a high magnetic permeability factor . Also, a special inverse problem for predicting the external applied magnetic field, upon a copper layer having a defect with various sizes, has been performed.
本文致力于在电磁理论的框架内研究具有缺陷的磁弹性层的正散射和反散射问题。根据缺陷轮廓上的位移场,针对透磁和不透磁缺陷,建立了边界积分方程的耦合系统。为了识别缺陷的位置和大小,利用拟牛顿迭代方法开发了一种有效的数值算法。为了检验磁场对层散射波的影响,给出了一系列具有不同噪声水平的数值例子。结果表明,当外部磁场增加时,磁场对识别过程有敏感影响,尤其是对于具有高磁导率因子的材料。此外,已经执行了用于预测在具有各种尺寸的缺陷的铜层上施加的外部磁场的特殊逆问题。
{"title":"Inverse scattering problem for detecting a defect in a magnetoelastic layer","authors":"Khaled M. Elmorabie, Rania R. Yahya","doi":"10.1080/17415977.2021.1884246","DOIUrl":"https://doi.org/10.1080/17415977.2021.1884246","url":null,"abstract":"ABSTRACT This work is devoted to studying a direct and inverse scattering problem for a magnetoelastic layer having a defect, in the frame of the electromagnetic theory. In terms of the displacement field over the defect's contour, a coupled system of boundary integral equations is formulated, for magnetically permeable and impermeable defects. To identify the position and size of the defect, an efficient numerical algorithm is developed by using the quasi-Newton iterative method. In order to check the influence of the magnetic field upon the scattering waves from the layer, a series of numerical examples is presented with different noise levels. The results showed that the magnetic field has a sensitive effect on the identification process when the external magnetic field increases, especially for the materials having a high magnetic permeability factor . Also, a special inverse problem for predicting the external applied magnetic field, upon a copper layer having a defect with various sizes, has been performed.","PeriodicalId":54926,"journal":{"name":"Inverse Problems in Science and Engineering","volume":"29 1","pages":"1864 - 1894"},"PeriodicalIF":1.3,"publicationDate":"2021-02-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1080/17415977.2021.1884246","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43317548","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 3
Inverse eigenvalue problems for discrete gyroscopic systems 离散陀螺系统的特征值反问题
IF 1.3 4区 工程技术 Q3 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2021-02-08 DOI: 10.1080/17415977.2021.1879804
Hairui Zhang, Yongxin Yuan
A discrete gyroscopic system is characterized by first-order ordinary differential equations defined by one symmetric and one skew-symmetric, which system describes the motion of a spinning body containing elastic parts. In this paper, we consider the inverse problems of such system: Given partial spectral data, find a system such that it is of the desired spectral data. The general solution of the problem is given and the best approximation solution to a pair of matrices is provided by QR-decomposition and matrix derivation. In addition, we also consider a special case in which the system operates below the lowest critical speed. The numerical examples show that the proposed method is effective.
离散陀螺系统的特征是由一个对称和一个斜对称定义的一阶常微分方程,该系统描述了包含弹性部件的旋转体的运动。在本文中,我们考虑了这种系统的逆问题:给定部分光谱数据,找到一个系统,使其具有所需的光谱数据。给出了该问题的一般解,并通过QR分解和矩阵推导给出了一对矩阵的最佳逼近解。此外,我们还考虑了一种特殊情况,即系统在最低临界速度以下运行。算例表明,该方法是有效的。
{"title":"Inverse eigenvalue problems for discrete gyroscopic systems","authors":"Hairui Zhang, Yongxin Yuan","doi":"10.1080/17415977.2021.1879804","DOIUrl":"https://doi.org/10.1080/17415977.2021.1879804","url":null,"abstract":"A discrete gyroscopic system is characterized by first-order ordinary differential equations defined by one symmetric and one skew-symmetric, which system describes the motion of a spinning body containing elastic parts. In this paper, we consider the inverse problems of such system: Given partial spectral data, find a system such that it is of the desired spectral data. The general solution of the problem is given and the best approximation solution to a pair of matrices is provided by QR-decomposition and matrix derivation. In addition, we also consider a special case in which the system operates below the lowest critical speed. The numerical examples show that the proposed method is effective.","PeriodicalId":54926,"journal":{"name":"Inverse Problems in Science and Engineering","volume":"29 1","pages":"1746 - 1763"},"PeriodicalIF":1.3,"publicationDate":"2021-02-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1080/17415977.2021.1879804","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42770920","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
A blocking scheme for dimension-robust Gibbs sampling in large-scale image deblurring 一种用于大规模图像去模糊的维鲁棒Gibbs采样的分块方案
IF 1.3 4区 工程技术 Q3 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2021-02-05 DOI: 10.1080/17415977.2021.1880398
Jesse Adams, M. Morzfeld, K. Joyce, M. Howard, A. Luttman
ABSTRACT Among the most significant challenges with using Markov chain Monte Carlo (MCMC) methods for sampling from the posterior distributions of Bayesian inverse problems is the rate at which the sampling becomes computationally intractable, as a function of the number of estimated parameters. In image deblurring, there are many MCMC algorithms in the literature, but few attempt reconstructions for images larger than pixels ( parameters). In quantitative X-ray radiography, used to diagnose dynamic materials experiments, the images can be much larger, leading to problems with millions of parameters. We address this issue and construct a Gibbs sampler via a blocking scheme that leads to a sparse and highly structured posterior precision matrix. The Gibbs sampler naturally exploits the special matrix structure during sampling, making it ‘dimension-robust’, so that its mixing properties are nearly independent of the image size, and generating one sample is computationally feasible. The dimension-robustness enables the characterization of posteriors for large-scale image deblurring problems on modest computational platforms. We demonstrate applicability of this approach by deblurring radiographs of size pixels ( parameters) taken at the Cygnus Dual Beam X-ray Radiography Facility at the U.S. Department of Energy's Nevada National Security Site.
使用马尔可夫链蒙特卡罗(MCMC)方法从贝叶斯反问题的后验分布中抽样的最重要挑战之一是抽样在计算上变得难以处理的速度,作为估计参数数量的函数。在图像去模糊方面,文献中有许多MCMC算法,但很少尝试对大于像素(参数)的图像进行重建。在用于诊断动态材料实验的定量x射线摄影中,图像可能要大得多,导致数百万个参数出现问题。我们解决了这个问题,并通过阻塞方案构建吉布斯采样器,该方案导致稀疏和高度结构化的后验精度矩阵。吉布斯采样器在采样过程中自然地利用了特殊的矩阵结构,使其具有“维鲁棒性”,因此其混合特性几乎与图像大小无关,并且生成一个样本在计算上是可行的。维鲁棒性可以在适度的计算平台上对大规模图像去模糊问题的后验进行表征。我们通过在美国能源部内华达国家安全基地的天鹅座双束x射线摄影设备上对大小像素(参数)的x射线照片进行去模糊处理,证明了这种方法的适用性。
{"title":"A blocking scheme for dimension-robust Gibbs sampling in large-scale image deblurring","authors":"Jesse Adams, M. Morzfeld, K. Joyce, M. Howard, A. Luttman","doi":"10.1080/17415977.2021.1880398","DOIUrl":"https://doi.org/10.1080/17415977.2021.1880398","url":null,"abstract":"ABSTRACT Among the most significant challenges with using Markov chain Monte Carlo (MCMC) methods for sampling from the posterior distributions of Bayesian inverse problems is the rate at which the sampling becomes computationally intractable, as a function of the number of estimated parameters. In image deblurring, there are many MCMC algorithms in the literature, but few attempt reconstructions for images larger than pixels ( parameters). In quantitative X-ray radiography, used to diagnose dynamic materials experiments, the images can be much larger, leading to problems with millions of parameters. We address this issue and construct a Gibbs sampler via a blocking scheme that leads to a sparse and highly structured posterior precision matrix. The Gibbs sampler naturally exploits the special matrix structure during sampling, making it ‘dimension-robust’, so that its mixing properties are nearly independent of the image size, and generating one sample is computationally feasible. The dimension-robustness enables the characterization of posteriors for large-scale image deblurring problems on modest computational platforms. We demonstrate applicability of this approach by deblurring radiographs of size pixels ( parameters) taken at the Cygnus Dual Beam X-ray Radiography Facility at the U.S. Department of Energy's Nevada National Security Site.","PeriodicalId":54926,"journal":{"name":"Inverse Problems in Science and Engineering","volume":"29 1","pages":"1789 - 1810"},"PeriodicalIF":1.3,"publicationDate":"2021-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1080/17415977.2021.1880398","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47491214","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 3
Solution of the symmetric band partial inverse eigenvalue problem for the damped mass spring system 阻尼质量弹簧系统对称带偏反特征值问题的求解
IF 1.3 4区 工程技术 Q3 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2021-02-03 DOI: 10.1080/17415977.2021.1876688
S. Rakshit, B. Datta
ABSTRACT The structured partial quadratic inverse eigenvalue problem (SPQIEP) is to construct the structured quadratic matrix polynomial using the partial eigendata. The structures arising in physical applications include symmetry, band (tridiagonal, diagonal, pentagonal) etc. The construction of the structured matrix polynomial is the most difficult aspect of this problem and the research on structured inverse eigenvalue problem is rare. In this paper, the symmetric band partial quadratic inverse eigenvalue problem (SBPQIEP) for the damped mass spring system is considered. This problem concerns in finding the symmetric band matrices , and C with bandwidth p from m ( ) prescribed eigenpairs so that the corresponding quadratic matrix polynomial has the given eigenpairs as its eigenvalues and eigenvectors. In general, SBPQIEP is very hard to solve due to the additional band structure constraint. We propose a novel method, based on the matrix-vectorization and Kronecker product of matrices for solving this problem. Furthermore, explicit expressions for general solutions are presented. Numerical experiments on a spring mass problem are presented to illustrate the applicability and the practical usefulness of the proposed method.
结构化部分二次特征值反问题(SPQIEP)是利用部分特征数据构造结构化二次矩阵多项式的问题。在物理应用中出现的结构包括对称、带状(三对角线、对角线、五角线)等。结构化矩阵多项式的构造是这一问题的难点,对结构化特征值反问题的研究较少。研究了阻尼质量弹簧系统的对称带偏二次特征值反问题。这个问题涉及到从m()个规定的特征对中找到对称带矩阵和带宽为p的C,使相应的二次矩阵多项式以给定的特征对作为其特征值和特征向量。一般来说,由于额外的频带结构约束,SBPQIEP很难求解。我们提出了一种基于矩阵矢量化和矩阵的Kronecker积的新方法来解决这一问题。进一步给出了一般解的显式表达式。通过一个弹簧质量问题的数值实验,说明了该方法的适用性和实用性。
{"title":"Solution of the symmetric band partial inverse eigenvalue problem for the damped mass spring system","authors":"S. Rakshit, B. Datta","doi":"10.1080/17415977.2021.1876688","DOIUrl":"https://doi.org/10.1080/17415977.2021.1876688","url":null,"abstract":"ABSTRACT The structured partial quadratic inverse eigenvalue problem (SPQIEP) is to construct the structured quadratic matrix polynomial using the partial eigendata. The structures arising in physical applications include symmetry, band (tridiagonal, diagonal, pentagonal) etc. The construction of the structured matrix polynomial is the most difficult aspect of this problem and the research on structured inverse eigenvalue problem is rare. In this paper, the symmetric band partial quadratic inverse eigenvalue problem (SBPQIEP) for the damped mass spring system is considered. This problem concerns in finding the symmetric band matrices , and C with bandwidth p from m ( ) prescribed eigenpairs so that the corresponding quadratic matrix polynomial has the given eigenpairs as its eigenvalues and eigenvectors. In general, SBPQIEP is very hard to solve due to the additional band structure constraint. We propose a novel method, based on the matrix-vectorization and Kronecker product of matrices for solving this problem. Furthermore, explicit expressions for general solutions are presented. Numerical experiments on a spring mass problem are presented to illustrate the applicability and the practical usefulness of the proposed method.","PeriodicalId":54926,"journal":{"name":"Inverse Problems in Science and Engineering","volume":"11 1","pages":"1497 - 1518"},"PeriodicalIF":1.3,"publicationDate":"2021-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1080/17415977.2021.1876688","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"59998275","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
Comparative studies on the criteria for regularization parameter selection based on moving force identification 基于移动力识别的正则化参数选择准则的比较研究
IF 1.3 4区 工程技术 Q3 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2021-02-01 DOI: 10.1080/17415977.2020.1781848
Zhen Chen, Zhen Wang, Zhihao Wang, T. Chan
The studies on inverse problems exist extensively in aerospace, mechanical, identification, detection, scanning imaging and other fields. Its ill-posed characteristics often lead to large oscillations in the solution of the inverse problem. In this study, the truncated generalized singular value decomposition (TGSVD) method is introduced to identify two kinds of moving forces, single and multi-axial forces. The truncating point is the most influential regularization parameter of TGSVD, which is initially selected by two classic regularization parameter selection criteria, namely, the L-curve criterion and the generalized cross-validation (GCV) criterion. Due to numerical non-uniqueness and noise disturbance in moving force identification (MFI), numerical simulation results show that neither of the two criteria can effectively help select the optimal truncating point of TGSVD. Hence, a relative percentage error (RPE) criterion is proposed for selecting the truncating point of TGSVD. Comparative studies show that the RPE criterion can be used to select the optimal truncating point of TGSVD more accurately against the GCV criterion and L-curve criterion. Moreover, the RPE criterion can be used to reflect the connections between certain properties and the ill-posedness problem existing in MFI, which should be adopted priority for the optimal truncating point selection of TGSVD.
反问题研究广泛存在于航空航天、机械、识别、检测、扫描成像等领域。它的不适定特性经常导致反问题求解中的大振荡。本文引入截断广义奇异值分解(TGSVD)方法来识别单轴力和多轴力这两种运动力。截断点是TGSVD最具影响力的正则化参数,它最初是由两个经典的正则参数选择标准选择的,即L-曲线标准和广义交叉验证(GCV)标准。由于运动力识别中的数值非唯一性和噪声干扰,数值模拟结果表明,这两个准则都不能有效地帮助选择TGSVD的最优截断点。因此,提出了一种选择TGSVD截断点的相对百分比误差(RPE)准则。对比研究表明,与GCV准则和L-曲线准则相比,RPE准则可以更准确地选择TGSVD的最佳截断点。此外,RPE准则可以用来反映某些性质与MFI中存在的不适定性问题之间的联系,这应该优先用于TGSVD的最优截断点选择。
{"title":"Comparative studies on the criteria for regularization parameter selection based on moving force identification","authors":"Zhen Chen, Zhen Wang, Zhihao Wang, T. Chan","doi":"10.1080/17415977.2020.1781848","DOIUrl":"https://doi.org/10.1080/17415977.2020.1781848","url":null,"abstract":"The studies on inverse problems exist extensively in aerospace, mechanical, identification, detection, scanning imaging and other fields. Its ill-posed characteristics often lead to large oscillations in the solution of the inverse problem. In this study, the truncated generalized singular value decomposition (TGSVD) method is introduced to identify two kinds of moving forces, single and multi-axial forces. The truncating point is the most influential regularization parameter of TGSVD, which is initially selected by two classic regularization parameter selection criteria, namely, the L-curve criterion and the generalized cross-validation (GCV) criterion. Due to numerical non-uniqueness and noise disturbance in moving force identification (MFI), numerical simulation results show that neither of the two criteria can effectively help select the optimal truncating point of TGSVD. Hence, a relative percentage error (RPE) criterion is proposed for selecting the truncating point of TGSVD. Comparative studies show that the RPE criterion can be used to select the optimal truncating point of TGSVD more accurately against the GCV criterion and L-curve criterion. Moreover, the RPE criterion can be used to reflect the connections between certain properties and the ill-posedness problem existing in MFI, which should be adopted priority for the optimal truncating point selection of TGSVD.","PeriodicalId":54926,"journal":{"name":"Inverse Problems in Science and Engineering","volume":"29 1","pages":"153 - 173"},"PeriodicalIF":1.3,"publicationDate":"2021-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1080/17415977.2020.1781848","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42279517","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 6
期刊
Inverse Problems in Science and Engineering
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1