首页 > 最新文献

Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique最新文献

英文 中文
The Morozov's principle applied to data assimilation problems 莫罗佐夫原理应用于数据同化问题
IF 1.9 3区 数学 Q2 Mathematics Pub Date : 2022-07-08 DOI: 10.1051/m2an/2022061
L. Bourgeois, J. Dardé
This paper is focused on the Morozov’s principle applied to an abstract data assimilation framework, with particular attention to three simple examples: the data assimilation problem for the Laplace equation, the Cauchy problem for the Laplace equation and the data assimilation problem for the heat equation. Those ill-posed problems are regularized with the help of a mixed type formulation which is proved to be equivalent to a Tikhonov regularization applied to a well-chosen operator. The main issue is that such operator may not have a dense range, which makes it necessary to extend well-known results related to the Morozov’s choice of the regularization parameter to that unusual situation. The solution which satisfies the Morozov’s principle is computed with the help of the duality in optimization, possibly by forcing the solution to satisfy given a priori constraints. Some numerical results in two dimensions are proposed in the case of the data assimilation problem for the Laplace equation.
本文重点讨论了Morozov原理在抽象数据同化框架中的应用,并特别注意了三个简单的例子:拉普拉斯方程的数据同化问题、拉普拉斯方程的Cauchy问题和热方程的数据同化问题。这些不适定问题在混合型公式的帮助下被正则化,该公式被证明等同于应用于选定算子的吉洪诺夫正则化。主要的问题是,这样的算子可能没有一个密集的范围,这使得有必要将众所周知的与Morozov选择正则化参数相关的结果扩展到这种不寻常的情况。在优化的对偶性的帮助下,可能通过强迫解满足给定的先验约束来计算满足Morozov原理的解。对于拉普拉斯方程的数据同化问题,给出了二维上的一些数值结果。
{"title":"The Morozov's principle applied to data assimilation problems","authors":"L. Bourgeois, J. Dardé","doi":"10.1051/m2an/2022061","DOIUrl":"https://doi.org/10.1051/m2an/2022061","url":null,"abstract":"This paper is focused on the Morozov’s principle applied to an abstract data assimilation framework, with particular attention to three simple examples: the data assimilation problem for the Laplace equation, the Cauchy problem for the Laplace equation and the data assimilation problem for the heat equation. Those ill-posed problems are regularized with the help of a mixed type formulation which is proved to be equivalent to a Tikhonov regularization applied to a well-chosen operator. The main issue is that such operator may not have a dense range, which makes it necessary to extend well-known results related to the Morozov’s choice of the regularization parameter to that unusual situation. The solution which satisfies the Morozov’s principle is computed with the help of the duality in optimization, possibly by forcing the solution to satisfy given a priori constraints. Some numerical results in two dimensions are proposed in the case of the data assimilation problem for the Laplace equation.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-07-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"74187566","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Acoustic waveguide with a dissipative inclusion 具有耗散内含物的声波导
IF 1.9 3区 数学 Q2 Mathematics Pub Date : 2022-07-01 DOI: 10.1051/m2an/2023070
L. Chesnel, J. Heleine, S. Nazarov, J. Taskinen
We consider the propagation of acoustic waves in a waveguide containing a penetrable dissipative inclusion. We prove that as soon as the dissipation, characterized by some coefficient η , is non zero, the scattering solutions are uniquely defined. Additionally, we give an asymptotic expansion of the corresponding scattering matrix when η → 0 + (small dissipation) and when η → +∞ (large dissipation). Surprisingly, at the limit η → +∞, we show that no energy is absorbed by the inclusion. This is due to the so-called skin-effect phenomenon and can be explained by the fact that the field no longer penetrates into the highly dissipative inclusion. These results guarantee that in monomode regime, the amplitude of the reflection coefficient has a global minimum with respect to η . The situation where this minimum is zero, that is when the device acts as a perfect absorber, is particularly interesting for certain applications. However it does not happen in general. In this work, we show how to perturb the geometry of the waveguide to create 2D perfect absorbers in monomode regime. Asymptotic expansions are justified by error estimates and theoretical results are supported by numerical illustrations.
我们考虑声波在含有可穿透的耗散包体的波导中的传播。我们证明了只要耗散不为零,散射解是唯一定义的。此外,我们给出了η→0 +(小耗散)和η→+∞(大耗散)时对应的散射矩阵的渐近展开式。令人惊讶的是,在η→+∞极限处,我们发现包体没有能量被吸收。这是由于所谓的皮肤效应现象,并可以解释的事实,即场不再穿透到高耗散包涵体。这些结果保证了在单模状态下,反射系数的振幅相对于η具有全局最小值。这个最小值为零的情况,即当器件作为一个完美的吸收器时,在某些应用中特别有趣。然而,这种情况通常不会发生。在这项工作中,我们展示了如何扰动波导的几何形状以在单模状态下创建二维完美吸收器。用误差估计证明了渐近展开式的合理性,并用数值实例支持了理论结果。
{"title":"Acoustic waveguide with a dissipative inclusion","authors":"L. Chesnel, J. Heleine, S. Nazarov, J. Taskinen","doi":"10.1051/m2an/2023070","DOIUrl":"https://doi.org/10.1051/m2an/2023070","url":null,"abstract":"We consider the propagation of acoustic waves in a waveguide containing a penetrable dissipative inclusion. We prove that as soon as the dissipation, characterized by some coefficient η , is non zero, the scattering solutions are uniquely defined. Additionally, we give an asymptotic expansion of the corresponding scattering matrix when η → 0 + (small dissipation) and when η → +∞ (large dissipation). Surprisingly, at the limit η → +∞, we show that no energy is absorbed by the inclusion. This is due to the so-called skin-effect phenomenon and can be explained by the fact that the field no longer penetrates into the highly dissipative inclusion. These results guarantee that in monomode regime, the amplitude of the reflection coefficient has a global minimum with respect to η . The situation where this minimum is zero, that is when the device acts as a perfect absorber, is particularly interesting for certain applications. However it does not happen in general. In this work, we show how to perturb the geometry of the waveguide to create 2D perfect absorbers in monomode regime. Asymptotic expansions are justified by error estimates and theoretical results are supported by numerical illustrations.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"85649209","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Sparse Grid reconstructions for Particle-In-Cell methods 基于粒子单元法的稀疏网格重建
IF 1.9 3区 数学 Q2 Mathematics Pub Date : 2022-06-26 DOI: 10.1051/m2an/2022055
C. Guillet, F. Deluzet, G. Fubiani, L. Garrigues, J. Narski
In this article, we propose and analyse Particle-In-Cell (PIC)methods embedding sparse grid reconstruction as those introduced in [1, 2].The sparse grid reconstructions offer a significant improvement on the sta-tistical error of PIC schemes as well as a reduction in the complexity of theproblem providing the electric field. Main results on the convergence of theelectric field interpolant and conservation properties are provided in this pa-per. Besides, tailored sparse grid reconstructions, in the frame of the offsetcombination technique, are proposed to introduce PIC methods with improvedefficiency. The methods are assessed numerically and compared to existing PICschemes thanks to classical benchmarks with remarkable prospects for threedimensional computations.
在本文中,我们提出并分析了嵌入稀疏网格重建的PIC方法,如[1,2]中所介绍的方法。稀疏网格重建不仅显著改善了PIC方案的统计误差,而且降低了提供电场问题的复杂性。本文给出了电场插值的收敛性和守恒性的主要结果。此外,在偏移组合技术的框架下,提出了定制的稀疏网格重建,以提高PIC方法的效率。这些方法进行了数值评估,并与现有的pics方案进行了比较,这要归功于具有显着前景的三维计算的经典基准。
{"title":"Sparse Grid reconstructions for Particle-In-Cell methods","authors":"C. Guillet, F. Deluzet, G. Fubiani, L. Garrigues, J. Narski","doi":"10.1051/m2an/2022055","DOIUrl":"https://doi.org/10.1051/m2an/2022055","url":null,"abstract":"In this article, we propose and analyse Particle-In-Cell (PIC)\u0000\u0000methods embedding sparse grid reconstruction as those introduced in [1, 2].\u0000\u0000The sparse grid reconstructions offer a significant improvement on the sta-\u0000\u0000tistical error of PIC schemes as well as a reduction in the complexity of the\u0000\u0000problem providing the electric field. Main results on the convergence of the\u0000\u0000electric field interpolant and conservation properties are provided in this pa-\u0000\u0000per. Besides, tailored sparse grid reconstructions, in the frame of the offset\u0000\u0000combination technique, are proposed to introduce PIC methods with improved\u0000\u0000efficiency. The methods are assessed numerically and compared to existing PIC\u0000\u0000schemes thanks to classical benchmarks with remarkable prospects for three\u0000\u0000dimensional computations.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-06-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"79353074","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 5
A study of the local dynamics of modified Patankar Dec and higher order modified Patankar-RK methods 修正Patankar- rk方法和高阶修正Patankar- rk方法的局部动力学研究
IF 1.9 3区 数学 Q2 Mathematics Pub Date : 2022-06-15 DOI: 10.1051/m2an/2023053
Thomas Izgin, Philipp Offner
Patankar schemes have attracted increasing interest in recent years because they preserve the positivity of the analytical solution of a production-destruction system (PDS) irrespective of the chosen time step size. Although they are now of great interest, for a long time it was not clear what stability properties such schemes have. Recently a new stability approach based on Lyapunov stability with an extension of the center manifold theorem has been proposed to study the stability properties of positivity preserving time integrators. In this work, we study the stability properties of the classical modified Patankar–Runge–Kutta schemes (MPRK) and the modified Patankar Deferred Correction (MPDeC) approaches. We prove that most of the considered MPRK schemes are stable for any time step size and compute the stability function of MPDeC. We investigate its properties numerically revealing that also most MPDeC are stable irrespective of the chosen time step size. Finally, we verify our theoretical results with numerical simulations.
Patankar方案近年来引起了人们越来越多的兴趣,因为它保持了生产-破坏系统(PDS)的解析解的正性,而与所选择的时间步长无关。尽管它们现在引起了极大的兴趣,但长期以来,人们并不清楚这种方案具有什么样的稳定性。为了研究保正时间积分器的稳定性,本文提出了一种基于李雅普诺夫稳定性和中心流形定理推广的稳定性方法。本文研究了经典修正Patankar - runge - kutta格式(MPRK)和修正Patankar延迟修正(MPDeC)方法的稳定性。我们证明了大多数被考虑的MPRK方案在任何时间步长下都是稳定的,并计算了MPDeC的稳定性函数。我们对其性质进行了数值研究,表明大多数MPDeC与所选择的时间步长无关,都是稳定的。最后,用数值模拟验证了理论结果。
{"title":"A study of the local dynamics of modified Patankar Dec and higher order modified Patankar-RK methods","authors":"Thomas Izgin, Philipp Offner","doi":"10.1051/m2an/2023053","DOIUrl":"https://doi.org/10.1051/m2an/2023053","url":null,"abstract":"Patankar schemes have attracted increasing interest in recent years because they preserve the positivity of the analytical solution of a production-destruction system (PDS) irrespective of the chosen time step size. Although they are now of great interest, for a long time it was not clear what stability properties such schemes have. Recently a new stability approach based on Lyapunov stability with an extension of the center manifold theorem has been proposed to study the stability properties of positivity preserving time integrators. In this work, we study the stability properties of the classical modified Patankar–Runge–Kutta schemes (MPRK) and the modified Patankar Deferred Correction (MPDeC) approaches. We prove that most of the considered MPRK schemes are stable for any time step size and compute the stability function of MPDeC. We investigate its properties numerically revealing that also most MPDeC are stable irrespective of the chosen time step size. Finally, we verify our theoretical results with numerical simulations.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-06-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81105229","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 4
A fast second-order discretization scheme for the linearized Green-Naghdi system with absorbing boundary conditions 具有吸收边界条件的线性化Green-Naghdi系统的快速二阶离散化格式
IF 1.9 3区 数学 Q2 Mathematics Pub Date : 2022-05-27 DOI: 10.1051/m2an/2022051
Gang Pang, Songsong Ji, X. Antoine
In this paper, we present a fully discrete second-order finite-difference scheme with fast evaluation of the convolution involved in the absorbing boundary conditions to solve the one-dimensional linearized Green-Naghdi system. The Pad´e expansion of the square-root function in the complex plane is used to implement the fast convolution. By introducing a constant damping parameter into the governing equations, the convergence analysis is developed when the damping term fulfills some conditions. In addition, the scheme is stable and leads to a highly reduced computational cost and low memory storage. A numerical example is provided to support the theoretical analysis and to illustrate the performance of the fast numerical scheme.
本文给出了一种完全离散二阶有限差分格式,该格式能快速计算吸收边界条件中涉及的卷积,从而求解一维线性化的Green-Naghdi系统。采用复平面平方根函数的Pad´e展开式实现快速卷积。通过在控制方程中引入一个恒定的阻尼参数,给出了阻尼项满足一定条件时的收敛性分析。此外,该方案稳定,大大降低了计算成本和低内存存储。给出了一个数值算例来支持理论分析和说明快速数值格式的性能。
{"title":"A fast second-order discretization scheme for the linearized Green-Naghdi system with absorbing boundary conditions","authors":"Gang Pang, Songsong Ji, X. Antoine","doi":"10.1051/m2an/2022051","DOIUrl":"https://doi.org/10.1051/m2an/2022051","url":null,"abstract":"In this paper, we present a fully discrete second-order finite-difference scheme with fast evaluation of the convolution involved in the absorbing boundary conditions to solve the one-dimensional linearized Green-Naghdi system. The Pad´e expansion of the square-root function in the complex plane is used to implement the fast convolution. By introducing a constant damping parameter into the governing equations, the convergence analysis is developed when the damping term fulfills some conditions. In addition, the scheme is stable and leads to a highly reduced computational cost and low memory storage. A numerical example is provided to support the theoretical analysis and to illustrate the performance of the fast numerical scheme.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-05-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"74500495","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
Computing effective diffusivities in 3D time-dependent chaotic flows with a convergent Lagrangian numerical method 用收敛拉格朗日数值方法计算三维时变混沌流的有效扩散系数
IF 1.9 3区 数学 Q2 Mathematics Pub Date : 2022-05-17 DOI: 10.1051/m2an/2022049
Zhongjian Wang, J. Xin, Zhiwen Zhang
In this paper, we study the convergence analysis for a robust stochastic structure-preserving Lagrangian numerical scheme in computing effective diffusivity of time-dependent chaotic flows, which are modeled by stochastic differential equations (SDEs). Our numerical scheme is based on a splitting method to solve the corresponding SDEs in which the deterministic subproblem is discretized using structure-preserving schemes while the random subproblem is discretized using the Euler-Maruyama scheme. We obtain a sharp and uniform-in-time convergence analysis for the proposed numerical scheme that allows us to accurately compute long-time solutions of the SDEs. As such, we can compute the effective diffusivity for time-dependent chaotic flows. Finally, we present numerical results to demonstrate the accuracy and efficiency of the proposed method in computing effective diffusivity for the time-dependent Arnold-Beltrami-Childress (ABC) flow and Kolmogorov flow in three-dimensional space.
本文研究了一种鲁棒随机保结构拉格朗日数值格式在计算随机微分方程(SDEs)模拟的时变混沌流的有效扩散率时的收敛性分析。我们的数值方案是基于分裂方法来求解相应的SDEs,其中确定性子问题使用结构保持格式离散,随机子问题使用Euler-Maruyama格式离散。对于所提出的数值格式,我们得到了一个清晰的和一致的时间收敛分析,使我们能够准确地计算出SDEs的长期解。因此,我们可以计算时变混沌流的有效扩散系数。最后,我们给出了数值结果,证明了该方法在计算三维空间中随时间变化的Arnold-Beltrami-Childress (ABC)流和Kolmogorov流的有效扩散系数时的准确性和有效性。
{"title":"Computing effective diffusivities in 3D time-dependent chaotic flows with a convergent Lagrangian numerical method","authors":"Zhongjian Wang, J. Xin, Zhiwen Zhang","doi":"10.1051/m2an/2022049","DOIUrl":"https://doi.org/10.1051/m2an/2022049","url":null,"abstract":"In this paper, we study the convergence analysis for a robust stochastic structure-preserving Lagrangian numerical scheme in computing effective diffusivity of time-dependent chaotic flows, which are modeled by stochastic differential equations (SDEs). Our numerical scheme is based on a splitting method to solve the corresponding SDEs in which the deterministic subproblem is discretized using structure-preserving schemes while the random subproblem is discretized using the Euler-Maruyama scheme. We obtain a sharp and uniform-in-time convergence analysis for the proposed numerical scheme that allows us to accurately compute long-time solutions of the SDEs. As such, we can compute the effective diffusivity for time-dependent chaotic flows. Finally, we present numerical results to demonstrate the accuracy and efficiency of the proposed method in computing effective diffusivity for the time-dependent Arnold-Beltrami-Childress (ABC) flow and Kolmogorov flow in three-dimensional space.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91198363","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 3
Truncation errors and modified equations for the lattice Boltzmann method via the corresponding Finite Difference schemes 晶格玻尔兹曼方法的截断误差和相应有限差分格式的修正方程
IF 1.9 3区 数学 Q2 Mathematics Pub Date : 2022-05-05 DOI: 10.1051/m2an/2023008
T. Bellotti
Lattice Boltzmann schemes are efficient numerical methods to solve a broad range of problems under the form of conservation laws. However, they suffer from a chronic lack of clear theoretical foundations. In particular, the consistency analysis and the derivation of the modified equations are still open issues. This has prevented, until today, to have an analogous of the Lax equivalence theorem for Lattice Boltzmann schemes. We propose a rigorous consistency study and the derivation of the modified equations for any lattice Boltzmann scheme under acoustic and diffusive scalings. This is done by passing from a kinetic (lattice Boltzmann) to a macroscopic (Finite Difference) point of view at a fully discrete level in order to eliminate the non-conserved moments relaxing away from the equilibrium. We rewrite the lattice Boltzmann scheme as a multi-step Finite Difference scheme on the conserved variables, as introduced in our previous contribution. We then perform the usual analyses for Finite Difference by exploiting its precise characterization using  matrices of Finite Difference operators. Though we present the derivation of the modified equations until second-order under acoustic scaling, we provide all the elements to extend it to higher orders, since the kinetic-macroscopic connection is conducted at the fully discrete level. Finally, we show that our strategy yields, in a more rigorous setting, the same results as previous works in the literature.
晶格玻尔兹曼格式是一种有效的数值方法,可以在守恒律的形式下解决广泛的问题。然而,它们长期缺乏明确的理论基础。特别是,修正方程的一致性分析和推导仍然是一个有待解决的问题。直到今天,这都阻碍了晶格玻尔兹曼格式的Lax等价定理的类比。本文对任意晶格玻尔兹曼格式在声学标度和扩散标度下的一致性进行了严格的研究,并推导了修正方程。这是通过在完全离散的水平上从动力学(晶格玻尔兹曼)过渡到宏观(有限差分)的观点来完成的,以便消除从平衡状态放松的非守恒力矩。我们将晶格玻尔兹曼格式重写为守恒变量上的多步有限差分格式,如我们之前的贡献所介绍的那样。然后,我们通过利用有限差分算子矩阵的精确表征,对有限差分进行通常的分析。虽然我们提出了在声学尺度下直到二阶的修正方程的推导,但我们提供了将其扩展到更高阶的所有元素,因为动力学-宏观连接是在完全离散的水平上进行的。最后,我们表明,在更严格的设置下,我们的策略产生了与文献中先前工作相同的结果。
{"title":"Truncation errors and modified equations for the lattice Boltzmann method via the corresponding Finite Difference schemes","authors":"T. Bellotti","doi":"10.1051/m2an/2023008","DOIUrl":"https://doi.org/10.1051/m2an/2023008","url":null,"abstract":"Lattice Boltzmann schemes are efficient numerical methods to solve a broad range of problems under the form of conservation laws. However, they suffer from a chronic lack of clear theoretical foundations. In particular, the consistency analysis and the derivation of the modified equations are still open issues. This has prevented, until today, to have an analogous of the Lax equivalence theorem for Lattice Boltzmann schemes. We propose a rigorous consistency study and the derivation of the modified equations for any lattice Boltzmann scheme under acoustic and diffusive scalings. This is done by passing from a kinetic (lattice Boltzmann) to a macroscopic (Finite Difference) point of view at a fully discrete level in order to eliminate the non-conserved moments relaxing away from the equilibrium. We rewrite the lattice Boltzmann scheme as a multi-step Finite Difference scheme on the conserved variables, as introduced in our previous contribution. We then perform the usual analyses for Finite Difference by exploiting its precise characterization using  matrices of Finite Difference operators. Though we present the derivation of the modified equations until second-order under acoustic scaling, we provide all the elements to extend it to higher orders, since the kinetic-macroscopic connection is conducted at the fully discrete level. Finally, we show that our strategy yields, in a more rigorous setting, the same results as previous works in the literature.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-05-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81522390","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 2
A counterexample to analyticity in frictional dynamics 摩擦动力学分析的一个反例
IF 1.9 3区 数学 Q2 Mathematics Pub Date : 2022-04-18 DOI: 10.1051/m2an/2022033
Christopher Roger Dance
We consider the motion of a particle acted on by dry friction and a force that is an analytic function of time. We give a counterexample to the claim that such motions are given by analytic functions of time. Several published arguments concerning existence and uniqueness in unilateral dynamics with friction rely on the analyticity of such motions. The counterexample invalidates those arguments for motions in three or more dimensions.
我们考虑一个质点在干摩擦作用下的运动和一个力是时间的解析函数。对于这种运动是由时间解析函数给出的说法,我们给出一个反例。一些已发表的关于带有摩擦的单边动力学的存在性和唯一性的论点依赖于这种运动的分析性。反例在三维或多维运动中使这些论点无效。
{"title":"A counterexample to analyticity in frictional dynamics","authors":"Christopher Roger Dance","doi":"10.1051/m2an/2022033","DOIUrl":"https://doi.org/10.1051/m2an/2022033","url":null,"abstract":"We consider the motion of a particle acted on by dry friction and a force that is an analytic function of time. We give a counterexample to the claim that such motions are given by analytic functions of time. Several published arguments concerning existence and uniqueness in unilateral dynamics with friction rely on the analyticity of such motions. The counterexample invalidates those arguments for motions in three or more dimensions.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-04-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"89022148","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Optimal error estimates to smooth solutions of the central discontinuous Galerkin methods for nonlinear scalar conservation laws 非线性标量守恒律中心不连续伽辽金方法光滑解的最优误差估计
IF 1.9 3区 数学 Q2 Mathematics Pub Date : 2022-04-15 DOI: 10.1051/m2an/2022037
Mengjiao Jiao, Yan Jiang, Chi-Wang Shu, Mengping Zhang
In this paper, we study the error estimates to sufficiently smooth solutions of the nonlinear scalar conservation laws for the semi-discrete central discontinuous Galerkin (DG) nite element methods on uniform Cartesian meshes. A general approach with an explicitly checkable condition is established for the proof of optimal L2 error estimates of the semi-discrete CDG schemes, and this condition is checked to be valid in one and two dimensions for polynomials of degree up to k = 8. Numerical experiments are given to verify the theoretical results.
本文研究了均匀笛卡尔网格上半离散中心不连续伽辽金(DG)有限元法非线性标量守恒律充分光滑解的误差估计。建立了半离散CDG格式的最优L2误差估计的一般证明方法,该方法具有显式可检验条件,并对k = 8次多项式在一维和二维上的有效性进行了检验。数值实验验证了理论结果。
{"title":"Optimal error estimates to smooth solutions of the central discontinuous Galerkin methods for nonlinear scalar conservation laws","authors":"Mengjiao Jiao, Yan Jiang, Chi-Wang Shu, Mengping Zhang","doi":"10.1051/m2an/2022037","DOIUrl":"https://doi.org/10.1051/m2an/2022037","url":null,"abstract":"In this paper, we study the error estimates to sufficiently smooth solutions of the nonlinear scalar conservation laws for the semi-discrete central discontinuous Galerkin (DG) nite element methods on uniform Cartesian meshes. A general approach with an explicitly checkable condition is established for the proof of optimal L2 error estimates of the semi-discrete CDG schemes, and this condition is checked to be valid in one and two dimensions for polynomials of degree up to k = 8. Numerical experiments are given to verify the theoretical results.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78185936","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
An ultraweak space-time variational formulation for the wave equation: Analysis and efficient numerical solution 波动方程的超弱时空变分公式:分析与有效数值解
IF 1.9 3区 数学 Q2 Mathematics Pub Date : 2022-04-11 DOI: 10.1051/m2an/2022035
K. Urban, Julian Henning, D. Palitta, V. Simoncini
We introduce an ultraweak space-time variational formulation for the wave equation, prove its well-posedness (even in the case of minimal regularity) and optimal inf-sup stability.  Then, we introduce a tensor product-style space-time Petrov-Galerkin discretization with optimal discrete inf-sup stability, obtained by a non-standard definition of the trial space. As a consequence, the numerical approximation error is equal to the residual, which is particularly useful for a posteriori error estimation.    For the arising discrete linear systems in space and time, we  introduce efficient numerical solvers that appropriately exploit the equation structure, either at the preconditioning level or in the approximation phase by using a tailored Galerkin projection. This Galerkin method shows competitive behavior concerning wall-clock time, accuracy and memory as compared with a standard time-stepping method in particular in low regularity cases. Numerical experiments with a 3D (in space) wave equation illustrate our findings.
我们引入了波动方程的一个超弱时空变分公式,证明了它的适定性(即使在最小正则性的情况下)和最优支撑稳定性。然后,通过非标准的试验空间定义,引入具有最优离散稳定性的张量积型时空Petrov-Galerkin离散化。因此,数值近似误差等于残差,这对于后验误差估计特别有用。对于空间和时间中出现的离散线性系统,我们引入了有效的数值求解器,通过使用定制的伽辽金投影,在预处理水平或近似阶段适当地利用方程结构。与标准时间步进方法相比,这种Galerkin方法在挂钟时间、精度和内存方面表现出竞争行为,特别是在低规律性情况下。三维(空间)波动方程的数值实验说明了我们的发现。
{"title":"An ultraweak space-time variational formulation for the wave equation: Analysis and efficient numerical solution","authors":"K. Urban, Julian Henning, D. Palitta, V. Simoncini","doi":"10.1051/m2an/2022035","DOIUrl":"https://doi.org/10.1051/m2an/2022035","url":null,"abstract":"We introduce an ultraweak space-time variational formulation for the wave equation, prove its well-posedness (even in the case of minimal regularity) and optimal inf-sup stability.  Then, we introduce a tensor product-style space-time Petrov-Galerkin discretization with optimal discrete inf-sup stability, obtained by a non-standard definition of the trial space. As a consequence, the numerical approximation error is equal to the residual, which is particularly useful for a posteriori error estimation.\u0000    For the arising discrete linear systems in space and time, we  introduce efficient numerical solvers that appropriately exploit the equation structure, either at the preconditioning level or in the approximation phase by using a tailored Galerkin projection. This Galerkin method shows competitive behavior concerning wall-clock time, accuracy and memory as compared with a standard time-stepping method in particular in low regularity cases. Numerical experiments with a 3D (in space) wave equation illustrate our findings.","PeriodicalId":50499,"journal":{"name":"Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2022-04-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75265513","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 15
期刊
Esaim-Mathematical Modelling and Numerical Analysis-Modelisation Mathematique et Analyse Numerique
全部 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