Computer Physics Communications最新文献_第8页

Self-consistent field theory for loop-containing polymers: A general algorithm for path-determination 含环聚合物的自洽场理论：路径确定的通用算法

IF 7.2 2区物理与天体物理

Computer Physics Communications Pub Date : 2024-08-02 DOI: 10.1016/j.cpc.2024.109330

引用次数: 0

A well-balanced finite volume solver for the 2D shallow water magnetohydrodynamic equations with topography 具有地形的二维浅水磁流体动力学方程的均衡有限体积求解器

IF 7.2 2区物理与天体物理

Computer Physics Communications Pub Date : 2024-08-02 DOI: 10.1016/j.cpc.2024.109328

{"title":"A well-balanced finite volume solver for the 2D shallow water magnetohydrodynamic equations with topography","authors":"","doi":"10.1016/j.cpc.2024.109328","DOIUrl":"10.1016/j.cpc.2024.109328","url":null,"abstract":"<div><p>In this paper, a second-order finite volume Non-Homogeneous Riemann Solver is used to obtain an approximate solution for the two-dimensional shallow water magnetohydrodynamic (SWMHD) equations considering non-flat bottom topography. We investigate the stability of a perturbed steady state, as well as the stability of energy in these equations after a perturbation of a steady state using a dispersive analysis. To address the elliptic constraint <span><math><mi>∇</mi><mo>⋅</mo><mi>h</mi><mi>B</mi><mo>=</mo><mn>0</mn></math></span>, the GLM (Generalized Lagrange Multiplier) method designed specifically for finite volume schemes, is used. The proposed solver is implemented on unstructured meshes and verifies the exact conservation property. Several numerical results are presented to validate the high accuracy of our schemes, the well-balanced, and the ability to resolve smooth and discontinuous solutions. The developed finite volume Non-Homogeneous Riemann Solver and the GLM method offer a reliable approach for solving the SWMHD equations, preserving numerical and physical equilibrium, and ensuring stability in the presence of perturbations.</p></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":null,"pages":null},"PeriodicalIF":7.2,"publicationDate":"2024-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141941231","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

Melting phenomena of self-organized magnetic structures investigated by variational autoencoder 利用变异自动编码器研究自组织磁性结构的熔化现象

IF 7.2 2区物理与天体物理

Computer Physics Communications Pub Date : 2024-07-30 DOI: 10.1016/j.cpc.2024.109329

引用次数: 0

libEMMI_MGFD: A program of marine controlled-source electromagnetic modelling and inversion using frequency-domain multigrid solver libEMMI_MGFD: 使用频域多网格求解器进行海洋可控源电磁建模和反演的程序

IF 7.2 2区物理与天体物理

Computer Physics Communications Pub Date : 2024-07-29 DOI: 10.1016/j.cpc.2024.109327

引用次数: 0

Efficient calculation of magnetocrystalline anisotropy energy using symmetry-adapted Wannier functions 利用对称适配的万尼尔函数高效计算磁晶各向异性能

IF 7.2 2区物理与天体物理

Computer Physics Communications Pub Date : 2024-07-29 DOI: 10.1016/j.cpc.2024.109325

引用次数: 0

KelbgLIP: Program implementation of the high-temperature Kelbg density matrix for path integral and molecular dynamics simulations with long-range Coulomb interaction KelbgLIP：高温 Kelbg 密度矩阵的程序实现，用于具有长程库仑相互作用的路径积分和分子动力学模拟

IF 7.2 2区物理与天体物理

Computer Physics Communications Pub Date : 2024-07-29 DOI: 10.1016/j.cpc.2024.109326

{"title":"KelbgLIP: Program implementation of the high-temperature Kelbg density matrix for path integral and molecular dynamics simulations with long-range Coulomb interaction","authors":"","doi":"10.1016/j.cpc.2024.109326","DOIUrl":"10.1016/j.cpc.2024.109326","url":null,"abstract":"<div><p>In this paper, we present the <span>KelbgLIP</span> code to implement the previously obtained analytical density matrix that includes Coulomb long-range interactions. The method is based on the work of G. Kelbg, who derived a high temperature density matrix for the Coulomb potential. To include all long-range interactions in the density matrix, we use the Ewald technique, specifically the angular-averaged Ewald potential (AAEP). The solution of the Blöch equation within the AAEP has a direct analytic form that can be easily implemented in classical and quantum Monte Carlo or molecular dynamics codes, including exchange effects. The potential part of the density matrix remains finite at small distances, preventing the collapse of a two-component system. Using <span>KelbgLIP</span>, one can calculate the diagonal Kelbg-AAE pseudopotential and the pair density matrix. In the case of a hydrogen plasma, the code is able to calculate action, kinetic and potential energy in the path integral representation. We validated our approach by simulating a nondegenerate weakly coupled hydrogen plasma and obtained the thermodynamic limit in agreement with the Debye-Hückel approximation. In addition, we observe the agreement with classical simulations using the unbounded from below AAEP, which is possible in the weakly-coupled regime.</p></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":null,"pages":null},"PeriodicalIF":7.2,"publicationDate":"2024-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141941234","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

