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

Comparison of rigorous scattering models to accurately replicate the behaviour of scattered electromagnetic waves in optical surface metrology 比较严格的散射模型，以准确复制光学表面计量中的散射电磁波行为

IF 3.8 2区物理与天体物理

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

Helia Hooshmand , Tobias Pahl , Poul-Erik Hansen , Liwei Fu , Alexander Birk , Mirza Karamehmedović , Peter Lehmann , Stephan Reichelt , Richard Leach , Samanta Piano

{"title":"Comparison of rigorous scattering models to accurately replicate the behaviour of scattered electromagnetic waves in optical surface metrology","authors":"Helia Hooshmand , Tobias Pahl , Poul-Erik Hansen , Liwei Fu , Alexander Birk , Mirza Karamehmedović , Peter Lehmann , Stephan Reichelt , Richard Leach , Samanta Piano","doi":"10.1016/j.jcp.2024.113519","DOIUrl":"10.1016/j.jcp.2024.113519","url":null,"abstract":"<div><div>Rigorous scattering models are based on Maxwell's equations and can provide high-accuracy solutions to model electromagnetic wave scattering from objects. Being able to calculate the scattered field from any surface geometry and considering the effect of the polarisation of the incident light, make rigorous models the most promising tools for complex light-matter interaction problems. The total intensity of the electric near-field scattering from a silicon cylinder illuminated by the transverse electric and transverse magnetic polarisation of the incident light is obtained using various rigorous models including, the local field Fourier modal method, boundary element method and finite element method. The intensity of the total electric near-field obtained by these rigorous models is compared using the Mie solution as a reference for both polarisation modes of the incident light. Additionally, the intensity of the total electric near-field scattered from a silicon sinusoid profile using the same rigorous models is analysed. The results are discussed in detail, and for the cylinder, the deviations in the intensity of the total electric field from the exact Mie solution are investigated.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"521 ","pages":"Article 113519"},"PeriodicalIF":3.8,"publicationDate":"2024-10-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142587080","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

Corrigendum to “Trace preserving quantum dynamics using a novel reparametrization-neutral summation-by-parts difference operator” [J.Comput.Phys. 425 (2021) 109917] 使用新型重参数化中性逐部求和差分算子的痕量保全量子动力学"[J.Comput.Phys. 425 (2021) 109917] 更正

IF 3.8 2区物理与天体物理

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

Oskar Ålund , Yukinao Akamatsu , Fredrik Laurén , Takahiro Miura , Jan Nordström , Alexander Rothkopf

引用次数: 0

Second order conservative Lagrangian DG schemes for compressible flow and their application in preserving spherical symmetry in two-dimensional cylindrical geometry 可压缩流动的二阶保守拉格朗日 DG 方案及其在二维圆柱几何中保持球面对称性的应用

IF 3.8 2区物理与天体物理

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

Wenjing Feng , Juan Cheng , Chi-Wang Shu

{"title":"Second order conservative Lagrangian DG schemes for compressible flow and their application in preserving spherical symmetry in two-dimensional cylindrical geometry","authors":"Wenjing Feng , Juan Cheng , Chi-Wang Shu","doi":"10.1016/j.jcp.2024.113530","DOIUrl":"10.1016/j.jcp.2024.113530","url":null,"abstract":"<div><div>In this paper, we construct a class of second-order cell-centered Lagrangian discontinuous Galerkin (DG) schemes for the two-dimensional compressible Euler equations on quadrilateral meshes. This Lagrangian DG scheme is based on the physical coordinates rather than the fixed reference coordinates, hence it does not require studying the evolution of the Jacobian matrix for the flow mapping between the different coordinates. The conserved variables are solved directly, and the scheme can preserve the conservation property for mass, momentum and total energy. The strong stability preserving (SSP) Runge-Kutta (RK) method is used for the time discretization. Furthermore, there are two main contributions. Firstly, differently from the previous work, we design a new Lagrangian DG scheme which is truly second-order accurate for all the variables such as density, momentum, total energy, pressure and velocity, while the similar DG schemes in the literature may lose second-order accuracy for certain variables, as shown in numerical experiments. Secondly, as an extension and application, we develop a particular Lagrangian DG scheme in the cylindrical geometry, which is designed to be able to preserve one-dimensional spherical symmetry for all the linear polynomials in two-dimensional cylindrical coordinates when computed on an equal-angle-zoned initial grid. The distinguished feature is that it can maintain both the spherical symmetry and conservation properties, which is very important for many applications such as implosion problems. A series of numerical experiments in the two-dimensional Cartesian and cylindrical coordinates are given to demonstrate the good performance of the Lagrangian DG schemes in terms of accuracy, symmetry and non-oscillation.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"521 ","pages":"Article 113530"},"PeriodicalIF":3.8,"publicationDate":"2024-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142560650","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

Simplified inverse distance weighting-immersed boundary method for simulation of fluid-structure interaction 用于模拟流固耦合的简化反距离加权浸润边界法

IF 3.8 2区物理与天体物理

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

Buchen Wu , Yinjie Du , Chang Shu

{"title":"Simplified inverse distance weighting-immersed boundary method for simulation of fluid-structure interaction","authors":"Buchen Wu , Yinjie Du , Chang Shu","doi":"10.1016/j.jcp.2024.113524","DOIUrl":"10.1016/j.jcp.2024.113524","url":null,"abstract":"<div><div>When simulating fluid-structure interaction (FSI) problems involving moving objects, the implicit inverse distance weighting-immersed boundary method (IDW-IBM) developed by Du et al. <span><span>[1]</span></span> has to construct a large square correlation matrix and solve its inversion at each time step. In this work, a simplified inverse distance weighting-immersed boundary method (SIDW-IBM) is proposed to eliminate the intrinsic limitations in the original implicit IDW-IBM. Through error analysis using Taylor series expansion, a second order approximation can be derived, which allows us to approximate the large square correlation matrix into a diagonal matrix; thereby, we proposed the SIDW-IBM based on this second order approximation to explicitly evaluate the velocity corrections, where the needs to assemble the large correlation matrix and inverse it are circumvented. Owing to the fact that the inverse distance weighting interpolation removes the limitations in the Dirac delta function, the proposed SIDW-IBM has been successfully implemented on the non-uniform meshes to further improve the computational efficiency. The proposed SIDW-IBM is integrated with the reconstructed lattice Boltzmann flux solver (RLBFS) <span><span>[2]</span></span> to simulate some classic incompressible viscous flows, including flow past an in-line oscillating cylinder, flow past a heaving airfoil, and flow past a three-dimensional flexible plate. The good agreement between the present results and reference data demonstrates the capability and feasibility of the SIDW-IBM for simulating FSI problems with moving boundaries and large deformations.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"521 ","pages":"Article 113524"},"PeriodicalIF":3.8,"publicationDate":"2024-10-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142553768","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

On the use of fast projection methods with unsteady velocity boundary conditions 关于使用非稳态速度边界条件的快速投影方法

IF 3.8 2区物理与天体物理

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

Maher Eid , Mokbel Karam , Tony Saad

{"title":"On the use of fast projection methods with unsteady velocity boundary conditions","authors":"Maher Eid , Mokbel Karam , Tony Saad","doi":"10.1016/j.jcp.2024.113529","DOIUrl":"10.1016/j.jcp.2024.113529","url":null,"abstract":"<div><div>This article aims at clarifying and amending a previously published result on the use of pseudo-pressure approximations for fast incompressible Navier-Stokes solvers with time-dependent boundary conditions. Fast projection methods were proposed to speed up high-order incompressible Navier-Stokes solvers by replacing pressure projections with simple approximations. A generalized approximation theory was proposed in <span><span>[1]</span></span>, in which the time-dependent boundary conditions were ignored. In this comment, we investigate the effect of unsteady boundary conditions on the order of accuracy and demonstrate that the pressure approximations derived in <span><span>[1]</span></span> hold for all RK2 integrators, while RK3 approximations are only third order for certain sets of coefficients. We show that third order is lost for other cases and shed light on why order is broken. This work serves as a reference for the treatment of unsteady boundary conditions for fast-projection methods. Numerical validation is performed on a well established channel flow problem with a variety of unsteady inlet conditions for a wide range of explicit RK2 and RK3 integrators.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"521 ","pages":"Article 113529"},"PeriodicalIF":3.8,"publicationDate":"2024-10-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142554229","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

Polytopic autoencoders with smooth clustering for reduced-order modeling of flows 具有平滑聚类功能的多拓扑自编码器，用于简化流量建模

IF 3.8 2区物理与天体物理

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

Jan Heiland , Yongho Kim

引用次数: 0

Symmetric form of the magnetohydrodynamics equations 磁流体动力学方程的对称形式

IF 3.8 2区物理与天体物理

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

S.K. Godunov

引用次数: 0

Interesting class of quasilinear systems 有趣的类准线性系统

IF 3.8 2区物理与天体物理

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

S.K. Godunov

引用次数: 0

A fully-coupled algorithm with implicit surface tension treatment for interfacial flows with large density ratios 针对大密度比界面流的隐式表面张力处理全耦合算法

IF 3.8 2区物理与天体物理

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

Romain Janodet , Berend van Wachem , Fabian Denner

{"title":"A fully-coupled algorithm with implicit surface tension treatment for interfacial flows with large density ratios","authors":"Romain Janodet , Berend van Wachem , Fabian Denner","doi":"10.1016/j.jcp.2024.113520","DOIUrl":"10.1016/j.jcp.2024.113520","url":null,"abstract":"<div><div>The stability of most surface-tension-driven interfacial flow simulations is governed by the capillary time-step constraint. This concerns particularly small-scale flows and, more generally, highly-resolved liquid-gas simulations with moderate inertia. To date, the majority of interfacial-flow simulations are performed using an explicit surface-tension treatment, which restrains the performance of such simulations. Recently, an implicit treatment of surface tension able to breach the capillary time-step constraint using the volume-of-fluid (VOF) method was proposed, based on a fully-coupled pressure-based finite-volume algorithm. To this end, the interface-advection equation is incorporated implicitly into the linear flow solver, resulting in a tight coupling between all implicit solution variables (colour function, pressure, velocity). However, this algorithm is limited to uniform density and viscosity fields. Here, we present a fully-coupled algorithm for interfacial flows with implicit surface tension applicable to interfacial flows with large density and viscosity ratios. This is achieved by solving the continuity and momentum equations in conservative form, whereby the density is treated implicitly with respect to the colour function, and the advection term of the interface-advection equation is discretised using the THINC/QQ algebraic VOF scheme, yielding a consistent discretisation of the advective terms. This new algorithm is tested by considering representative surface-tension-dominated interfacial flows, including the Laplace equilibrium of a stationary droplet and the three-dimensional Rayleigh-Plateau instability of a liquid filament. The presented results demonstrate that interfacial flows with large density and viscosity ratios can be simulated and energy conservation is ensured, even with a time step larger than the capillary time-step constraint, provided that other time-step restrictions are satisfied.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"520 ","pages":"Article 113520"},"PeriodicalIF":3.8,"publicationDate":"2024-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142536280","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

Memoirs of finite difference schemes 有限差分方案回忆录

IF 3.8 2区物理与天体物理

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

S.K. Godunov

引用次数: 0