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

A simplified fast multipole method based on strong recursive skeletonization

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-12-30 DOI: 10.1016/j.jcp.2024.113707

Anna Yesypenko , Chao Chen , Per-Gunnar Martinsson

{"title":"A simplified fast multipole method based on strong recursive skeletonization","authors":"Anna Yesypenko , Chao Chen , Per-Gunnar Martinsson","doi":"10.1016/j.jcp.2024.113707","DOIUrl":"10.1016/j.jcp.2024.113707","url":null,"abstract":"<div><div>This work introduces a kernel-independent, multilevel, adaptive algorithm for efficiently evaluating a discrete convolution kernel with a given source distribution. The method is based on linear algebraic tools such as low rank approximation and “skeleton representations” to approximate far-field interactions. While this work is related to previous linear algebraic formulations of the fast multipole method, the proposed algorithm is distinguished by relying on simpler data structures.</div><div>The proposed algorithm eliminates the need for explicit interaction lists by restructuring computations to operate exclusively on the near-neighbor list at each level of the tree, thereby simplifying both implementation and data structures. This work also introduces novel translation operators that significantly simplify the handling of adaptive point distributions. As a kernel-independent approach, it only requires evaluation of the kernel function, making it easily adaptable to a variety of kernels. By using operations on the neighbor list (of size at most 27 in 3D) rather than the interaction list (of size up to 189 in 3D), the algorithm is particularly well-suited for parallel implementation on modern hardware.</div><div>Numerical experiments on uniform and non-uniform point distributions in 2D and 3D demonstrate the effectiveness of the proposed parallel algorithm for Laplace and (low-frequency) Helmholtz kernels. The algorithm constructs a tailored skeleton representation for the given geometry during a precomputation stage. After precomputation, the fast summation achieves high efficiency on the GPU using batched linear algebra operations.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"524 ","pages":"Article 113707"},"PeriodicalIF":3.8,"publicationDate":"2024-12-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143131708","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

Physics-informed data-driven cavitation model for a specific Mie–Grüneisen equation of state

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-12-30 DOI: 10.1016/j.jcp.2024.113703

Minsheng Huang , Chengbao Yao , Pan Wang , Lidong Cheng , Wenjun Ying

{"title":"Physics-informed data-driven cavitation model for a specific Mie–Grüneisen equation of state","authors":"Minsheng Huang , Chengbao Yao , Pan Wang , Lidong Cheng , Wenjun Ying","doi":"10.1016/j.jcp.2024.113703","DOIUrl":"10.1016/j.jcp.2024.113703","url":null,"abstract":"<div><div>Unsteady cavitation, as observed in phenomena like underwater explosions, entails dynamically evolving boundaries and developing dimensions of cavitation before collapse. The classic one-fluid models often fail to accurately simulate this, as they either rely solely on physics or experimental data. In this study, we introduce a novel one-fluid cavitation model tailored for a specific Mie-Grüneisen equation of state (EOS) known as polynomial EOS, employing an artificial neural network. Our approach integrates physics-informed equations with experimental data through an optimization problem. This physics-informed data-driven model accurately predicts pressure within the cavitation region, where pressure tends towards zero along with density. We apply this model to complex compressible multi-phase flow simulations, including nuclear and underwater explosions. Validation against experimental data and comparison with existing models precedes its application in one- and two-dimensional cases. The simulation of three-dimensional cavitation phenomena near submarines during underwater explosions yields a high-quality numerical solution. The outcome underscores the significant engineering value of the novel cavitation model.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"524 ","pages":"Article 113703"},"PeriodicalIF":3.8,"publicationDate":"2024-12-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143131706","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

A new flow dynamic approach for Wasserstein gradient flows

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-12-27 DOI: 10.1016/j.jcp.2024.113696

Qing Cheng , Qianqian Liu , Wenbin Chen , Jie Shen

引用次数: 0

Diffuse-interface modeling and energy-stable numerical framework for the heat transfer-coupled two-phase fluids in contact with solids

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-12-27 DOI: 10.1016/j.jcp.2024.113699

Fang Zhu , Keyue Sun , Guangtao Zhang , Junxiang Yang

{"title":"Diffuse-interface modeling and energy-stable numerical framework for the heat transfer-coupled two-phase fluids in contact with solids","authors":"Fang Zhu , Keyue Sun , Guangtao Zhang , Junxiang Yang","doi":"10.1016/j.jcp.2024.113699","DOIUrl":"10.1016/j.jcp.2024.113699","url":null,"abstract":"<div><div>To efficiently simulate the heat transfer-coupled two-phase fluid flows with wetting condition in irregular domains, we develop a diffuse-interface heat fluid system. A traditional ternary Cahn–Hilliard model is modified to approximate the damping of solid on fluid. Based on equilibrium interface assumption and Young's equality, an extra term reflecting wetting contact line is derived and added into the diffuse-interface model. The heat transfer and fluid dynamics are described by coupling the penalized incompressible Navier–Stokes equations and a diffusion conduction equation with variable coefficients. The proposed model can efficiently describe complex heat fluid flows in contact with solids because the computations are implemented in regular rectangular domains. The complex techniques for the treatment of fluid-solid boundary are not necessary. Moreover, the proposed model also leads to an energy dissipation law. To satisfy this basic physical property in simulation, we propose linear, totally decoupled, and second-order energy-stable scheme to update the solutions. The time-discretized energy law is analytically estimated. In each time step, the solutions can be easily obtained by solving several linear elliptic-type equations in a step-by-step manner. Extensive numerical experiments in two- and three-dimensional spaces are implemented to validate the accuracy and stability of our method. These results also indicate that the proposed method has good potential in simulating complex fluid flows with heat transfer.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"524 ","pages":"Article 113699"},"PeriodicalIF":3.8,"publicationDate":"2024-12-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143132674","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

Improved gas kinetic flux solver with locally rescaled Gauss-Hermite quadrature for supersonic continuum and rarefied flows

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-12-27 DOI: 10.1016/j.jcp.2024.113700

Zhe Li

引用次数: 0

RF-PINNs: Reactive flow physics-informed neural networks for field reconstruction of laminar and turbulent flames using sparse data

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-12-27 DOI: 10.1016/j.jcp.2024.113698

Vikas Yadav, Mario Casel, Abdulla Ghani

{"title":"RF-PINNs: Reactive flow physics-informed neural networks for field reconstruction of laminar and turbulent flames using sparse data","authors":"Vikas Yadav, Mario Casel, Abdulla Ghani","doi":"10.1016/j.jcp.2024.113698","DOIUrl":"10.1016/j.jcp.2024.113698","url":null,"abstract":"<div><div>Physics-Informed Neural Networks (PINNs) have emerged as a promising tool to model flow fields by embedding physical laws into neural networks and thereby reducing the dependency on data. While PINNs have shown substantial success in modeling non-reactive flows, their application for chemically reacting flows is less well understood. The overarching objective of this work is to propose reactive flow physics-informed neural networks (RF-PINNs) that reconstruct all flow quantities of interest solely based on sparse velocity profiles. We present RF-PINNs to reconstruct steady-state mean fields of premixed laminar and turbulent flames with sparse data. First, we extensively study a prototype flame to determine the best combination of governing equations to find a balance between prediction accuracy and optimization speed. We then demonstrate the RF-PINN prediction capabilities by utilizing measured velocity profiles to reconstruct the entire velocity, temperature, and density mean fields of laminar and turbulent flames. To highlight the bidirectional reconstruction capability, we utilize in the second step measured temperature profiles for both flames and, again, successfully reconstruct the flow fields of interest. For both scenarios, the subsequent comparison of predicted field quantities based on only 0.2% of the entire data set is sufficient for accurate field reconstruction. This study underscores the potential of RF-PINNs to complement comprehensive field data of laminar and turbulent reacting flows using both sparse and noisy experimental data.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"524 ","pages":"Article 113698"},"PeriodicalIF":3.8,"publicationDate":"2024-12-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143132671","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

Approximately pressure-equilibrium-preserving scheme for fully conservative simulations of compressible multi-species and real-fluid interfacial flows

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-12-27 DOI: 10.1016/j.jcp.2024.113701

H. Terashima , N. Ly , M. Ihme

引用次数: 0

Structure preserving hybrid Finite Volume Finite Element method for compressible MHD

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-12-20 DOI: 10.1016/j.jcp.2024.113691

Francesco Fambri , Eric Sonnendrücker

{"title":"Structure preserving hybrid Finite Volume Finite Element method for compressible MHD","authors":"Francesco Fambri , Eric Sonnendrücker","doi":"10.1016/j.jcp.2024.113691","DOIUrl":"10.1016/j.jcp.2024.113691","url":null,"abstract":"<div><div>In this manuscript we present a novel and efficient numerical method for the compressible viscous and resistive MHD equations for all Mach number regimes. The time-integration strategy is a semi-implicit splitting, combined with a hybrid Finite Volume and Finite Element (FE) discretization in space. The nonlinear convection is solved by a robust explicit FV scheme, while the magneto-acoustic terms are treated implicitly in time. The resulting CFL stability condition depends only on the fluid velocity, and not on the Alfvénic and acoustic modes. The magneto-acoustic terms are discretized by compatible FE based on a continuous and a discrete de Rham complexes designed using Finite Element Exterior Calculus (FEEC). Thanks to the use of FEEC, energy stability, magnetic-helicity conservation and the divergence-free conditions can be preserved also at the discrete level. A very efficient splitting approach is used to separate the acoustic and the Alfvénic modes in such a fashion that the original symmetries of the PDE governing equations are preserved. In this way, the algorithm relies on the solution of linear, symmetric and positive-definite algebraic systems, that are very efficiently handled by the simple matrix-free conjugate-gradient method. The resulting algorithm showed to be robust and accurate in low and high Mach regimes even at large Courant numbers. Non-trivial tests are solved in one-, two- and three-space dimensions to confirm the robustness, accuracy, and the low-dissipative and conserving properties of the final algorithm. While the formulation of the method is very general, numerical results for a second-order accurate FV-FE scheme will be presented.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"523 ","pages":"Article 113691"},"PeriodicalIF":3.8,"publicationDate":"2024-12-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143149916","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

An immersed boundary method using online sequential data assimilation

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-12-20 DOI: 10.1016/j.jcp.2024.113697

Miguel M. Valero, Marcello Meldi

{"title":"An immersed boundary method using online sequential data assimilation","authors":"Miguel M. Valero, Marcello Meldi","doi":"10.1016/j.jcp.2024.113697","DOIUrl":"10.1016/j.jcp.2024.113697","url":null,"abstract":"<div><div>A data-driven strategy relying on the Ensemble Kalman Filter (EnKF) is here used to augment the accuracy of a continuous Immersed Boundary Method (IBM). The latter is a classical penalty method accounting for the presence of the immersed body via a volume source term, which is included in the Navier–Stokes equations. The model coefficients of the penalisation method, which are usually selected by the user, are optimised using the data-driven strategy. The parametric inference is governed by the physical knowledge of local and global features of the flow, such as the no-slip condition and the shear stress at the wall. The C++ library CONES (Coupling OpenFOAM with Numerical EnvironmentS) developed by the team is used to perform an online investigation, coupling on-the-fly data from synthetic sensors with results from an ensemble of coarse-grained numerical simulations. The analysis is performed for a classical test case, namely the turbulent plane channel flow with <span><math><mi>R</mi><msub><mrow><mi>e</mi></mrow><mrow><mi>τ</mi></mrow></msub><mo>=</mo><mn>550</mn></math></span>. The results, which are compared with high-fidelity Direct Numerical Simulation (DNS), show that the data-driven procedure exhibits remarkable accuracy despite the ensemble members' relatively low grid resolution. The findings present open perspectives of application in dynamic complex systems, such as digital twins.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"524 ","pages":"Article 113697"},"PeriodicalIF":3.8,"publicationDate":"2024-12-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143132666","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

Novel local characteristic decomposition based path-conservative central-upwind schemes

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2024-12-20 DOI: 10.1016/j.jcp.2024.113692

Shaoshuai Chu , Michael Herty , Alexander Kurganov

引用次数: 0