Journal of Computational Physics最新文献_第10页

A novel energy-bounded Boussinesq model and a well balanced and stable numerical discretisation 新颖的能量约束布森斯克模型和平衡稳定的数值离散化方法

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-10-18 DOI: 10.1016/j.jcp.2024.113516

Magnus Svärd, Henrik Kalisch

引用次数: 0

Parallel Eulerian-Lagrangian coupling method on hierarchical meshes 分层网格上的并行欧拉-拉格朗日耦合方法

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-10-18 DOI: 10.1016/j.jcp.2024.113509

Tim Wegmann , Ansgar Niemöller , Matthias Meinke , Wolfgang Schröder

{"title":"Parallel Eulerian-Lagrangian coupling method on hierarchical meshes","authors":"Tim Wegmann , Ansgar Niemöller , Matthias Meinke , Wolfgang Schröder","doi":"10.1016/j.jcp.2024.113509","DOIUrl":"10.1016/j.jcp.2024.113509","url":null,"abstract":"<div><div>An Eulerian-Lagrangian coupling method based on hierarchical meshes is presented, which allows an efficient parallelization on high-performance computing hardware. It features an interleaved execution pattern with non-blocking communication, where the hierarchical mesh structure facilitates the redistribution of the computational load. The Lagrangian and Eulerian solvers use hierarchical Cartesian meshes which share a common coarse mesh level. The domain decomposition is based on a space-filling curve defined on the joint computational mesh, where the load is projected to a coarse mesh level used for the partitioning. The performance of the coupled method is evaluated for the problem of spray modeling in turbulent flow. A solution adaptive mesh is utilized for the large-eddy simulation of the flow field and the Lagrangian tracking method is used for the spray particles. Static and dynamic workload estimators are compared with respect to the alleviation of load imbalances. Liquid fuel spray injection in a constant pressure chamber and in an internal combustion engine serves as applications with varying scale resolution and localized computational load. The parallel efficiency of the approach on high performance systems is demonstrated for meshes with up to <span><math><mn>2.8</mn><mo>⋅</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>9</mn></mrow></msup></math></span> cells and <span><math><mn>21</mn><mo>⋅</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>6</mn></mrow></msup></math></span> particles. Detailed performance analyses show a performance gain of the novel algorithm of approx. 20% compared to a non-interleaved time step execution for two-way coupled spray injection simulations. Results of strong scaling experiments at different injection phases show a good parallel performance with an efficiency of up to 81% using 262000 MPI processes.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"521 ","pages":"Article 113509"},"PeriodicalIF":3.8,"publicationDate":"2024-10-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142554352","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

The localized radial basis function collocation method for dendritic solidification, solid phase sintering and wetting phenomenon based on phase field 基于相场的树枝状凝固、固相烧结和润湿现象的局部径向基函数搭配法

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-10-18 DOI: 10.1016/j.jcp.2024.113515

Pengfei Jiang , Hui Zheng , Jingang Xiong , Timon Rabczuk

引用次数: 0

Spurious vorticity in Eulerian and Lagrangian methods 欧拉和拉格朗日方法中的虚假涡度

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-10-18 DOI: 10.1016/j.jcp.2024.113510

David Sidilkover

{"title":"Spurious vorticity in Eulerian and Lagrangian methods","authors":"David Sidilkover","doi":"10.1016/j.jcp.2024.113510","DOIUrl":"10.1016/j.jcp.2024.113510","url":null,"abstract":"<div><div>Lagrangian methods for computational continuum mechanics, since their inception, traditionally relied on staggered meshes. This feature, while facilitating their robustness and reliability, presented some difficulties. The latter motivated the search for collocated Lagrangian schemes. One of the attempts to develop such a scheme was the CAVEAT method/code. Numerical solutions produced by this method suffered sometimes from large vorticity errors, which could lead to mesh entanglement and premature run termination. The efforts to devise a more robust collocated scheme began to bear fruit a couple of decades later starting from the groundbreaking method GLACE, closely followed by EUCCLHYD and later on by CCH and others.</div><div>One of the aims of this paper is to present a novel Lagrangian collocated <em>factorizable</em> scheme. The notion of a <em>factorizable</em> method was introduced more than two decades ago within the Eulerian approach. It designates a numerical scheme that reflects/preserves the mixed character of the Euler equations, i.e. does not introduce non-physical coupling between the different factors of the system of equations - advection and acoustics operators.</div><div>Another aim of this paper is to explore the connection between the <em>factorizability</em> property of a Lagrangian method and whether or not it suffers from spurious vorticity. Several existing schemes are surveyed for this purpose. A conjecture summarizing our findings is formulated.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"520 ","pages":"Article 113510"},"PeriodicalIF":3.8,"publicationDate":"2024-10-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535509","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

An interface tracking method with triangle edge cuts 使用三角形边缘切割的界面跟踪方法

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-10-17 DOI: 10.1016/j.jcp.2024.113504

Mengdi Wang , Matthew Cong , Bo Zhu

引用次数: 0

Energy bounds for discontinuous Galerkin spectral element approximations of well-posed overset grid problems for hyperbolic systems 非连续伽勒金谱元近似双曲系统的好求过载网格问题的能量边界

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-10-17 DOI: 10.1016/j.jcp.2024.113508

David A. Kopriva , Andrew R. Winters , Jan Nordström

引用次数: 0

Bound-preserving OEDG schemes for Aw–Rascle–Zhang traffic models on networks 网络上 Aw-Rascle-Zhang 流量模型的保界 OEDG 方案

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-10-16 DOI: 10.1016/j.jcp.2024.113507

Wei Chen , Shumo Cui , Kailiang Wu , Tao Xiong

{"title":"Bound-preserving OEDG schemes for Aw–Rascle–Zhang traffic models on networks","authors":"Wei Chen , Shumo Cui , Kailiang Wu , Tao Xiong","doi":"10.1016/j.jcp.2024.113507","DOIUrl":"10.1016/j.jcp.2024.113507","url":null,"abstract":"<div><div>Physical solutions to the widely used Aw–Rascle–Zhang (ARZ) traffic model and the adapted pressure ARZ model should satisfy the positivity of density, the minimum and maximum principles with respect to the velocity <em>v</em> and other Riemann invariants. Many numerical schemes suffer from instabilities caused by violating these bounds, and the only existing bound-preserving (BP) numerical scheme (for ARZ model) is random, only first-order accurate, and not strictly conservative. This paper introduces arbitrarily high-order provably BP discontinuous Galerkin (DG) schemes for these two models, preserving all the aforementioned bounds except the maximum principle of <em>v</em>, which has been rigorously proven to conflict with the consistency and conservation of numerical schemes. Although the maximum principle of <em>v</em> is not directly enforced, we find that the strictly preserved maximum principle of another Riemann invariant <em>w</em> actually enforces an alternative upper bound on <em>v</em>. At the core of this work, analyzing and rigorously proving the BP property is a particularly nontrivial task: the Lax–Friedrichs (LF) splitting property, usually expected for hyperbolic conservation laws and employed to construct BP schemes, does not hold for these two models. To overcome this challenge, we formulate a generalized version of the LF splitting property, and prove it via the geometric quasilinearization approach (Wu and Shu, 2023 <span><span>[47]</span></span>). To suppress spurious oscillations in the DG solutions, we incorporate the oscillation-eliminating technique, recently proposed in (Peng et al., 2024 <span><span>[34]</span></span>), which is based on the solution operator of a novel damping equation. Several numerical examples are included to demonstrate the effectiveness, accuracy, and BP properties of our schemes, with applications to traffic simulations on road networks.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"520 ","pages":"Article 113507"},"PeriodicalIF":3.8,"publicationDate":"2024-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535524","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Non-radiating sources in elastodynamics and their applications in the exterior cloaking 弹性动力学中的非辐射源及其在外部隐形中的应用

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-10-16 DOI: 10.1016/j.jcp.2024.113505

Shuxiang Chen , Jue Wang , Lei Zhang

{"title":"Non-radiating sources in elastodynamics and their applications in the exterior cloaking","authors":"Shuxiang Chen , Jue Wang , Lei Zhang","doi":"10.1016/j.jcp.2024.113505","DOIUrl":"10.1016/j.jcp.2024.113505","url":null,"abstract":"<div><div>In this work, we develop a mathematical framework on constructed general non-radiating sources of elastic waves governed by the Navier equation via the approach of Helmholtz decomposition and potential theory in elastodynamics. Our study offers a rather comprehensive analysis. We first provide a rigorous justification of the general non-radiating sources. Based on the complete destructive interference of external elastic fields generated by specific radiating sources, a general non-radiating elastic source is derived and shown to possess a hidden interior wave field. For an incident wave, targets remain invisible within non-radiating source regions, and the geometry and boundary conditions of obstacles can be very general, which holds significant practical implications. Moreover, we introduce an effective novel method for designing such generalized non-radiating sources. To avoid the complex structure, we propose to use radiating source overlay construction on specific nodes at the boundary of non-radiating regions construction and derive sharp error estimates to evaluate the cloaking performance. The proposed scheme is capable of nearly cloaking arbitrary obstacles with a high accuracy. Numerical verifications validate the precision of our analytical findings.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"520 ","pages":"Article 113505"},"PeriodicalIF":3.8,"publicationDate":"2024-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142441918","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Preconditioning elliptic operators in high-performance all-scale atmospheric models on unstructured meshes 在非结构网格上对高性能全尺度大气模型中的椭圆算子进行预处理

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-10-15 DOI: 10.1016/j.jcp.2024.113503

Mike Gillard , Joanna Szmelter , Francesco Cocetta

{"title":"Preconditioning elliptic operators in high-performance all-scale atmospheric models on unstructured meshes","authors":"Mike Gillard , Joanna Szmelter , Francesco Cocetta","doi":"10.1016/j.jcp.2024.113503","DOIUrl":"10.1016/j.jcp.2024.113503","url":null,"abstract":"<div><div>Effective simulation of all-scale atmospheric flows – e.g., cloud-resolving global weather – involves semi-implicit integration of the non-hydrostatic compressible Euler equations under gravity on a rotating sphere. Such integrations depend on complex non-symmetric elliptic solvers. The condition number of the underlying sparse linear operator is <span><math><mi>O</mi><mo>(</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>10</mn></mrow></msup><mo>)</mo></math></span>, which necessitates bespoke operator preconditioning. This paper highlights the development and implementation on unstructured meshes of specialised preconditioners for the non-symmetric Krylov-subspace solver. These developments are set in the context of a massively-parallel high-performance computing environment, aimed at architectures evolving towards exascale.</div><div>The baroclinic instability benchmark bearing representative features relevant to numerical weather prediction (NWP) has been selected to study the performance of the preconditioning options. The reported results illustrate the improved performance with the new preconditioning options. In particular, the Jacobi based option, for the computational meshes tested in this study, provides an excellent time to solution improvement.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"520 ","pages":"Article 113503"},"PeriodicalIF":3.8,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535523","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Structural topology optimization based on deep learning 基于深度学习的结构拓扑优化

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-10-15 DOI: 10.1016/j.jcp.2024.113506

Yingning Gao, Sizhu Zhou, Meiqiu Li

{"title":"Structural topology optimization based on deep learning","authors":"Yingning Gao, Sizhu Zhou, Meiqiu Li","doi":"10.1016/j.jcp.2024.113506","DOIUrl":"10.1016/j.jcp.2024.113506","url":null,"abstract":"<div><div>In the mechanical design of structures, traditional topology optimization methods involve numerous finite element iterative analyses, leading to a significant expenditure of computational resources. Therefore, the improved multi-scale gradient generative adversarial networks topology optimization technique is proposed. The topology optimization condition parameters are compressed into a low-dimensional latent space feature representation using the encoder, allowing the model to better extract features from these parameters. To speed up model training, the generator and discriminator networks use lightweight residual convolutional blocks. The hybrid attention mechanism extracts prominent region features from the topology optimization structure map. The model training process is guided by a multi-dimensional fusion loss function to enhance the quality of generated model samples. Finally, transferring the parameters of the low-resolution topology optimization model to the high-resolution model enables complete training on a limited amount of high-resolution topology optimization datasets. The experimental data on the low- and high-resolution topology optimization datasets demonstrate that, when compared to alternative methods, this method produces better-quality topology optimization structure maps. Additionally, it can generate high-resolution topology optimization structure maps in minimal time, enabling real-time topology optimization.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"520 ","pages":"Article 113506"},"PeriodicalIF":3.8,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142442033","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0