Adaptive unified gas-kinetic scheme for diatomic gases with rotational and vibrational nonequilibrium 具有旋转和振动非平衡的双原子气体的自适应统一气体动力学方案

IF 7.2 2区物理与天体物理

Computer Physics Communications Pub Date : 2024-07-26 DOI: 10.1016/j.cpc.2024.109324

{"title":"Adaptive unified gas-kinetic scheme for diatomic gases with rotational and vibrational nonequilibrium","authors":"","doi":"10.1016/j.cpc.2024.109324","DOIUrl":"10.1016/j.cpc.2024.109324","url":null,"abstract":"<div><p>Multiscale nonequilibrium physics at large variations of local Knudsen number are encountered in applications of aerospace engineering and micro-electro-mechanical systems, such as high-speed flying vehicles and low pressure of the encapsulation. An accurate description of flow physics in all flow regimes within a single computation requires a genuinely multiscale method. The adaptive unified gas-kinetic scheme (AUGKS) is developed for such multiscale flow simulation. The AUGKS applies discretized velocity space to accurately capture the non-equilibrium physics in the multiscale UGKS, and adaptively employs continuous distribution functions following Chapman–Enskog expansion to efficiently recover near-equilibrium flow region in GKS. The UGKS and GKS are dynamically connected at the cell interface through the fluxes from the discretized and continuous gas distribution functions, which avoids any buffer zone between them. In this study, the AUGKS with rotation and vibration non-equilibrium is developed based on a multiple temperature relaxation model. The real gas effect in different flow regimes has been properly captured. To capture aerodynamic heating accurately, the heat flux modifications from the rotation and vibration modes are also included in the current scheme. Unstructured discrete particle velocity space is adopted to further improve the computational performance of the AUGKS. Numerical tests, including Sod tube, normal shock structure, high-speed flow around the two-dimensional cylinder and three-dimensional sphere and space vehicles, and an unsteady nozzle plume flow from the continuum flow to the background vacuum, have been conducted to validate the current scheme. In comparison with the original UGKS, the current scheme speeds up the computation, reduces the memory requirement, and maintains the equivalent accuracy for multiscale flow simulation, which provides an effective tool for nonequilibrium flow simulations, especially for the flows at low and medium speed.</p></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":null,"pages":null},"PeriodicalIF":7.2,"publicationDate":"2024-07-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141941235","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 second-order particle Fokker-Planck model for rarefied gas flows 稀薄气体流的二阶粒子福克-普朗克模型

IF 7.2 2区物理与天体物理

Computer Physics Communications Pub Date : 2024-07-25 DOI: 10.1016/j.cpc.2024.109323

{"title":"A second-order particle Fokker-Planck model for rarefied gas flows","authors":"","doi":"10.1016/j.cpc.2024.109323","DOIUrl":"10.1016/j.cpc.2024.109323","url":null,"abstract":"<div><p>The direct simulation Monte Carlo (DSMC) method has become a powerful tool for studying rarefied gas flows. However, for the DSMC method to be effective, the cell size must be smaller than the mean free path, and the time step smaller than the mean collision time. These constraints make it difficult to use the DSMC method in multiscale rarefied gas flows. Over the past decade, the particle Fokker-Planck (FP) method has been studied to address computational cost issues in the near-continuum regime. To capture the main features of the Boltzmann equation, various FP models have been proposed, such as the quadratic entropic FP (Quad-EFP) and the ellipsoidal statistical FP (ESFP). Nevertheless, few studies have clearly demonstrated that the FP method offers a computational advantage over the DSMC method without sacrificing accuracy. This is because conventional particle FP methods have employed first-order accuracy schemes. The present study proposes a unified stochastic particle ESFP (USP-ESFP) model. This model improves the accuracy of shear stress and heat flux predictions. Additionally, a spatial interpolation scheme is introduced to the particle FP method. The numerical test cases include relaxation problem, Couette flows, Poiseuille flows, velocity perturbation, and hypersonic flows around a cylinder. The results show that the USP-ESFP model agrees well with both analytical and DSMC results. Furthermore, the USP-ESFP model is found to be less sensitive to cell size and time step than the DSMC method, resulting in a factor of four speed-up for the considered hypersonic flow around a cylinder.</p></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":null,"pages":null},"PeriodicalIF":7.2,"publicationDate":"2024-07-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0010465524002467/pdfft?md5=ddd38f537108ba1509755299dddad47b&pid=1-s2.0-S0010465524002467-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141844852","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

