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

A fully-decoupled second-order-in-time and unconditionally energy stable scheme for the Cahn-Hilliard-Navier-Stokes equations with variable density 针对密度可变的卡恩-希利亚德-纳维尔-斯托克斯方程的完全解耦的二阶时间和无条件能量稳定方案

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2025-03-20 DOI: 10.1016/j.jcp.2025.113943

Jinpeng Zhang , Li Luo , Xiaoping Wang

{"title":"A fully-decoupled second-order-in-time and unconditionally energy stable scheme for the Cahn-Hilliard-Navier-Stokes equations with variable density","authors":"Jinpeng Zhang , Li Luo , Xiaoping Wang","doi":"10.1016/j.jcp.2025.113943","DOIUrl":"10.1016/j.jcp.2025.113943","url":null,"abstract":"<div><div>In this paper, we develop a second-order, fully decoupled, and energy-stable numerical scheme for the Cahn-Hilliard-Navier-Stokes model for two phase flow with variable density and viscosity. We propose a new decoupling Constant Scalar Auxiliary Variable (D-CSAV) method which is easy to generalize to schemes with high order accuracy in time. The method is designed using the “zero-energy-contribution” property while maintaining conservative time discretization for the “non-zero-energy-contribution” terms. A new set of scalar auxiliary variables is introduced to develop second-order-in-time, unconditionally energy stable, and decoupling-type numerical schemes. We also introduce a stabilization parameter <em>α</em> to improve the stability of the scheme by slowing down the dynamics of the scalar auxiliary variables. Our algorithm simplifies to solving three independent linear elliptic systems per time step, two of them with constant coefficients. The update of all scalar auxiliary variables is explicit and decoupled from solving the phase field variable and velocity field. We rigorously prove unconditional energy stability of the scheme and perform extensive benchmark simulations to demonstrate accuracy and efficiency of the method.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"532 ","pages":"Article 113943"},"PeriodicalIF":3.8,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143696116","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

High-order inverse Lax-Wendroff procedure for compressible fluid-structure interaction problems

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2025-03-20 DOI: 10.1016/j.jcp.2025.113942

Zepeng Liu , Yan Jiang , Chi-Wang Shu , Mengping Zhang

{"title":"High-order inverse Lax-Wendroff procedure for compressible fluid-structure interaction problems","authors":"Zepeng Liu , Yan Jiang , Chi-Wang Shu , Mengping Zhang","doi":"10.1016/j.jcp.2025.113942","DOIUrl":"10.1016/j.jcp.2025.113942","url":null,"abstract":"<div><div>In this study, we propose a global high-order approach for fluid-structure interaction (FSI) problems involving compressible inviscid flows and deformable elastic solids. A partitioned coupling strategy is employed to solve the fluid and solid equations. The compressible Euler equations in the fluid domain are solved using a high-order finite difference weighted essentially non-oscillatory (WENO) method on fixed Cartesian Eulerian grids. In the solid domain, the linear elastodynamic equations are discretized via the Lagrangian discontinuous Galerkin (DG) finite element method on unstructured meshes. To handle the moving interface between the fluid and solid domains, we develop a high-order treatment derived from the inverse Lax-Wendroff (ILW) boundary scheme. This approach avoids the need for mesh generation and sub-iterations at each time step, simplifying implementation. Furthermore, the specialized interface treatment ensures stability in challenging cases, such as those involving light solids coupled with heavy fluids. Stability analysis for linear systems further demonstrates the robustness of the method. We validate the proposed approach through numerical tests on one- and two-dimensional problems. The results demonstrate that our method could achieve third-order accuracy for smooth solutions, handle shock induced FSI problems without oscillation, and remain stable across a wide range of material parameters.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"531 ","pages":"Article 113942"},"PeriodicalIF":3.8,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143686397","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 accurate immersed boundary method using radial-basis functions for incompressible flows

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2025-03-19 DOI: 10.1016/j.jcp.2025.113928

Hamayun Farooq , Imran Akhtar , Arman Hemmati , Muhammad Saif Ullah Khalid

{"title":"An accurate immersed boundary method using radial-basis functions for incompressible flows","authors":"Hamayun Farooq , Imran Akhtar , Arman Hemmati , Muhammad Saif Ullah Khalid","doi":"10.1016/j.jcp.2025.113928","DOIUrl":"10.1016/j.jcp.2025.113928","url":null,"abstract":"<div><div>In this work, we introduce a parallel computational solver based on the sharp-interface immersed boundary method for simulating three-dimensional incompressible flows over stationary and moving boundaries. Despite the robustness and ease of implementation of conventional body-conformal grid methods, they are limited to relatively simple immersed geometries, leading to challenges in grid generation and quality. Our approach employs a multi-dimensional ghost-cell methodology and radial basis functions interpolation/splines to achieve accurate boundary condition and superior efficiency. We utilize unstructured triangular elements for geometric surface discretization and non-uniform Cartesian grids for constructing flow domains around the immersed boundaries. Furthermore, full parallelization using domain decomposition ensures scalability on distributed memory platforms, facilitated through message-passing interface libraries. Additionally, we introduce a flow smoothing strategy to mitigate spurious pressure oscillations near immersed bodies. Through simulations of two- and three-dimensional fluid-structure interaction problems, we demonstrate the effectiveness, accuracy, and efficiency of our computational solver.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"531 ","pages":"Article 113928"},"PeriodicalIF":3.8,"publicationDate":"2025-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143686454","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

Non-isothermal filtration problem: Two-temperature computational model

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2025-03-19 DOI: 10.1016/j.jcp.2025.113941

Maksim I. Ivanov, Igor A. Kremer, Yuri M. Laevsky

{"title":"Non-isothermal filtration problem: Two-temperature computational model","authors":"Maksim I. Ivanov, Igor A. Kremer, Yuri M. Laevsky","doi":"10.1016/j.jcp.2025.113941","DOIUrl":"10.1016/j.jcp.2025.113941","url":null,"abstract":"<div><div>The article proposes a new computational model of non-isothermal filtration of a two-phase incompressible fluid. From the point of view of applications, we are talking about the displacement of oil by water when hot water enters from an injection well. The specificity of the proposed model is its two-temperature formulation. The two-temperature formulation is understood as thermal heterogeneity, in which at each point of the domain under consideration the temperatures of the two-phase fluid and porous blocks are determined, and the thermal interaction between the two continua is indicated. The mathematical model is presented in a mixed formulation in the form of a system of first-order equations in terms of four scalar functions (liquid pressure, water saturation, liquid and porous medium temperatures) and three vector functions (total liquid velocity and heat conductive fluxes of the liquid and porous medium). The spatial approximation is based on a combination of mixed FEM and centered FVM. The time approximation consists of using an explicit-implicit scheme with upwinding. In particular, IMPES-type method is used for the filtration equations, and the energy equations explicitly consider convective transfer with the choice of time step according to the CFL condition. For heat exchange between the fluid and the porous medium, both explicit and implicit approximations are used. It is shown that the stability condition of the explicit scheme is significantly weaker than the CFL conditions for convective flows in the mass and energy conservation laws at an accuracy coinciding with the accuracy of the implicit scheme. Also, the two-temperature model made it possible to study the role of heat conductive transfer in a liquid.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"531 ","pages":"Article 113941"},"PeriodicalIF":3.8,"publicationDate":"2025-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143686534","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

Constructing stable, high-order finite-difference operators on point clouds over complex geometries

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2025-03-19 DOI: 10.1016/j.jcp.2025.113940

Jason Hicken, Ge Yan , Sharanjeet Kaur

{"title":"Constructing stable, high-order finite-difference operators on point clouds over complex geometries","authors":"Jason Hicken, Ge Yan , Sharanjeet Kaur","doi":"10.1016/j.jcp.2025.113940","DOIUrl":"10.1016/j.jcp.2025.113940","url":null,"abstract":"<div><div>High-order difference operators with the summation-by-parts (SBP) property can be used to build stable discretizations of hyperbolic conservation laws; however, most high-order SBP operators require a conforming, high-order mesh for the domain of interest. To circumvent this requirement, we present an algorithm for building high-order, diagonal-norm, first-derivative SBP operators on point clouds over level-set geometries. The algorithm is <em>not</em> mesh-free, since it uses a Cartesian cut-cell mesh to define the sparsity pattern of the operators and to provide intermediate quadrature rules; however, the mesh is generated automatically and can be discarded once the SBP operators have been constructed. Using this temporary mesh, we construct local, cell-based SBP difference operators that are assembled into global SBP operators. We identify conditions for the existence of a positive-definite diagonal mass matrix, and we compute the diagonal norm by solving a sparse system of linear inequalities using an interior-point algorithm. We also describe an artificial dissipation operator that complements the first-derivative operators when solving hyperbolic problems, although the dissipation is not required for stability. The numerical results confirm the conditions under which a diagonal norm exists and study the distribution of the norm's entries. In addition, the results verify the accuracy and stability of the point-cloud SBP operators using the linear advection equation.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"532 ","pages":"Article 113940"},"PeriodicalIF":3.8,"publicationDate":"2025-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143697991","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 level set immersed finite element method for parabolic problems on surfaces with moving interfaces

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2025-03-18 DOI: 10.1016/j.jcp.2025.113939

Jiaqi Chen , Xufeng Xiao , Xinlong Feng , Dongwoo Sheen

{"title":"A level set immersed finite element method for parabolic problems on surfaces with moving interfaces","authors":"Jiaqi Chen , Xufeng Xiao , Xinlong Feng , Dongwoo Sheen","doi":"10.1016/j.jcp.2025.113939","DOIUrl":"10.1016/j.jcp.2025.113939","url":null,"abstract":"<div><div>This paper addresses the challenge of solving parabolic moving interface problems on surfaces. These problems have diverse applications, including the Stefan problem, solidification of dendrites on solid surfaces, and flow patterns on soap bubbles. The main difficulties lie in accurately discretizing complex surfaces, efficiently processing interface jump conditions, and tracking the moving interface. Existing numerical methods for interface problems on surfaces have limitations, such as handling only homogeneous jump conditions, having first-order accuracy, or requiring body-fitted nodes. To overcome these limitations, this paper proposes a second-order accurate immersed finite element method (IFEM) for solving parabolic moving interface problems on surfaces. The method is extended to handle non-homogeneous flux jump conditions by enriching the basis functions on interface elements. Furthermore, a novel computational framework is proposed by combining the IFEM with the level set method to track the moving interface. This framework simulates the heat conduction process involving moving interfaces in different velocity fields. The innovation of this paper lies in its ability to handle moving interface problems on surfaces with improved accuracy, efficiency, and versatility compared to existing methods. Verified through numerical simulation, the proposed method and computational framework enable the simulation of a wider range of heat conduction with moving interfaces on surfaces.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"531 ","pages":"Article 113939"},"PeriodicalIF":3.8,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143686401","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

Energy and entropy conserving compatible finite elements with upwinding for the thermal shallow water equations

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2025-03-18 DOI: 10.1016/j.jcp.2025.113937

Tamara A. Tambyah , David Lee , Santiago Badia

{"title":"Energy and entropy conserving compatible finite elements with upwinding for the thermal shallow water equations","authors":"Tamara A. Tambyah , David Lee , Santiago Badia","doi":"10.1016/j.jcp.2025.113937","DOIUrl":"10.1016/j.jcp.2025.113937","url":null,"abstract":"<div><div>In this work, we develop a new compatible finite element formulation of the thermal shallow water equations that conserves energy and mathematical entropies given by buoyancy–related quadratic tracer variances. Our approach relies on restating the governing equations to enable discontinuous approximations of thermodynamic variables and a variational continuous time integration. A key novelty is the inclusion of centred and upwinded fluxes. The proposed semi-discrete system conserves discrete entropy for centred fluxes, monotonically damps entropy for upwinded fluxes, and conserves energy. The fully discrete scheme preserves entropy conservation at the continuous level. The ability of a new linearised Jacobian, which accounts for both centred and upwinded fluxes, to capture large variations in buoyancy and simulate thermally unstable flows for long periods of time is demonstrated for two different transient case studies. The first involves a thermogeostrophic instability where including upwinded fluxes is shown to suppress spurious oscillations while successfully conserving energy and monotonically damping entropy. The second is a double vortex where a constrained fully discrete formulation is shown to achieve exact entropy conservation in time.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"531 ","pages":"Article 113937"},"PeriodicalIF":3.8,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143644166","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

Optimal transport-based displacement interpolation with data augmentation for reduced order modeling of nonlinear dynamical systems

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2025-03-18 DOI: 10.1016/j.jcp.2025.113938

Moaad Khamlich , Federico Pichi , Michele Girfoglio , Annalisa Quaini , Gianluigi Rozza

{"title":"Optimal transport-based displacement interpolation with data augmentation for reduced order modeling of nonlinear dynamical systems","authors":"Moaad Khamlich , Federico Pichi , Michele Girfoglio , Annalisa Quaini , Gianluigi Rozza","doi":"10.1016/j.jcp.2025.113938","DOIUrl":"10.1016/j.jcp.2025.113938","url":null,"abstract":"<div><div>We present a novel reduced-order Model (ROM) that leverages optimal transport (OT) theory and displacement interpolation to enhance the representation of nonlinear dynamics in complex systems. While traditional ROM techniques face challenges in this scenario, especially when data (i.e., observational snapshots) are limited, our method addresses these issues by introducing a data augmentation strategy based on OT principles. The proposed framework generates interpolated solutions tracing geodesic paths in the space of probability distributions, enriching the training dataset for the ROM. A key feature of our approach is its ability to provide a continuous representation of the solution's dynamics by exploiting a virtual-to-real time mapping. This enables the reconstruction of solutions at finer temporal scales than those provided by the original data. To further improve prediction accuracy, we employ Gaussian Process Regression to learn the residual and correct the representation between the interpolated snapshots and the physical solution.</div><div>We demonstrate the effectiveness of our methodology with atmospheric mesoscale benchmarks characterized by highly nonlinear, advection-dominated dynamics. Our results show improved accuracy and efficiency in predicting complex system behaviors, indicating the potential of this approach for a wide range of applications in computational physics and engineering.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"531 ","pages":"Article 113938"},"PeriodicalIF":3.8,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143686533","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

Efficient stochastic response analysis of high-dimensional nonlinear systems subject to multiplicative noise via the DR-PDEE 通过 DR-PDEE 对受乘法噪声影响的高维非线性系统进行高效随机响应分析

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2025-03-17 DOI: 10.1016/j.jcp.2025.113929

Jianbing Chen , Tingting Sun , Pol D. Spanos , Jie Li

{"title":"Efficient stochastic response analysis of high-dimensional nonlinear systems subject to multiplicative noise via the DR-PDEE","authors":"Jianbing Chen , Tingting Sun , Pol D. Spanos , Jie Li","doi":"10.1016/j.jcp.2025.113929","DOIUrl":"10.1016/j.jcp.2025.113929","url":null,"abstract":"<div><div>Significant challenges persist for the reliable probabilistic analyses of high-dimensional nonlinear dynamical systems subject to multiplicative white noise, particularly at the level of probability density. To address this issue, an efficient method is proposed by employing the dimension-reduced probability density evolution equation (DR-PDEE) to determine the probability density of the responses in such systems. The starting point of the method is that, in many cases, only a limited number of quantities in a system are of interest. Thus, the corresponding DR-PDEE is a one- or two-dimensional partial differential equation (PDE) that governs the instantaneous probability density function (PDF) of the quantity/response(s) of interest in high-dimensional stochastic dynamical systems. This is with the stipulation that the response meets the path continuity condition, and there is no restriction on the excitations being multiplicative or additive. The intrinsic drift and diffusion functions in the DR-PDEE are conditional expectation functions of these responses of the original high-dimensional systems that can be reliably estimated, where assessing the latter is specific to multiplicative noise problems. Interestingly, for a wide class of systems subject to multiplicative local noise, the intrinsic diffusion functions are analytically determinable. Subsequently, the instantaneous PDF of the quantity of interest can be efficiently obtained by numerically integrating the one- or two-dimensional DR-PDEE. The accuracy and efficiency of the DR-PDEE are verified by several typical nonlinear high-dimensional dynamical systems. Particularly, the DR-PDEE captures accurately the refined traits that are easily overlooked, and the tail range of response PDFs.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"531 ","pages":"Article 113929"},"PeriodicalIF":3.8,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143686398","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-oscillatory entropy stable DG schemes for hyperbolic conservation law

IF 3.8 2区物理与天体物理

Journal of Computational Physics Pub Date : 2025-03-17 DOI: 10.1016/j.jcp.2025.113926

Yuchang Liu , Wei Guo , Yan Jiang , Mengping Zhang

引用次数: 0