首页 > 最新文献

Communications in Computational Physics最新文献

英文 中文
A Causality-DeepONet for Causal Responses of Linear Dynamical Systems 线性动力系统因果响应的因果深网
IF 3.7 3区 物理与天体物理 Q1 PHYSICS, MATHEMATICAL Pub Date : 2024-06-01 DOI: 10.4208/cicp.oa-2023-0078
Lizuo Liu,Kamaljyoti Nath, Wei Cai
In this paper, we propose a DeepONet structure with causality to representcausal linear operators between Banach spaces of time-dependent signals. The theorem of universal approximations to nonlinear operators proposed in [5] is extendedto operators with causalities, and the proposed Causality-DeepONet implements thephysical causality in its framework. The proposed Causality-DeepONet considerscausality (the state of the system at the current time is not affected by that of the future, but only by its current state and past history) and uses a convolution-type weightin its design. To demonstrate its effectiveness in handling the causal response of aphysical system, the Causality-DeepONet is applied to learn the operator representingthe response of a building due to earthquake ground accelerations. Extensive numerical tests and comparisons with some existing variants of DeepONet are carried out,and the Causality-DeepONet clearly shows its unique capability to learn the retardeddynamic responses of the seismic response operator with good accuracy.
在本文中,我们提出了一种具有因果性的 DeepONet 结构,用于表示时间相关信号的巴拿赫空间之间的因果线性算子。本文将[5]中提出的非线性算子通用近似定理扩展到了具有因果性的算子,并在其框架中实现了物理因果性。所提出的 Causality-DeepONet 考虑了因果性(系统当前的状态不受未来状态的影响,只受当前状态和过去历史的影响),并在设计中使用了卷积型权重。为了证明其在处理物理系统因果响应方面的有效性,我们将因果-深度网络用于学习表示建筑物在地震地面加速度作用下的响应的算子。我们进行了广泛的数值测试,并与 DeepONet 的一些现有变体进行了比较,结果清楚地表明,因果-DeepONet 具有独特的能力,能够准确地学习地震响应算子的迟滞动态响应。
{"title":"A Causality-DeepONet for Causal Responses of Linear Dynamical Systems","authors":"Lizuo Liu,Kamaljyoti Nath, Wei Cai","doi":"10.4208/cicp.oa-2023-0078","DOIUrl":"https://doi.org/10.4208/cicp.oa-2023-0078","url":null,"abstract":"In this paper, we propose a DeepONet structure with causality to represent\u0000causal linear operators between Banach spaces of time-dependent signals. The theorem of universal approximations to nonlinear operators proposed in [5] is extended\u0000to operators with causalities, and the proposed Causality-DeepONet implements the\u0000physical causality in its framework. The proposed Causality-DeepONet considers\u0000causality (the state of the system at the current time is not affected by that of the future, but only by its current state and past history) and uses a convolution-type weight\u0000in its design. To demonstrate its effectiveness in handling the causal response of a\u0000physical system, the Causality-DeepONet is applied to learn the operator representing\u0000the response of a building due to earthquake ground accelerations. Extensive numerical tests and comparisons with some existing variants of DeepONet are carried out,\u0000and the Causality-DeepONet clearly shows its unique capability to learn the retarded\u0000dynamic responses of the seismic response operator with good accuracy.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"45 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141502901","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
JefiPIC: A 3-D Full Electromagnetic Particle-in-Cell Simulator Based on Jefimenko’s Equations on GPU JefiPIC:基于杰菲门科方程的 GPU 3-D 全电磁粒子池模拟器
IF 3.7 3区 物理与天体物理 Q1 PHYSICS, MATHEMATICAL Pub Date : 2024-06-01 DOI: 10.4208/cicp.oa-2023-0156
Jian-Nan Chen,Jun-Jie Zhang,Hai-Liang Qiao,Xue-Ming Li, Yong-Tao Zhao
This paper presents a novel 3-D full electromagnetic particle-in-cell (PIC)code called JefiPIC, which uses Jefimenko’s equations as the electromagnetic (EM) fieldsolver through a full-space integration method. Leveraging the power of state-of-the-art graphic processing units (GPUs), we have made the challenging integral task ofPIC simulations achievable. Our proposed code offers several advantages by utilizing the integral method. Firstly, it offers a natural solution for modeling non-neutralplasmas without the need for pre-processing such as solving Poisson’s equation. Secondly, it eliminates the requirement for designing elaborate boundary layers to absorbfields and particles. Thirdly, it maintains the stability of the plasma simulation regardless of the time step chosen. Lastly, it does not require strict charge-conservationparticle-to-grid apportionment techniques or electric field divergence amendment algorithms, which are commonly used in finite-difference time-domain (FDTD)-basedPIC simulations. To validate the accuracy and advantages of our code, we comparedthe evolutions of particles and fields in different plasma systems simulated by threeother codes. Our results demonstrate that the combination of Jefimenko’s equationsand the PIC method can produce accurate particle distributions and EM fields in open-boundary plasma systems. Additionally, our code is able to accomplish these computations within an acceptable execution time. This study highlights the effectivenessand efficiency of JefiPIC, showing its potential for advancing plasma simulations.
本文介绍了一种名为 JefiPIC 的新型三维全电磁粒子入胞(PIC)代码,该代码通过全空间积分方法使用 Jefimenko 方程作为电磁(EM)场求解器。利用最先进的图形处理器(GPU),我们实现了具有挑战性的 PIC 仿真积分任务。通过利用积分法,我们提出的代码具有多项优势。首先,它为非中性等离子体建模提供了一种自然的解决方案,无需进行诸如求解泊松方程等预处理。其次,它无需为吸积场和粒子设计复杂的边界层。第三,无论选择何种时间步长,它都能保持等离子体模拟的稳定性。最后,它不需要基于有限差分时域(FDTD)的 PIC 仿真中常用的严格的电荷保留粒子到网格分摊技术或电场发散修正算法。为了验证我们代码的准确性和优势,我们比较了其他三种代码模拟的不同等离子体系统中粒子和场的演变。结果表明,结合杰菲门科方程和 PIC 方法可以在开放边界等离子体系统中产生精确的粒子分布和电磁场。此外,我们的代码能够在可接受的执行时间内完成这些计算。这项研究突出了 JefiPIC 的有效性和效率,显示了它在推进等离子体模拟方面的潜力。
{"title":"JefiPIC: A 3-D Full Electromagnetic Particle-in-Cell Simulator Based on Jefimenko’s Equations on GPU","authors":"Jian-Nan Chen,Jun-Jie Zhang,Hai-Liang Qiao,Xue-Ming Li, Yong-Tao Zhao","doi":"10.4208/cicp.oa-2023-0156","DOIUrl":"https://doi.org/10.4208/cicp.oa-2023-0156","url":null,"abstract":"This paper presents a novel 3-D full electromagnetic particle-in-cell (PIC)\u0000code called JefiPIC, which uses Jefimenko’s equations as the electromagnetic (EM) field\u0000solver through a full-space integration method. Leveraging the power of state-of-the-art graphic processing units (GPUs), we have made the challenging integral task of\u0000PIC simulations achievable. Our proposed code offers several advantages by utilizing the integral method. Firstly, it offers a natural solution for modeling non-neutral\u0000plasmas without the need for pre-processing such as solving Poisson’s equation. Secondly, it eliminates the requirement for designing elaborate boundary layers to absorb\u0000fields and particles. Thirdly, it maintains the stability of the plasma simulation regardless of the time step chosen. Lastly, it does not require strict charge-conservation\u0000particle-to-grid apportionment techniques or electric field divergence amendment algorithms, which are commonly used in finite-difference time-domain (FDTD)-based\u0000PIC simulations. To validate the accuracy and advantages of our code, we compared\u0000the evolutions of particles and fields in different plasma systems simulated by three\u0000other codes. Our results demonstrate that the combination of Jefimenko’s equations\u0000and the PIC method can produce accurate particle distributions and EM fields in open-boundary plasma systems. Additionally, our code is able to accomplish these computations within an acceptable execution time. This study highlights the effectiveness\u0000and efficiency of JefiPIC, showing its potential for advancing plasma simulations.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"80 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141515596","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
A Model-Data Asymptotic-Preserving Neural Network Method Based on Micro-Macro Decomposition for Gray Radiative Transfer Equations 基于微宏分解的灰色辐射传输方程的模型-数据渐近保留神经网络方法
IF 3.7 3区 物理与天体物理 Q1 PHYSICS, MATHEMATICAL Pub Date : 2024-06-01 DOI: 10.4208/cicp.oa-2022-0315
Hongyan Li,Song Jiang,Wenjun Sun,Liwei Xu, Guanyu Zhou
We propose a model-data asymptotic-preserving neural network(MD-APNN) method to solve the nonlinear gray radiative transfer equations (GRTEs). Thesystem is challenging to be simulated with both the traditional numerical schemesand the vanilla physics-informed neural networks (PINNs) due to the multiscale characteristics. Under the framework of PINNs, we employ a micro-macro decomposition technique to construct a new asymptotic-preserving (AP) loss function, which includes the residual of the governing equations in the micro-macro coupled form, theinitial and boundary conditions with additional diffusion limit information, the conservation laws, and a few labeled data. A convergence analysis is performed for theproposed method, and a number of numerical examples are presented to illustrate theefficiency of MD-APNNs, and particularly, the importance of the AP property in theneural networks for the diffusion dominating problems. The numerical results indicate that MD-APNNs lead to a better performance than APNNs or pure Data-drivennetworks in the simulation of the nonlinear non-stationary GRTEs.
我们提出了一种模型-数据渐近保留神经网络(MD-APNN)方法来求解非线性灰色辐射传递方程(GRTEs)。由于该系统的多尺度特性,用传统数值方案和虚构物理信息神经网络(PINNs)模拟该系统都具有挑战性。在 PINNs 框架下,我们采用微宏分解技术构建了一个新的渐近保全(AP)损失函数,其中包括微宏耦合形式的治理方程残差、初始条件和边界条件以及额外的扩散极限信息、守恒定律和一些标记数据。对所提出的方法进行了收敛分析,并给出了一些数值示例,以说明 MD-APNN 的效率,特别是神经网络的 AP 特性对扩散主导问题的重要性。数值结果表明,在模拟非线性非平稳 GRTEs 时,MD-APNNs 比 APNNs 或纯数据驱动网络具有更好的性能。
{"title":"A Model-Data Asymptotic-Preserving Neural Network Method Based on Micro-Macro Decomposition for Gray Radiative Transfer Equations","authors":"Hongyan Li,Song Jiang,Wenjun Sun,Liwei Xu, Guanyu Zhou","doi":"10.4208/cicp.oa-2022-0315","DOIUrl":"https://doi.org/10.4208/cicp.oa-2022-0315","url":null,"abstract":"We propose a model-data asymptotic-preserving neural network(MD-APNN) method to solve the nonlinear gray radiative transfer equations (GRTEs). The\u0000system is challenging to be simulated with both the traditional numerical schemes\u0000and the vanilla physics-informed neural networks (PINNs) due to the multiscale characteristics. Under the framework of PINNs, we employ a micro-macro decomposition technique to construct a new asymptotic-preserving (AP) loss function, which includes the residual of the governing equations in the micro-macro coupled form, the\u0000initial and boundary conditions with additional diffusion limit information, the conservation laws, and a few labeled data. A convergence analysis is performed for the\u0000proposed method, and a number of numerical examples are presented to illustrate the\u0000efficiency of MD-APNNs, and particularly, the importance of the AP property in the\u0000neural networks for the diffusion dominating problems. The numerical results indicate that MD-APNNs lead to a better performance than APNNs or pure Data-driven\u0000networks in the simulation of the nonlinear non-stationary GRTEs.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"5 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141502900","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
Localization in the Incommensurate Systems: A Plane Wave Study via Effective Potentials 不相称系统中的定位:通过有效势能进行的平面波研究
IF 3.7 3区 物理与天体物理 Q1 PHYSICS, MATHEMATICAL Pub Date : 2024-05-01 DOI: 10.4208/cicp.oa-2023-0203
Ting Wang,Yuzhi Zhou, Aihui Zhou
In this paper, we apply the effective potentials in the localization landscapetheory (Filoche et al., 2012, Arnold et al., 2016) to study the spectral properties of theincommensurate systems. We uniquely develop a plane wave framework for the effective potentials of the incommensurate systems. And utilizing the effective potentialsrepresented by the plane wave, the location of the electron density can be inferred.Moreover, the spectral distribution can be obtained from the effective potential versionof Weyl’s law. We perform some numerical experiments on some typical incommensurate systems, showing that the effective potential provides an alternative tool forinvestigating the localization and spectral properties of the incommensurate systems,without solving the eigenvalue problem explicitly.
在本文中,我们应用局域景观理论(Filoche 等人,2012 年;Arnold 等人,2016 年)中的有效势来研究不相称系统的谱特性。我们独特地建立了一个平面波框架,用于研究不相称系统的有效势。利用平面波所代表的有效势,可以推断出电子密度的位置。此外,还可以从韦尔定律的有效势版本中获得谱分布。我们对一些典型的非共相系统进行了数值实验,结果表明有效电势为研究非共相系统的定位和谱特性提供了另一种工具,而无需明确求解特征值问题。
{"title":"Localization in the Incommensurate Systems: A Plane Wave Study via Effective Potentials","authors":"Ting Wang,Yuzhi Zhou, Aihui Zhou","doi":"10.4208/cicp.oa-2023-0203","DOIUrl":"https://doi.org/10.4208/cicp.oa-2023-0203","url":null,"abstract":"In this paper, we apply the effective potentials in the localization landscape\u0000theory (Filoche et al., 2012, Arnold et al., 2016) to study the spectral properties of the\u0000incommensurate systems. We uniquely develop a plane wave framework for the effective potentials of the incommensurate systems. And utilizing the effective potentials\u0000represented by the plane wave, the location of the electron density can be inferred.\u0000Moreover, the spectral distribution can be obtained from the effective potential version\u0000of Weyl’s law. We perform some numerical experiments on some typical incommensurate systems, showing that the effective potential provides an alternative tool for\u0000investigating the localization and spectral properties of the incommensurate systems,\u0000without solving the eigenvalue problem explicitly.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"29 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140940520","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
Finite-Volume TENO Scheme with a New Cell-Interface Flux Evaluation Strategy for Unstructured Meshes 针对非结构网格采用新的单元界面流量评估策略的有限体积 TENO 方案
IF 3.7 3区 物理与天体物理 Q1 PHYSICS, MATHEMATICAL Pub Date : 2024-05-01 DOI: 10.4208/cicp.oa-2023-0289
Tian Liang, Lin Fu
The development of high-order shock-capturing schemes is critical for compressible fluid simulations, in particular for cases where both shock waves and small-scale turbulence structures present. As one of the state-of-the-art high-order numericalschemes, the family of high-order targeted ENO (TENO) schemes proposed by Fu etal. [Journal of Computational Physics 305 (2016): 333-359] has been demonstrated toperform well for compressible gas dynamics on structured meshes and recently extended to unstructured meshes by Ji et al. [Journal of Scientific Computing 92(2022):1-39]. In this paper, with the observation that the TENO scheme not only providesthe high-order reconstructed data at the cell interface but also features the potential toseparate the local flow scales in the wavenumber space, we propose a low-dissipationfinite-volume TENO scheme with a new cell-interface flux evaluation strategy for unstructured meshes. The novelty originates from the fact that the local flow scales areclassified, following a strong scale separation in the reconstruction process, as “verysmooth” or not. When the corresponding large central-biased stencil for the targetedcell interface is judged to be “very smooth”, a low-dissipation Riemann solver, eventhe non-dissipative central flux scheme, is employed for the cell-interface flux computing. Otherwise, a dissipative approximate Riemann solver is employed to avoid spurious oscillations and achieve stable shock-capturing. Such a strategy provides separatecontrol over the numerical dissipation of the high-order reconstruction process andthe cell-interface flux calculation within a unified framework and leads to a resultantfinite-volume method with extremely low-dissipation properties and good numericalrobustness. Without parameter tuning case by case, a set of canonical benchmark simulations has been conducted to assess the performance of the proposed scheme.
开发高阶冲击捕获方案对于可压缩流体模拟至关重要,尤其是在同时存在冲击波和小尺度湍流结构的情况下。作为最先进的高阶数值方案之一,Fu 等人提出的高阶目标 ENO(TENO)方案系列[《计算物理学报》305 (2016): 3-37] [Journal.Computing Physics 305 (2016: 3-37]] 是目前最先进的高阶数值方案。[Journal of Computational Physics 305 (2016): 333-359]提出的高阶目标ENO(TENO)方案族在结构网格的可压缩气体动力学中表现良好,最近Ji等人[Journal of Scientific Computing 92(2022):1-39]将其扩展到非结构网格。在本文中,我们观察到 TENO 方案不仅能提供单元界面的高阶重建数据,而且还具有在波数空间中分离局部流动尺度的潜力,因此我们提出了一种低耗散有限体积 TENO 方案,并为非结构网格提出了一种新的单元界面通量评估策略。该方案的新颖之处在于,在重建过程中对局部流尺度进行强尺度分离后,将其划分为 "非常平滑 "或 "不平滑"。当目标细胞界面的相应大型中心偏置模版被判定为 "非常平滑 "时,细胞界面流量计算将采用低耗散黎曼求解器,即非耗散中心流量方案。否则,采用耗散近似黎曼求解器,以避免虚假振荡,实现稳定的冲击捕捉。这种策略在一个统一的框架内对高阶重构过程和细胞界面通量计算的数值耗散进行了单独控制,从而产生了一种具有极低耗散特性和良好数值稳健性的有限体积方法。在不逐个调整参数的情况下,进行了一组典型基准模拟,以评估所提出方案的性能。
{"title":"Finite-Volume TENO Scheme with a New Cell-Interface Flux Evaluation Strategy for Unstructured Meshes","authors":"Tian Liang, Lin Fu","doi":"10.4208/cicp.oa-2023-0289","DOIUrl":"https://doi.org/10.4208/cicp.oa-2023-0289","url":null,"abstract":"The development of high-order shock-capturing schemes is critical for compressible fluid simulations, in particular for cases where both shock waves and small-scale turbulence structures present. As one of the state-of-the-art high-order numerical\u0000schemes, the family of high-order targeted ENO (TENO) schemes proposed by Fu et\u0000al. [Journal of Computational Physics 305 (2016): 333-359] has been demonstrated to\u0000perform well for compressible gas dynamics on structured meshes and recently extended to unstructured meshes by Ji et al. [Journal of Scientific Computing 92(2022):\u00001-39]. In this paper, with the observation that the TENO scheme not only provides\u0000the high-order reconstructed data at the cell interface but also features the potential to\u0000separate the local flow scales in the wavenumber space, we propose a low-dissipation\u0000finite-volume TENO scheme with a new cell-interface flux evaluation strategy for unstructured meshes. The novelty originates from the fact that the local flow scales are\u0000classified, following a strong scale separation in the reconstruction process, as “very\u0000smooth” or not. When the corresponding large central-biased stencil for the targeted\u0000cell interface is judged to be “very smooth”, a low-dissipation Riemann solver, even\u0000the non-dissipative central flux scheme, is employed for the cell-interface flux computing. Otherwise, a dissipative approximate Riemann solver is employed to avoid spurious oscillations and achieve stable shock-capturing. Such a strategy provides separate\u0000control over the numerical dissipation of the high-order reconstruction process and\u0000the cell-interface flux calculation within a unified framework and leads to a resultant\u0000finite-volume method with extremely low-dissipation properties and good numerical\u0000robustness. Without parameter tuning case by case, a set of canonical benchmark simulations has been conducted to assess the performance of the proposed scheme.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"191 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140940447","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
A Comparative Study of Hydrodynamic Lattice Boltzmann Equation in Phase-Field-Based Multiphase Flow Models 基于相场的多相流模型中的流体力学晶格玻尔兹曼方程比较研究
IF 3.7 3区 物理与天体物理 Q1 PHYSICS, MATHEMATICAL Pub Date : 2024-05-01 DOI: 10.4208/cicp.re-2023-0066
Qiang He,Yiqian Cheng,Fengming Hu,Weifeng Huang, Decai Li
In recent years, phase-field-based models for multiphase flows have gainedsignificant popularity, particularly within the lattice Boltzmann (LB) community. Thesemodels typically use two lattice Boltzmann equations (LBEs), one for interface trackingand the other for solving hydrodynamic properties. However, for the purposes of thispaper, we focus only on the LB model for hydrodynamics. Our goal is to undertake acomparative investigation into the differences between three classical hydrodynamicLB models proposed by Lee et al. [1], Liang et al. [2] and Fakhari et al. [3]. The interface-tracking equation used in this study is based on the conservative phase-field model.We provide a detailed derivation of the governing equations in each model using theChapman-Enskog analysis. Additionally, three discretization methods for the interaction forces are introduced, and a modified method for the gradient term is proposedbased on the nonequilibrium distribution method. The accuracy of three LB modelsin combination with four discretization methods is examined in this study. Based onthe results, it appears that different combinations of models and methods are appropriate for different types of problems. However, some suggestions for the selection ofhydrodynamic models and discrete methods for the gradient term are provided in thispaper.
近年来,基于相场的多相流模型大受欢迎,尤其是在晶格玻尔兹曼(LB)领域。这些模型通常使用两个晶格玻尔兹曼方程(LBE),一个用于界面跟踪,另一个用于求解流体动力学特性。不过,在本文中,我们只关注流体力学的 LB 模型。我们的目标是比较研究 Lee 等人[1]、Liang 等人[2]和 Fakhari 等人[3]提出的三种经典流体力学 LB 模型之间的差异。本研究中使用的界面跟踪方程是基于保守相场模型的。我们利用查普曼-恩斯科格分析法详细推导了每个模型的控制方程。此外,我们还介绍了三种相互作用力离散化方法,并在非平衡分布法的基础上提出了梯度项的修正方法。本研究考察了三种 LB 模型与四种离散化方法相结合的精度。从结果来看,不同的模型和方法组合适用于不同类型的问题。不过,本文对梯度项流体力学模型和离散方法的选择提出了一些建议。
{"title":"A Comparative Study of Hydrodynamic Lattice Boltzmann Equation in Phase-Field-Based Multiphase Flow Models","authors":"Qiang He,Yiqian Cheng,Fengming Hu,Weifeng Huang, Decai Li","doi":"10.4208/cicp.re-2023-0066","DOIUrl":"https://doi.org/10.4208/cicp.re-2023-0066","url":null,"abstract":"In recent years, phase-field-based models for multiphase flows have gained\u0000significant popularity, particularly within the lattice Boltzmann (LB) community. These\u0000models typically use two lattice Boltzmann equations (LBEs), one for interface tracking\u0000and the other for solving hydrodynamic properties. However, for the purposes of this\u0000paper, we focus only on the LB model for hydrodynamics. Our goal is to undertake a\u0000comparative investigation into the differences between three classical hydrodynamic\u0000LB models proposed by Lee et al. [1], Liang et al. [2] and Fakhari et al. [3]. The interface-tracking equation used in this study is based on the conservative phase-field model.\u0000We provide a detailed derivation of the governing equations in each model using the\u0000Chapman-Enskog analysis. Additionally, three discretization methods for the interaction forces are introduced, and a modified method for the gradient term is proposed\u0000based on the nonequilibrium distribution method. The accuracy of three LB models\u0000in combination with four discretization methods is examined in this study. Based on\u0000the results, it appears that different combinations of models and methods are appropriate for different types of problems. However, some suggestions for the selection of\u0000hydrodynamic models and discrete methods for the gradient term are provided in this\u0000paper.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"3 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140886665","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
The Positivity-Preserving Finite Volume Coupled with Finite Element Method for the Keller-Segel-Navier-Stokes Model 用于凯勒-西格尔-纳维尔-斯托克斯模型的保正有限体积与有限元方法
IF 3.7 3区 物理与天体物理 Q1 PHYSICS, MATHEMATICAL Pub Date : 2024-05-01 DOI: 10.4208/cicp.oa-2023-0309
Ping Zeng, Guanyu Zhou
We propose a linear decoupled positivity-preserving scheme for thechemotaxis-fluid system, which models the interaction between aerobic bacteria andthe fluid flow surrounding them. This scheme comprises the finite element method(FEM) for the fluid equations on a regular triangulation and an upwind finite volumemethod (FVM) for the chemotaxis system on two types of dual mesh. The discretecellular density and chemical concentration are represented as the piecewise constantfunctions on the dual mesh. They can also be equivalently expressed as the piecewise linear functions on the triangulation in the sense of mass-lumping. These discrete solutions are obtained by the upwind finite volume approximation satisfyingthe laws of positivity preservation and mass conservation. The finite element methodis used to compute the numerical velocity in the triangulation, which is then usedto determine the upwind-style numerical flux in the dual mesh. We analyze the $M$-property of the matrices from the discrete system and prove the well-posedness andthe positivity-preserving property. By using the $L^p$-estimate of the discrete Laplaceoperators, semigroup analysis, and induction method, we are able to establish the optimal error estimates for chemical concentration, cellular density, and velocity field in $(l^∞(W^{1,p}), l^∞(L^p),l^∞(W^{1,p}))$-norm. Several numerical examples are presented to verify the theoretical results.
我们提出了一种针对趋化-流体系统的线性解耦正向保留方案,该方案模拟了好氧细菌与其周围流体之间的相互作用。该方案包括在规则三角网上对流体方程采用有限元法(FEM),以及在两种双网格上对趋化系统采用上风有限体积法(FVM)。离散细胞密度和化学浓度在双重网格上表示为片断常数函数。它们也可以等价地表示为质量结块意义上的三角形上的分片线性函数。这些离散解是通过上风有限体积近似得到的,满足正定守恒定律和质量守恒定律。有限元法用于计算三角网中的数值速度,然后用它来确定对偶网格中的上风式数值通量。我们分析了离散系统中矩阵的 $M$ 特性,并证明了其好求性和正保性特性。利用离散拉普拉斯运算器的 $L^p$估计、半群分析和归纳法,我们能够在 $(l^∞(W^{1,p}), l^∞(L^p),l^∞(W^{1,p}))$ 规范下建立化学浓度、细胞密度和速度场的最优误差估计。本文列举了几个数值例子来验证理论结果。
{"title":"The Positivity-Preserving Finite Volume Coupled with Finite Element Method for the Keller-Segel-Navier-Stokes Model","authors":"Ping Zeng, Guanyu Zhou","doi":"10.4208/cicp.oa-2023-0309","DOIUrl":"https://doi.org/10.4208/cicp.oa-2023-0309","url":null,"abstract":"We propose a linear decoupled positivity-preserving scheme for the\u0000chemotaxis-fluid system, which models the interaction between aerobic bacteria and\u0000the fluid flow surrounding them. This scheme comprises the finite element method\u0000(FEM) for the fluid equations on a regular triangulation and an upwind finite volume\u0000method (FVM) for the chemotaxis system on two types of dual mesh. The discrete\u0000cellular density and chemical concentration are represented as the piecewise constant\u0000functions on the dual mesh. They can also be equivalently expressed as the piecewise linear functions on the triangulation in the sense of mass-lumping. These discrete solutions are obtained by the upwind finite volume approximation satisfying\u0000the laws of positivity preservation and mass conservation. The finite element method\u0000is used to compute the numerical velocity in the triangulation, which is then used\u0000to determine the upwind-style numerical flux in the dual mesh. We analyze the $M$-property of the matrices from the discrete system and prove the well-posedness and\u0000the positivity-preserving property. By using the $L^p$-estimate of the discrete Laplace\u0000operators, semigroup analysis, and induction method, we are able to establish the optimal error estimates for chemical concentration, cellular density, and velocity field in $(l^∞(W^{1,p}), l^∞(L^p),l^∞(W^{1,p}))$-norm. Several numerical examples are presented to verify the theoretical results.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"308 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140940615","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
A Characteristic Mapping Method for Vlasov-Poisson with Extreme Resolution Properties 具有极限分辨率特性的 Vlasov-Poisson 特征映射方法
IF 3.7 3区 物理与天体物理 Q1 PHYSICS, MATHEMATICAL Pub Date : 2024-05-01 DOI: 10.4208/cicp.oa-2024-0012
Philipp Krah,Xi-Yuan Yin,Julius Bergmann,Jean-Christophe Nave, Kai Schneider
We propose an efficient semi-Lagrangian characteristic mapping method forsolving the one+one-dimensional Vlasov-Poisson equations with high precision on acoarse grid. The flow map is evolved numerically and exponential resolution in linear time is obtained. Global third-order convergence in space and time is shown andconservation properties are assessed. For benchmarking, we consider linear and nonlinear Landau damping and the two-stream instability. We compare the results with aFourier pseudo-spectral method and results from the literature. The extreme fine-scaleresolution features are illustrated showing the method’s capabilities to efficiently treatfilamentation in fusion plasma simulations.
我们提出了一种高效的半拉格朗日特征映射方法,用于在粗网格上高精度地求解一维+一维 Vlasov-Poisson 方程。流图以数值方式演化,并在线性时间内获得指数分辨率。在空间和时间上显示了全局三阶收敛性,并评估了保存特性。作为基准,我们考虑了线性和非线性朗道阻尼和双流不稳定性。我们将结果与傅里叶伪谱方法和文献结果进行了比较。极端精细分辨率特征说明了该方法在核聚变等离子体模拟中有效处理细丝化的能力。
{"title":"A Characteristic Mapping Method for Vlasov-Poisson with Extreme Resolution Properties","authors":"Philipp Krah,Xi-Yuan Yin,Julius Bergmann,Jean-Christophe Nave, Kai Schneider","doi":"10.4208/cicp.oa-2024-0012","DOIUrl":"https://doi.org/10.4208/cicp.oa-2024-0012","url":null,"abstract":"We propose an efficient semi-Lagrangian characteristic mapping method for\u0000solving the one+one-dimensional Vlasov-Poisson equations with high precision on a\u0000coarse grid. The flow map is evolved numerically and exponential resolution in linear time is obtained. Global third-order convergence in space and time is shown and\u0000conservation properties are assessed. For benchmarking, we consider linear and nonlinear Landau damping and the two-stream instability. We compare the results with a\u0000Fourier pseudo-spectral method and results from the literature. The extreme fine-scale\u0000resolution features are illustrated showing the method’s capabilities to efficiently treat\u0000filamentation in fusion plasma simulations.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"41 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140940617","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
Geometric Decomposition and Efficient Implementation of High Order Face and Edge Elements 高阶面和边元素的几何分解与高效实现
IF 3.7 3区 物理与天体物理 Q1 PHYSICS, MATHEMATICAL Pub Date : 2024-05-01 DOI: 10.4208/cicp.oa-2023-0249
Chunyu Chen,Long Chen,Xuehai Huang, Huayi Wei
This study investigates high-order face and edge elements in finite elementmethods, with a focus on their geometric attributes, indexing management, and practical application. The exposition begins by a geometric decomposition of Lagrangefinite elements, setting the foundation for further analysis. The discussion then extends to $H$(div)-conforming and $H$(curl)-conforming finite element spaces, adoptingvariable frames across differing sub-simplices. The imposition of tangential or normalcontinuity is achieved through the strategic selection of corresponding bases. The paper concludes with a focus on efficient indexing management strategies for degreesof freedom, offering practical guidance to researchers and engineers. It serves as acomprehensive resource that bridges the gap between theory and practice.
本研究探讨了有限元方法中的高阶面元和边元,重点是它们的几何属性、索引管理和实际应用。论述从拉格朗日有限元的几何分解开始,为进一步分析奠定了基础。然后,讨论延伸到$H$(div)-conforming和$H$(curl)-conforming有限元空间,在不同的子简约中采用不同的框架。通过策略性地选择相应的基,可以实现切向或法向连续性。论文最后重点介绍了自由度的高效索引管理策略,为研究人员和工程师提供了实用指导。它是一种综合性资源,在理论与实践之间架起了一座桥梁。
{"title":"Geometric Decomposition and Efficient Implementation of High Order Face and Edge Elements","authors":"Chunyu Chen,Long Chen,Xuehai Huang, Huayi Wei","doi":"10.4208/cicp.oa-2023-0249","DOIUrl":"https://doi.org/10.4208/cicp.oa-2023-0249","url":null,"abstract":"This study investigates high-order face and edge elements in finite element\u0000methods, with a focus on their geometric attributes, indexing management, and practical application. The exposition begins by a geometric decomposition of Lagrange\u0000finite elements, setting the foundation for further analysis. The discussion then extends to $H$(div)-conforming and $H$(curl)-conforming finite element spaces, adopting\u0000variable frames across differing sub-simplices. The imposition of tangential or normal\u0000continuity is achieved through the strategic selection of corresponding bases. The paper concludes with a focus on efficient indexing management strategies for degrees\u0000of freedom, offering practical guidance to researchers and engineers. It serves as a\u0000comprehensive resource that bridges the gap between theory and practice.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"29 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140940521","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
A Data-Driven Scale-Invariant Weighted Compact Nonlinear Scheme for Hyperbolic Conservation Laws 双曲守恒定律的数据驱动尺度不变加权紧凑非线性方案
IF 3.7 3区 物理与天体物理 Q1 PHYSICS, MATHEMATICAL Pub Date : 2024-05-01 DOI: 10.4208/cicp.oa-2023-0162
Zixuan Zhang,Yidao Dong,Yuanyang Zou,Hao Zhang, Xiaogang Deng
With continuous developments in various techniques, machine learning isbecoming increasingly viable and promising in the field of fluid mechanics. In thisarticle, we present a machine learning approach for enhancing the resolution and robustness of the weighted compact nonlinear scheme (WCNS). We employ a neuralnetwork as a weighting function in the WCNS scheme and follow a data-driven approach to train this neural network. Neural networks can learn a new smoothnessmeasure and calculate a weight function inherently. To facilitate the machine learning task and train with fewer data, we integrate the prior knowledge into the learningprocess, such as a Galilean invariant input layer and CNS polynomials. The normalization in the Delta layer (the so-called Delta layer is used to calculate input features)ensures that the WCNS3-NN schemes achieve a scale-invariant property (Si-property)with an arbitrary scale of a function, and an essentially non-oscillatory approximation of a discontinuous function (ENO-property). The Si-property and ENO-propertyof the data-driven WCNS schemes are validated numerically. Several one- and two-dimensional benchmark examples, including strong shocks and shock-density waveinteractions, are presented to demonstrate the advantages of the proposed method.
随着各种技术的不断发展,机器学习在流体力学领域的可行性和前景越来越广阔。在本文中,我们提出了一种机器学习方法,用于提高加权紧凑非线性方案(WCNS)的分辨率和鲁棒性。我们在 WCNS 方案中采用神经网络作为加权函数,并采用数据驱动的方法来训练该神经网络。神经网络可以学习新的平滑度度量,并计算出固有的权重函数。为了简化机器学习任务并使用更少的数据进行训练,我们将先验知识整合到了学习过程中,例如伽利略不变输入层和 CNS 多项式。三角洲层的归一化(所谓的三角洲层用于计算输入特征)确保了 WCNS3-NN 方案在函数的任意标度下实现标度不变特性(Si-property),以及非连续函数的基本非振荡近似(ENO-property)。数据驱动的 WCNS 方案的 Si-property 和 ENO-property 得到了数值验证。介绍了几个一维和二维基准示例,包括强冲击和冲击-密度波相互作用,以展示所提方法的优势。
{"title":"A Data-Driven Scale-Invariant Weighted Compact Nonlinear Scheme for Hyperbolic Conservation Laws","authors":"Zixuan Zhang,Yidao Dong,Yuanyang Zou,Hao Zhang, Xiaogang Deng","doi":"10.4208/cicp.oa-2023-0162","DOIUrl":"https://doi.org/10.4208/cicp.oa-2023-0162","url":null,"abstract":"With continuous developments in various techniques, machine learning is\u0000becoming increasingly viable and promising in the field of fluid mechanics. In this\u0000article, we present a machine learning approach for enhancing the resolution and robustness of the weighted compact nonlinear scheme (WCNS). We employ a neural\u0000network as a weighting function in the WCNS scheme and follow a data-driven approach to train this neural network. Neural networks can learn a new smoothness\u0000measure and calculate a weight function inherently. To facilitate the machine learning task and train with fewer data, we integrate the prior knowledge into the learning\u0000process, such as a Galilean invariant input layer and CNS polynomials. The normalization in the Delta layer (the so-called Delta layer is used to calculate input features)\u0000ensures that the WCNS3-NN schemes achieve a scale-invariant property (Si-property)\u0000with an arbitrary scale of a function, and an essentially non-oscillatory approximation of a discontinuous function (ENO-property). The Si-property and ENO-property\u0000of the data-driven WCNS schemes are validated numerically. Several one- and two-dimensional benchmark examples, including strong shocks and shock-density wave\u0000interactions, are presented to demonstrate the advantages of the proposed method.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"27 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140941880","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
期刊
Communications in Computational Physics
全部 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