The final version of the 5D histogram package NORA 5D 直方图软件包 NORA 的最终版本

IF 7.2 2区物理与天体物理

Computer Physics Communications Pub Date : 2024-07-25 DOI: 10.1016/j.cpc.2024.109320

{"title":"The final version of the 5D histogram package NORA","authors":"","doi":"10.1016/j.cpc.2024.109320","DOIUrl":"10.1016/j.cpc.2024.109320","url":null,"abstract":"<div><p>The presented software package is an advanced analogue of the famous HBOOK and, in part, ROOT packages.</p><p>The main features of the package are as follows.</p><p>Standard operations (accumulation, simulation, transformation) are extended up to 5D objects.</p><p>A new type of transformations of objects has been introduced, the x- transformation, which includes convolution of distributions.</p><p>Objects are accessed mainly by alphabetical name.</p><p>The formation of the object space is carried out using a data file, where the user can choose the form of setting attributes.</p><p>Automatic adjustment of object attributes is possible (only the number of channels is set).</p><p>The number of channels for each of the axes of the object is unlimited.</p><p>Operations between two objects (addition, subtraction, etc.) are possible with mismatched attributes. Only the ranks of the objects should match.</p><p>Data output is carried out in two forms, for graphics and for fit.</p><p>Group operations are provided for visualization, outputting files for graphics.</p><p>In the last version the program for x-converting 2D to 2D has been added.</p><p>A group of small programs has been added to simplify the user interface.</p><p>Improved representation of 2D objects, while removing the limit on the number of channels for all axes.</p><p>A number of bugs are fixed.</p><p>All programmes are written in FORTRAN-90.</p><p>The investigation has been performed at the Veksler and Baldin Laboratory of High Energy Physics, JINR.</p></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":null,"pages":null},"PeriodicalIF":7.2,"publicationDate":"2024-07-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141844064","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 novel model for direct numerical simulation of suspension dynamics with arbitrarily shaped convex particles 用于直接数值模拟任意形状凸颗粒悬浮动力学的新型模型

IF 7.2 2区物理与天体物理

Computer Physics Communications Pub Date : 2024-07-18 DOI: 10.1016/j.cpc.2024.109321

{"title":"A novel model for direct numerical simulation of suspension dynamics with arbitrarily shaped convex particles","authors":"","doi":"10.1016/j.cpc.2024.109321","DOIUrl":"10.1016/j.cpc.2024.109321","url":null,"abstract":"<div><p>This study presents an innovative direct numerical simulation approach for complex particle systems with irregular shapes and large numbers. Using partially saturated methods, it accurately models arbitrary shapes, albeit at considerable computational cost when integrating a compatible contact model. The introduction of a novel parallelization strategy significantly improves the performance of the contact model, enabling efficient four-way coupled simulations. Through hindered settling studies, the criticality of the explicit contact model for maintaining simulation accuracy is highlighted, especially at high particle volume fractions and low Archimedes numbers. The feasibility of simulating thousands of arbitrarily shaped convex particles is demonstrated with up to 1934 surface-resolved particles. The study also confirms the grid independence and linear convergence of the method. It shows for the first time that cube swarms settle 13 to 26% slower than swarms of volume-equivalent spheres across different Archimedes numbers (500 to 2000) and particle volume fractions (10 to 30%). These findings emphasize the shape dependence of particle systems and suggest avenues for exploring their nuanced dynamics.</p></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":null,"pages":null},"PeriodicalIF":7.2,"publicationDate":"2024-07-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0010465524002443/pdfft?md5=aefa678974806198e79fd9446b434c5b&pid=1-s2.0-S0010465524002443-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141729600","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