首页 > 最新文献

IMA Journal of Applied Mathematics最新文献

英文 中文
Identifying a response parameter in a model of brain tumor evolution under therapy 确定治疗下脑肿瘤进化模型中的反应参数
IF 1.2 4区 数学 Q3 Mathematics Pub Date : 2023-04-04 DOI: 10.1093/imamat/hxad013
A nonlinear conjugate gradient method is derived for the inverse problem of identifying a treatment parameter in a nonlinear model of reaction-diffusion type corresponding to the evolution of brain tumors under therapy. The treatment parameter is reconstructed from additional information about the tumour taken at a fixed instance of time. Well-posedness of the direct problems used in the iterative method is outlined as well as uniqueness of a solution to the inverse problem. Moreover, the parameter identification is recast as the minimization of a Tikhonov type functional and the existence of a minimizer to this functional is shown. Finite difference discretization of the space and time derivatives are employed for the numerical implementation. Numerical simulations on full 3-dimensional brain data is included showing that information about a spacewise dependent treatment parameter can be recovered in a stable way.
针对脑肿瘤治疗过程中反应扩散型非线性模型中治疗参数的反演问题,导出了一种非线性共轭梯度法。治疗参数是根据在固定时间实例中获取的关于肿瘤的附加信息重建的。概述了迭代法中使用的直接问题的适定性以及反问题解的唯一性。此外,参数辨识被重新定义为Tikhonov型泛函的最小化,并证明了该泛函的最小化存在性。采用空间导数和时间导数的有限差分离散化方法进行数值实现。对全三维脑数据进行了数值模拟,结果表明,空间相关治疗参数的信息可以稳定地恢复。
{"title":"Identifying a response parameter in a model of brain tumor evolution under therapy","authors":"","doi":"10.1093/imamat/hxad013","DOIUrl":"https://doi.org/10.1093/imamat/hxad013","url":null,"abstract":"\u0000 A nonlinear conjugate gradient method is derived for the inverse problem of identifying a treatment parameter in a nonlinear model of reaction-diffusion type corresponding to the evolution of brain tumors under therapy. The treatment parameter is reconstructed from additional information about the tumour taken at a fixed instance of time. Well-posedness of the direct problems used in the iterative method is outlined as well as uniqueness of a solution to the inverse problem. Moreover, the parameter identification is recast as the minimization of a Tikhonov type functional and the existence of a minimizer to this functional is shown. Finite difference discretization of the space and time derivatives are employed for the numerical implementation. Numerical simulations on full 3-dimensional brain data is included showing that information about a spacewise dependent treatment parameter can be recovered in a stable way.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.2,"publicationDate":"2023-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"61109929","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
Global threshold dynamics of a spatial chemotactic mosquito-borne disease model 空间趋化性蚊媒疾病模型的全局阈值动力学
IF 1.2 4区 数学 Q3 Mathematics Pub Date : 2023-04-04 DOI: 10.1093/imamat/hxad009
Kai Wang, Hao Wang, Hongyong Zhao
It is natural that mosquitoes move toward high human population density and environmental heterogeneity plays a pivotal role on disease transmission, and thus we formulate and analyze a mosquito-borne disease model with chemotaxis and spatial heterogeneity. The global existence and boundedness of solutions are proven to guarantee the solvability of the model and is challenging due to the model complexity. Under appropriate conditions, we demonstrate the disease-free equilibrium is globally asymptotically stable provided that the basic reproduction number $mathcal {R}_0$ is less than one, and the system is uniformly persistent and admits at least one endemic equilibrium if $mathcal {R}_0$ is greater than one. Furthermore, we numerically explore the impacts of chemotactic effect, spatial heterogeneity and dispersal rates of infected individuals to provide a clear picture on disease severity. In particular, the mosquito chemotaxis causes disease mild in some regions but severe in others, which suggests developing targeted strategies to control mosquitoes in specific locations and achieves a deep understanding on the chemotaxis.
蚊虫向高密度人群迁移是自然规律,环境异质性在疾病传播中起着关键作用,因此我们建立了具有趋化性和空间异质性的蚊媒疾病模型并进行了分析。证明了解的全局存在性和有界性,保证了模型的可解性,但由于模型的复杂性,对模型的可解性提出了挑战。在适当的条件下,证明了当基本繁殖数$mathcal {R}_0$小于1时,无病平衡点是全局渐近稳定的;当基本繁殖数$mathcal {R}_0$大于1时,系统是一致持久的,并且至少存在一个地方性平衡点。此外,我们在数值上探讨了趋化效应、空间异质性和感染个体扩散率的影响,以提供疾病严重程度的清晰图像。特别是,蚊子趋化性在一些地区引起轻度疾病,而在另一些地区引起严重疾病,这建议制定针对性的策略来控制特定地点的蚊子,并对趋化性有深入的了解。
{"title":"Global threshold dynamics of a spatial chemotactic mosquito-borne disease model","authors":"Kai Wang, Hao Wang, Hongyong Zhao","doi":"10.1093/imamat/hxad009","DOIUrl":"https://doi.org/10.1093/imamat/hxad009","url":null,"abstract":"\u0000 It is natural that mosquitoes move toward high human population density and environmental heterogeneity plays a pivotal role on disease transmission, and thus we formulate and analyze a mosquito-borne disease model with chemotaxis and spatial heterogeneity. The global existence and boundedness of solutions are proven to guarantee the solvability of the model and is challenging due to the model complexity. Under appropriate conditions, we demonstrate the disease-free equilibrium is globally asymptotically stable provided that the basic reproduction number $mathcal {R}_0$ is less than one, and the system is uniformly persistent and admits at least one endemic equilibrium if $mathcal {R}_0$ is greater than one. Furthermore, we numerically explore the impacts of chemotactic effect, spatial heterogeneity and dispersal rates of infected individuals to provide a clear picture on disease severity. In particular, the mosquito chemotaxis causes disease mild in some regions but severe in others, which suggests developing targeted strategies to control mosquitoes in specific locations and achieves a deep understanding on the chemotaxis.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.2,"publicationDate":"2023-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46031581","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
Scattering in a partially open waveguide: the forward problem 部分开放波导中的散射:前向问题
IF 1.2 4区 数学 Q3 Mathematics Pub Date : 2023-02-02 DOI: 10.1093/imamat/hxad004
L. Bourgeois, S. Fliss, Jean-François Fritsch, C. Hazard, A. Recoquillay
This paper is dedicated to an acoustic scattering problem in a two-dimensional partially open waveguide, in the sense that the left part of the waveguide is closed, that is with a bounded cross-section, while the right part is bounded in the transverse direction by some Perfectly Matched Layers that mimic the situation of an open waveguide, that is with an unbounded cross-section. We prove well-posedness of such scattering problem in the Fredholm sense (uniqueness implies existence) and exhibit the asymptotic behaviour of the solution in the longitudinal direction with the help of the Kondratiev approach. Having in mind the numerical computation of the solution, we also propose some transparent boundary conditions in such longitudinal direction, based on Dirichlet-to-Neumann operators. After proving that such artificial conditions actually enable us to approximate the exact solution, some numerical experiments illustrate the quality of such approximation.
本文致力于研究二维部分开放波导中的声散射问题,即波导的左侧部分是封闭的,即具有有界横截面,而右侧部分在横向方向上由一些模拟开放波导情况的完全匹配层限定,即具有无界横截面。我们证明了这种散射问题在Fredholm意义上的适定性(唯一性意味着存在),并借助Kondratiev方法展示了解在纵向上的渐近性。考虑到解的数值计算,我们还基于Dirichlet到Neumann算子,提出了在这种纵向方向上的一些透明边界条件。在证明了这种人为条件实际上使我们能够近似精确解之后,一些数值实验说明了这种近似的质量。
{"title":"Scattering in a partially open waveguide: the forward problem","authors":"L. Bourgeois, S. Fliss, Jean-François Fritsch, C. Hazard, A. Recoquillay","doi":"10.1093/imamat/hxad004","DOIUrl":"https://doi.org/10.1093/imamat/hxad004","url":null,"abstract":"\u0000 This paper is dedicated to an acoustic scattering problem in a two-dimensional partially open waveguide, in the sense that the left part of the waveguide is closed, that is with a bounded cross-section, while the right part is bounded in the transverse direction by some Perfectly Matched Layers that mimic the situation of an open waveguide, that is with an unbounded cross-section. We prove well-posedness of such scattering problem in the Fredholm sense (uniqueness implies existence) and exhibit the asymptotic behaviour of the solution in the longitudinal direction with the help of the Kondratiev approach. Having in mind the numerical computation of the solution, we also propose some transparent boundary conditions in such longitudinal direction, based on Dirichlet-to-Neumann operators. After proving that such artificial conditions actually enable us to approximate the exact solution, some numerical experiments illustrate the quality of such approximation.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.2,"publicationDate":"2023-02-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48882379","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
Correction to: A degenerating convection–diffusion system modelling froth flotation with drainage 修正:一个退化的对流-扩散系统,模拟了带排水的泡沫浮选
4区 数学 Q3 Mathematics Pub Date : 2023-02-01 DOI: 10.1093/imamat/hxad001
{"title":"Correction to: A degenerating convection–diffusion system modelling froth flotation with drainage","authors":"","doi":"10.1093/imamat/hxad001","DOIUrl":"https://doi.org/10.1093/imamat/hxad001","url":null,"abstract":"","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136297796","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
On global in time self-similar solutions of Smoluchowski equation with multiplicative kernel 带乘核的Smoluchowski方程的全局时间自相似解
IF 1.2 4区 数学 Q3 Mathematics Pub Date : 2022-12-23 DOI: 10.1093/imamat/hxad012
G. Breschi, M. Fontelos
We study the similarity solutions (SS) of Smoluchowski coagulation equation with multiplicative kernel $K(x,y)=(xy)^{s}$ for $s
我们研究了具有乘法核$K(x,y)=(xy)^{s}$的Smoluchowski凝聚方程的相似解(SS)。当$s<0$时,SS由三个具有不同渐近行为的区域组成。适当的匹配产生了解的全局描述,该全局描述由伽玛分布尾、由对数正态分布描述的中间区域和解在原点附近快速衰减到零的区域组成。当$sinleft(0,frac{1}{2}right)$时,SS在原点是无界的。它还呈现了三个区域:伽马分布尾部、幂类(或帕累托分布)衰减的中间区域和靠近奇点发生的原点的区域。最后,对Smoluchowski方程的全数值模拟验证了我们的理论结果,并表明了自相似域解的收敛性。
{"title":"On global in time self-similar solutions of Smoluchowski equation with multiplicative kernel","authors":"G. Breschi, M. Fontelos","doi":"10.1093/imamat/hxad012","DOIUrl":"https://doi.org/10.1093/imamat/hxad012","url":null,"abstract":"\u0000 We study the similarity solutions (SS) of Smoluchowski coagulation equation with multiplicative kernel $K(x,y)=(xy)^{s}$ for $s<frac {1}{2}$. When $s<0$ , the SS consists of three regions with distinct asymptotic behaviours. The appropriate matching yields a global description of the solution consisting of a Gamma distribution tail, an intermediate region described by a lognormal distribution and a region of very fast decay of the solutions to zero near the origin. When $sin left ( 0,frac {1}{2}right ) $, the SS is unbounded at the origin. It also presents three regions: a Gamma distribution tail, an intermediate region of power-like (or Pareto distribution) decay and the region close to the origin where a singularity occurs. Finally, full numerical simulations of Smoluchowski equation serve to verify our theoretical results and show the convergence of solutions to the selfsimilar regime.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.2,"publicationDate":"2022-12-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44728371","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
Numerical methods and hypoexponential approximations for gamma distributed delay differential equations. 分布延迟微分方程的数值方法和次指数逼近。
IF 1.2 4区 数学 Q3 Mathematics Pub Date : 2022-12-13 eCollection Date: 2022-12-01 DOI: 10.1093/imamat/hxac027
Tyler Cassidy, Peter Gillich, Antony R Humphries, Christiaan H van Dorp

Gamma distributed delay differential equations (DDEs) arise naturally in many modelling applications. However, appropriate numerical methods for generic gamma distributed DDEs have not previously been implemented. Modellers have therefore resorted to approximating the gamma distribution with an Erlang distribution and using the linear chain technique to derive an equivalent system of ordinary differential equations (ODEs). In this work, we address the lack of appropriate numerical tools for gamma distributed DDEs in two ways. First, we develop a functional continuous Runge-Kutta (FCRK) method to numerically integrate the gamma distributed DDE without resorting to Erlang approximation. We prove the fourth-order convergence of the FCRK method and perform numerical tests to demonstrate the accuracy of the new numerical method. Nevertheless, FCRK methods for infinite delay DDEs are not widely available in existing scientific software packages. As an alternative approach to solving gamma distributed DDEs, we also derive a hypoexponential approximation of the gamma distributed DDE. This hypoexponential approach is a more accurate approximation of the true gamma distributed DDE than the common Erlang approximation but, like the Erlang approximation, can be formulated as a system of ODEs and solved numerically using standard ODE software. Using our FCRK method to provide reference solutions, we show that the common Erlang approximation may produce solutions that are qualitatively different from the underlying gamma distributed DDE. However, the proposed hypoexponential approximations do not have this limitation. Finally, we apply our hypoexponential approximations to perform statistical inference on synthetic epidemiological data to illustrate the utility of the hypoexponential approximation.

伽马分布延迟微分方程(DDEs)在许多建模应用中自然出现。然而,对于一般的伽马分布DDEs,以前没有实现适当的数值方法。因此,建模者采用厄朗分布近似伽马分布,并使用线性链技术推导出一个等效的常微分方程(ode)系统。在这项工作中,我们通过两种方式解决了伽马分布DDEs缺乏适当的数值工具的问题。首先,我们开发了一种泛函连续龙格-库塔(FCRK)方法来对伽马分布DDE进行数值积分,而不需要借助于Erlang近似。我们证明了FCRK方法的四阶收敛性,并通过数值试验验证了新数值方法的准确性。然而,对于无限延迟DDEs的FCRK方法在现有的科学软件包中并没有广泛使用。作为求解伽马分布DDE的另一种方法,我们还推导了伽马分布DDE的次指数近似。与常见的Erlang近似相比,这种次指数方法是对真正的伽玛分布DDE的更精确的近似,但是,与Erlang近似一样,可以将其表述为一个ODE系统,并使用标准ODE软件进行数值求解。使用我们的FCRK方法提供参考解决方案,我们表明,常见的Erlang近似可能产生与底层伽玛分布DDE在质量上不同的解决方案。然而,提出的次指数近似没有这种限制。最后,我们应用我们的次指数近似对合成流行病学数据进行统计推断,以说明次指数近似的效用。
{"title":"Numerical methods and hypoexponential approximations for gamma distributed delay differential equations.","authors":"Tyler Cassidy, Peter Gillich, Antony R Humphries, Christiaan H van Dorp","doi":"10.1093/imamat/hxac027","DOIUrl":"10.1093/imamat/hxac027","url":null,"abstract":"<p><p>Gamma distributed delay differential equations (DDEs) arise naturally in many modelling applications. However, appropriate numerical methods for generic gamma distributed DDEs have not previously been implemented. Modellers have therefore resorted to approximating the gamma distribution with an Erlang distribution and using the linear chain technique to derive an equivalent system of ordinary differential equations (ODEs). In this work, we address the lack of appropriate numerical tools for gamma distributed DDEs in two ways. First, we develop a functional continuous Runge-Kutta (FCRK) method to numerically integrate the gamma distributed DDE without resorting to Erlang approximation. We prove the fourth-order convergence of the FCRK method and perform numerical tests to demonstrate the accuracy of the new numerical method. Nevertheless, FCRK methods for infinite delay DDEs are not widely available in existing scientific software packages. As an alternative approach to solving gamma distributed DDEs, we also derive a hypoexponential approximation of the gamma distributed DDE. This hypoexponential approach is a more accurate approximation of the true gamma distributed DDE than the common Erlang approximation but, like the Erlang approximation, can be formulated as a system of ODEs and solved numerically using standard ODE software. Using our FCRK method to provide reference solutions, we show that the common Erlang approximation may produce solutions that are qualitatively different from the underlying gamma distributed DDE. However, the proposed hypoexponential approximations do not have this limitation. Finally, we apply our hypoexponential approximations to perform statistical inference on synthetic epidemiological data to illustrate the utility of the hypoexponential approximation.</p>","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.2,"publicationDate":"2022-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC9850366/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"10667425","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Extensions of the d’Alembert formulae to the half line and the finite interval obtained via the unified transform d’Alembert公式对由统一变换得到的半直线和有限区间的推广
IF 1.2 4区 数学 Q3 Mathematics Pub Date : 2022-10-29 DOI: 10.1093/imamat/hxac030
A. S. Fokas, K. Kalimeris
We derive the solution of the one dimensional wave equation for the Dirichlet and Robin initial-boundary value problems (IBVPs) formulated on the half line and the finite interval, with nonhomogeneous boundary conditions. Although explicit formulas already exist for these problems, the unified transform method provides a convenient framework for deriving different representations of the solutions for these and other types of IBVPs. Specifically, it provides solution formulas in the Fourier space or solutions which constitute the extension of the classical formula of d’Alembert of the initial value problem on the full line. We also derive the solution of the forced wave equation on the half line.
导出了在半直线和有限区间上具有非齐次边界条件的Dirichlet和Robin初边值问题(IBVPs)的一维波动方程的解。虽然这些问题的显式公式已经存在,但统一变换方法为导出这些和其他类型的ibvp解的不同表示提供了一个方便的框架。具体地说,它提供了傅里叶空间中的解公式或构成初值问题的经典达朗贝尔公式在全线上的推广的解。我们还推导了半线上强迫波动方程的解。
{"title":"Extensions of the d’Alembert formulae to the half line and the finite interval obtained via the unified transform","authors":"A. S. Fokas, K. Kalimeris","doi":"10.1093/imamat/hxac030","DOIUrl":"https://doi.org/10.1093/imamat/hxac030","url":null,"abstract":"\u0000 We derive the solution of the one dimensional wave equation for the Dirichlet and Robin initial-boundary value problems (IBVPs) formulated on the half line and the finite interval, with nonhomogeneous boundary conditions. Although explicit formulas already exist for these problems, the unified transform method provides a convenient framework for deriving different representations of the solutions for these and other types of IBVPs. Specifically, it provides solution formulas in the Fourier space or solutions which constitute the extension of the classical formula of d’Alembert of the initial value problem on the full line. We also derive the solution of the forced wave equation on the half line.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.2,"publicationDate":"2022-10-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46109417","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
Acceleration of Gossip Algorithms through the Euler–Poisson–Darboux Equation 利用Euler-Poisson-Darboux方程加速八卦算法
IF 1.2 4区 数学 Q3 Mathematics Pub Date : 2022-10-27 DOI: 10.1093/imamat/hxac029
Raphaël Berthier, Mufan (Bill) Li
Gossip algorithms and their accelerated versions have been studied exclusively in discrete time on graphs. In this work, we take a different approach, and consider the scaling limit of gossip algorithms in both large graphs and large number of iterations. These limits lead to well-known partial differential equations (PDEs) with insightful properties. On lattices, we prove that the non-accelerated gossip algorithm of (??) converges to the heat equation, and the accelerated Jacobi polynomial iteration of (??) converges to the Euler–Poisson–Darboux (EPD) equation — a damped wave equation. Remarkably, with appropriate parameters, the fundamental solution of the EPD equation has the ideal gossip behaviour: a uniform density over an ellipsoid, whose radius increases at a rate proportional to $t$ — the fastest possible rate for locally communicating gossip algorithms. This is in contrast with the heat equation where the density spreads on a typical scale of $sqrt {t}$. Additionally, we provide simulations demonstrating that the gossip algorithms are accurately approximated by their limiting PDEs.
本文研究了离散时间图上的八卦算法及其加速算法。在这项工作中,我们采取了不同的方法,并考虑了八卦算法在大图和大量迭代中的缩放极限。这些极限导致了众所周知的具有深刻性质的偏微分方程(PDEs)。在格上,我们证明了(??)的非加速八卦算法收敛于热方程,(??)的加速雅可比多项式迭代收敛于欧拉-泊松-达布(EPD)方程-一种阻尼波动方程。值得注意的是,在适当的参数下,EPD方程的基本解具有理想的八卦行为:椭球体上的均匀密度,其半径以与$t$成比例的速率增加-这是局部通信八卦算法的最快速率。这与热方程相反,其中密度以典型的尺度$sqrt {t}$扩散。此外,我们还提供了仿真,证明了八卦算法可以通过其极限偏微分方程精确地近似。
{"title":"Acceleration of Gossip Algorithms through the Euler–Poisson–Darboux Equation","authors":"Raphaël Berthier, Mufan (Bill) Li","doi":"10.1093/imamat/hxac029","DOIUrl":"https://doi.org/10.1093/imamat/hxac029","url":null,"abstract":"Gossip algorithms and their accelerated versions have been studied exclusively in discrete time on graphs. In this work, we take a different approach, and consider the scaling limit of gossip algorithms in both large graphs and large number of iterations. These limits lead to well-known partial differential equations (PDEs) with insightful properties. On lattices, we prove that the non-accelerated gossip algorithm of (??) converges to the heat equation, and the accelerated Jacobi polynomial iteration of (??) converges to the Euler–Poisson–Darboux (EPD) equation — a damped wave equation. Remarkably, with appropriate parameters, the fundamental solution of the EPD equation has the ideal gossip behaviour: a uniform density over an ellipsoid, whose radius increases at a rate proportional to $t$ — the fastest possible rate for locally communicating gossip algorithms. This is in contrast with the heat equation where the density spreads on a typical scale of $sqrt {t}$. Additionally, we provide simulations demonstrating that the gossip algorithms are accurately approximated by their limiting PDEs.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.2,"publicationDate":"2022-10-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138538076","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
Optimal analyticity estimates for non-linear active-dissipative evolution equations 非线性有源耗散演化方程的最优分析性估计
IF 1.2 4区 数学 Q3 Mathematics Pub Date : 2022-10-14 DOI: 10.1093/imamat/hxac028
D. Papageorgiou, Y. Smyrlis, R. Tomlin
Active-dissipative evolution equations emerge in a variety of physical and technological applications including liquid film flows, flame propagation, epitaxial film growth in materials manufacturing, to mention a few. They are characterised by three main ingredients: a term producing growth (active), a term providing damping at short length scales (dissipative), and a nonlinear term that transfers energy between modes and crucially produces a nonlinear saturation. The manifestation of these three mechanisms can produce large-time spatiotemporal chaos as evidenced by the Kuramoto-Sivashinsky equation (negative diffusion, fourth order dissipation, and a Burgers nonlinearity), which is arguably the simplest partial differential equation to produce chaos. The exact form of the terms (and in particular their Fourier symbol) determines the type of attractors that the equations possess. The present study considers the spatial analyticity of solutions under the assumption that the equations possess a global attractor. In particular we investigate the spatial analyticity of solutions of a class of one-dimensional evolutionary pseudo-differential equations with Burgers nonlinearity, which are periodic in space, thus generalising the Kuramoto-Sivashinsky equation motivated by both applications and their fundamental mathematical properties. Analyticity is examined by utilising a criterion involving the rate of growth of suitable norms of the $n$-th spatial derivative of the solution, with respect to the spatial variable, as $n$ tends to infinity. An estimate of the rate of growth of the $n$-th spatial derivative is obtained by fine-tuning the spectral method, developed elsewhere. We prove that the solutions are analytic if $gamma $, the order of dissipation of the pseudo-differential operator, is higher than one. We also present numerical evidence suggesting that this is optimal, i.e., if $gamma $ is not larger that one, then the solution is not in general analytic. Extensive numerical experiments are undertaken to confirm the analysis and also to compute the band of analyticity of solutions for a wide range of active-dissipative terms and large spatial periods that support chaotic solutions. These ideas can be applied to a wide class of active-dissipative-dispersive pseudo-differential equations.
主动耗散演化方程出现在各种物理和技术应用中,包括液膜流动、火焰传播、材料制造中的外延膜生长等。它们由三个主要成分组成:一个产生增长的项(主动),一个在短长度尺度上提供阻尼的项(耗散),以及一个在模式之间传递能量并关键地产生非线性饱和的非线性项。这三种机制的表现可以产生大的时空混沌,Kuramoto Sivashinsky方程(负扩散、四阶耗散和Burgers非线性)证明了这一点,它可以说是产生混沌的最简单的偏微分方程。这些项的确切形式(尤其是它们的傅立叶符号)决定了方程所具有的吸引子的类型。本研究考虑了在方程具有全局吸引子的假设下解的空间分析性。特别地,我们研究了一类具有Burgers非线性的一维演化伪微分方程解的空间分析性,该方程在空间中是周期性的,从而推广了Kuramoto Sivashinsky方程,其动机是应用及其基本数学性质。分析性是通过使用一个标准来检查的,该标准涉及解决方案的第$n$个空间导数的适当范数相对于空间变量的增长率,因为$n$趋于无穷大。通过微调其他地方开发的光谱方法,可以获得第$n$-个空间导数增长率的估计值。我们证明了如果伪微分算子的耗散阶$gamma$高于1,则解是解析的。我们还提供了数字证据,表明这是最优的,即,如果$gamma$不大于1,则该解不是一般的解析解。进行了大量的数值实验来证实分析,并计算支持混沌解的大范围主动耗散项和大空间周期的解的分析性带。这些思想可以应用于一类广泛的有源耗散色散伪微分方程。
{"title":"Optimal analyticity estimates for non-linear active-dissipative evolution equations","authors":"D. Papageorgiou, Y. Smyrlis, R. Tomlin","doi":"10.1093/imamat/hxac028","DOIUrl":"https://doi.org/10.1093/imamat/hxac028","url":null,"abstract":"\u0000 Active-dissipative evolution equations emerge in a variety of physical and technological applications including liquid film flows, flame propagation, epitaxial film growth in materials manufacturing, to mention a few. They are characterised by three main ingredients: a term producing growth (active), a term providing damping at short length scales (dissipative), and a nonlinear term that transfers energy between modes and crucially produces a nonlinear saturation. The manifestation of these three mechanisms can produce large-time spatiotemporal chaos as evidenced by the Kuramoto-Sivashinsky equation (negative diffusion, fourth order dissipation, and a Burgers nonlinearity), which is arguably the simplest partial differential equation to produce chaos. The exact form of the terms (and in particular their Fourier symbol) determines the type of attractors that the equations possess. The present study considers the spatial analyticity of solutions under the assumption that the equations possess a global attractor. In particular we investigate the spatial analyticity of solutions of a class of one-dimensional evolutionary pseudo-differential equations with Burgers nonlinearity, which are periodic in space, thus generalising the Kuramoto-Sivashinsky equation motivated by both applications and their fundamental mathematical properties. Analyticity is examined by utilising a criterion involving the rate of growth of suitable norms of the $n$-th spatial derivative of the solution, with respect to the spatial variable, as $n$ tends to infinity. An estimate of the rate of growth of the $n$-th spatial derivative is obtained by fine-tuning the spectral method, developed elsewhere. We prove that the solutions are analytic if $gamma $, the order of dissipation of the pseudo-differential operator, is higher than one. We also present numerical evidence suggesting that this is optimal, i.e., if $gamma $ is not larger that one, then the solution is not in general analytic. Extensive numerical experiments are undertaken to confirm the analysis and also to compute the band of analyticity of solutions for a wide range of active-dissipative terms and large spatial periods that support chaotic solutions. These ideas can be applied to a wide class of active-dissipative-dispersive pseudo-differential equations.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.2,"publicationDate":"2022-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45254697","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
Singular fourth-order Sturm–Liouville operators and acoustic black holes 奇异四阶Sturm–Liouville算子与声学黑洞
IF 1.2 4区 数学 Q3 Mathematics Pub Date : 2022-10-08 DOI: 10.1093/imamat/hxac021
B. Belinskiy, D. Hinton, R. Nichols
We derive conditions for a one-term fourth-order Sturm–Liouville operator on a finite interval with one singular endpoint to have essential spectrum equal to $[0,infty )$ or $varnothing $. Of particular usefulness are Kummer–Liouville transformations which have been a valuable tool in the study of second-order equations. Applications to a mechanical beam with a thickness tapering to zero at one of the endpoints are considered. When the thickness $2h$ satisfies $c_1x^{nu }leq h(x)leq c_2x^{nu }$, we show that the essential spectrum is empty if and only if $nu < 2$. As a final application, we consider a tapered beam on a Winkler foundation and derive sufficient conditions on the beam thickness and the foundational rigidity to guarantee the essential spectrum is equal to $[0,infty )$.
我们导出了具有一个奇异端点的有限区间上一项四阶Sturm–Liouville算子的本质谱等于$[0,infty)$或$varnote$。Kummer–Liouville变换特别有用,它是研究二阶方程的一个有价值的工具。考虑了在其中一个端点处厚度逐渐变为零的机械梁的应用。当厚度$2h$满足$c_1x^{nu}leq h(x)leq c_2x^{ nu}$,我们证明了当且仅当$nu<2$时,本质谱是空的。作为最后的应用,我们考虑了Winkler基础上的锥形梁,并推导了梁厚度和基础刚度的充分条件,以确保基本谱等于$[0,infty)$。
{"title":"Singular fourth-order Sturm–Liouville operators and acoustic black holes","authors":"B. Belinskiy, D. Hinton, R. Nichols","doi":"10.1093/imamat/hxac021","DOIUrl":"https://doi.org/10.1093/imamat/hxac021","url":null,"abstract":"\u0000 We derive conditions for a one-term fourth-order Sturm–Liouville operator on a finite interval with one singular endpoint to have essential spectrum equal to $[0,infty )$ or $varnothing $. Of particular usefulness are Kummer–Liouville transformations which have been a valuable tool in the study of second-order equations. Applications to a mechanical beam with a thickness tapering to zero at one of the endpoints are considered. When the thickness $2h$ satisfies $c_1x^{nu }leq h(x)leq c_2x^{nu }$, we show that the essential spectrum is empty if and only if $nu < 2$. As a final application, we consider a tapered beam on a Winkler foundation and derive sufficient conditions on the beam thickness and the foundational rigidity to guarantee the essential spectrum is equal to $[0,infty )$.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.2,"publicationDate":"2022-10-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49351707","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
期刊
IMA Journal of Applied Mathematics
全部 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