样式: 排序: IF: - GO 导出 标记为已读
-
Peak finding algorithm for cluster counting with domain adaptation Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-04-16 Guang Zhao, Linghui Wu, Francesco Grancagnolo, Nicola De Filippis, Mingyi Dong, Shengsen Sun
Cluster counting in drift chamber is the most promising breakthrough in particle identification (PID) technique in particle physics experiment. Reconstruction algorithm is one of the key challenges in cluster counting. In this paper, a semi-supervised domain adaptation (DA) algorithm is developed and applied on the peak finding problem in cluster counting. The algorithm uses optimal transport (OT)
-
Fast exact algorithm for neutrino oscillation in constant matter density Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-04-15 James Page
A recently published method for solving the neutrino evolution equation with constant matter density is further refined and used to lay out an exact algorithm for computing oscillation probabilities, which is moderately faster than previous methods when looping through neutrinos of different energies. In particular, the three examples of survival, survival and appearance probabilities are written in
-
AutoEFT: Automated operator construction for effective field theories Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-04-12 Robert V. Harlander, Magnus C. Schaaf
-
Nektar++: Development of the compressible flow solver for jet aeroacoustics Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-04-10 Daniel Lindblad, João Isler, Margarida Moragues Ginard, Spencer J. Sherwin, Chris D. Cantwell
A recently developed computational framework for jet noise predictions is presented. The framework consists of two main components, focusing on source prediction and noise propagation. To compute the noise sources, the turbulent jet is simulated using the compressible flow solver implemented in the open-source spectral/hp element framework Nektar++, which solves the unfiltered Navier-Stokes equations
-
Constraint free physics-informed machine learning for micromagnetic energy minimization Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-04-10 Sebastian Schaffer, Lukas Exl
We introduce a novel method for micromagnetic energy minimization which uses physics-informed neural networks to find a magnetic configuration which minimizes the Gibbs Free energy functional without the need of any constraint optimization framework. The Cayley transform is applied to a neural network to assure that the model output lives on the Lie group of rotation matrices . For the stray field
-
SHORYUKEN: An open-source software package for calculating nonlocal exchange interactions in nanowires Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-04-09 Yuan Chen, Simon N. Sandhofer, Bryan M. Wong
We present the open-source software package, SHORYUKEN (Streamlined High-level Operations in Real-space to Yield, Understand, and Keep Exchange in Nanowires), for calculating nonlocal exchange interactions in nanowires with arbitrary geometries, sizes, doping densities, and compositions. In addition to enabling new calculations of nonlocal exchange, the SHORYUKEN software package is a significant enhancement
-
Chrono DEM-Engine: A Discrete Element Method dual-GPU simulator with customizable contact forces and element shape Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-04-08 Ruochun Zhang, Bonaventura Tagliafierro, Colin Vanden Heuvel, Shlok Sabarwal, Luning Bakke, Yulong Yue, Xin Wei, Radu Serban, Dan Negruţ
This paper introduces DEM-Engine, a new submodule of Project Chrono, that is designed to carry out Discrete Element Method (DEM) simulations. Based on spherical primitive shapes, DEM-Engine can simulate polydisperse granular materials and handle complex shapes generated as assemblies of primitives, referred to as . DEM-Engine has a multi-tier parallelized structure that is optimized to operate simultaneously
-
ReMKiT1D - A framework for building reactive multi-fluid models of the tokamak scrape-off layer with coupled electron kinetics in 1D Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-04-05 Stefan Mijin, Dominic Power, Ryan Holden, William Hornsby, David Moulton, Fulvio Militello
-
A parallel and performance portable implementation of a full-field crystal plasticity model Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-04-04 Caleb O. Yenusah, Nathaniel R. Morgan, Ricardo A. Lebensohn, Miroslav Zecevic, Marko Knezevic
We have developed a parallel implementation of an Elasto-Viscoplastic Fast Fourier Transform-based (EVPFFT) micromechanical solver to enable computationally efficient crystal plasticity modeling for polycrystalline materials. Our primary focus lies in achieving performance portability, allowing a single EVPFFT implementation to run optimally on various homogeneous architectures, including multi-core
-
Comparison of two different integration methods for the (1+1)-dimensional Schrödinger-Poisson equation Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-29 Nico Schwersenz, Victor Loaiza, Tim Zimmermann, Javier Madroñero, Sandro Wimberger
We compare two different numerical methods to integrate in time spatially delocalized initial densities using the Schrödinger-Poisson equation system as the evolution law. The basic equation is a nonlinear Schrödinger equation with an auto-gravitating potential created by the wave function density itself. The latter is determined as a solution of Poisson's equation modelling, e.g., non-relativistic
-
Time-domain event detection using single-instruction, multiple-thread gpGPU architectures in single-molecule biophysical data Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-28 Boyan Penkov, David Niedzwiecki, Nicolae Lari, Marija Drndić, Kenneth Shepard
Discrete amplitude levels in ordered, time-domain data often represent different underlying latent states of the system that is being interrogated. Analysis and feature extraction from these data sets generally requires considering the order of each individual point; this approach cannot take advantage of contemporary general-purpose graphics processing units (gpGPU) and single-instruction multiple-data
-
Heterogeneous CPU-GPU parallelization for modeling supersonic reacting flows with detailed chemical kinetics Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-28 Sihang Rao, Bing Chen, Xu Xu
Accurate simulations of complex combustion phenomena associated with supersonic reacting flows require the use of detailed chemical kinetic mechanisms. However, detailed chemical mechanisms may consist of a large number of species and reactions resulting in extremely high computation cost. In order to accelerate the simulations of supersonic reacting flows with detailed chemical mechanisms, a heterogeneous
-
SMATool: Strength of materials analysis toolkit Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-28 C.E. Ekuma
The study of the strength of materials is a cornerstone in material science and engineering, playing a critical role in shaping the progress and application of materials in diverse industrial sectors. The strength of a material is meticulously examined to understand the behavior of the material under different stress conditions and environments, thereby guiding material selection and structural design
-
XLB: A differentiable massively parallel lattice Boltzmann library in Python Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-27 Mohammadmehdi Ataei, Hesam Salehipour
The lattice Boltzmann method (LBM) has emerged as a prominent technique for solving fluid dynamics problems due to its algorithmic potential for computational scalability. We introduce XLB library, a Python-based differentiable LBM library based on the JAX platform. The architecture of XLB is predicated upon ensuring accessibility, extensibility, and computational performance, enabling scaling effectively
-
Performance modeling and running strategy of parallel cdugksFOAM program Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-25 Yunlan Wang, Yufeng Liu, Rui Zhang, Tianhai Zhao, Sha Liu, Congshan Zhuo, Chengwen Zhong
The cdugksFOAM program realizes the physical space grid and the velocity space grid in parallel at the same time. Its distinguishing feature lies in its potential for large-scale parallelism. However, the running time of the cdugksFOAM program is significantly dependent on the number of physical and velocity space partitions. In order to find the optimal partitioning strategies for a specific CFD problem
-
A custom detector construction pattern for Geant4 applications Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-25 Mustafa Kandemir
Geant4 Detector Construction Pattern (G4DCP) is a template developed to flexibly construct complex detectors in Geant4 applications. The elements of G4DCP, including , form an elegant template for detector setups. We construct a sample detector geometry utilizing this template and make the developed code available to the public. G4DCP
-
Stokes-Einstein relation for binary mixtures Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-24 Yang Liu, Dietmar Block
The applicability of the Stokes-Einstein (SE) relation in two-dimensional finite binary mixtures is tested using the Langevin simulation method. Calculations of viscosity and self-diffusion coefficients are compared between monodisperse and binary systems. It is shown that adapted definitions of coupling strength Γ and screening parameter allow to describe the transport properties for both monodisperse
-
Local discontinuous Galerkin for the functional renormalisation group Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-24 Friederike Ihssen, Jan M. Pawlowski, Franz R. Sattler, Nicolas Wink
We apply a Local Discontinuous Galerkin discretisation to flow equations of the O(N)-model in the Local Potential Approximation. The improved stability is directly observed by solving the flow equation for various and space-time dimensions . A particular focus of this work is the numerical discretisation and its implementation. It is realised as a module within the high performance PDE framework DUNE
-
Low-mode deflation for twisted-mass and RHMC reweighting in lattice QCD Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-21 Simon Kuberski
We propose improved estimators to compute the reweighting factors which are needed for lattice QCD calculations that rely on twisted-mass reweighting for the light quark contribution and the Rational Hybrid Monte Carlo (RHMC) algorithm for non-degenerate quark masses. This is the case for a number of modern large-scale simulations based on improved Wilson fermions. We find a significant reduction of
-
Comparison of different discontinuous Galerkin methods based on various reformulations for gKdV equation: Soliton dynamics and blowup Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-21 Xue Hong, Qianrui Wei, Xiaofei Zhao
In this work, we explore the use of several discontinuous Galerkin (DG) methods for simulating soliton dynamics in the generalized Korteweg-de Vries (gKdV) equation. The presence of high-order nonlinearity in this model can lead to the finite-time blowup phenomenon, which challenges every aspect of a numerical scheme. The considered DG schemes encompass popular methods derived from the energy and/or
-
Arbitrary controlled re-orientation of a spinning body by evolving its tensor of inertia Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-21 Igor A. Ostanin, Matthias Sperl
Bodies with the nonspherical tensor of inertia (TOI) exhibit a variety of rotational motion patterns, including chaotic motion, stable periodic (quasi-periodic) rotation, unstable rotation around the direction close to the body's second principal axis, featuring a well-known tennis-racket (also known as Garriott-Dzhanibekov ) effect – series of seemingly spontaneous 180 degrees flips. These patterns
-
SCTBEM: A scaled coordinate transformation boundary element method with 99-line MATLAB code for solving Poisson's equation Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-20 Bo Yu, Ruijiang Jing
This study introduces the scaled coordinate transformation boundary element method (SCTBEM), a novel boundary-type method for solving three-dimensional potential problems. To apply the BEM to complex problems, it is common practice to use the fundamental solution corresponding to the partial governing equation operator to establish the integral equation. However, this approach introduces the domain
-
mqdtfit: A collection of Python functions for empirical multichannel quantum defect calculations Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-20 R.M. Potvliege
The Python functions distributed with this article can be used for calculating the parameters of multichannel quantum defect theory models describing excited bound states of complex atoms. These parameters are obtained by fitting a model to experimental data provided by the user. The two main formulations of the theory are supported, namely the one in which the parameters of the model are a set of
-
FiniteFieldSolve: Exactly solving large linear systems in high-energy theory Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-20 James Mangan
-
Simulating CDT quantum gravity Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-19 Joren Brunekreef, Renate Loll, Andrzej Görlich
We provide a hands-on introduction to Monte Carlo simulations in nonperturbative lattice quantum gravity, formulated in terms of Causal Dynamical Triangulations (CDT). We describe explicitly the implementation of Monte Carlo moves and the associated detailed-balance equations in two and three spacetime dimensions. We discuss how to optimize data storage and retrieval, which are nontrivial due to the
-
Supervised training of neural-network quantum states for the next-nearest neighbor Ising model Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-19 Zheyu Wu, Remmy Zen, Heitor P. Casagrande, Dario Poletti, Stéphane Bressan
Different neural network architectures can be unsupervisedly or supervisedly trained to represent quantum states. We explore and compare different strategies for the supervised training of feed forward neural network quantum states. We empirically and comparatively evaluate the performance of feed forward neural network quantum states in different phases of matter for variants of the architecture,
-
An improved methodology for modeling short pulse buried layer x-ray emission spectra Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-15 D.T. Cliche, M.E. Martin, R.A. London, H.A. Scott, M.V. Patel
Radiation-hydrodynamic and spectroscopic modeling are important aspects of high energy density experimental design. In this paper, we improve the performance and capabilities over those obtainable with a previous methodology used for simulating x-ray emission spectra from buried layer targets heated by short pulse lasers. The improvement incorporates post-processing HYDRA radiation-hydrodynamic output
-
On polytope intersection by half-spaces and hyperplanes for unsplit geometric volume of fluid methods on arbitrary grids Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-13 Joaquín López, Julio Hernández
Several tools with improved accuracy and computational efficiency, along with added capabilities, are presented for performing analytic and geometric operations involving intersections between half-spaces or hyperplanes and 2- or 3-polytopes (polygons or polyhedra, either convex or non-convex) that typically arise in advanced volume of fluid (VOF) methods on arbitrary grids. In particular, tools for
-
A physics-based tessellation algorithm for particle assemblies on arbitrary surfaces Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-12 Shah Wasif Sazzad, Sanjay Dharmavaram, Luigi E. Perotti
Interacting particle assemblies embedded on a surface are often used to model biophysical systems and to study new colloidal materials. The configurations resulting from the particles interacting with each other and with their substrate affect the system's physical properties, which depend on local symmetries and defects. It is therefore important to identify the nature and location of defects in the
-
SIMULATeQCD: A simple multi-GPU lattice code for QCD calculations Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-08 Lukas Mazur, Dennis Bollweg, David A. Clarke, Luis Altenkort, Olaf Kaczmarek, Rasmus Larsen, Hai-Tao Shu, Jishnu Goswami, Philipp Scior, Hauke Sandmeyer, Marius Neumann, Henrik Dick, Sajid Ali, Jangho Kim, Christian Schmidt, Peter Petreczky, Swagato Mukherjee, (HotQCD collaboration)
The rise of exascale supercomputers has fueled competition among GPU vendors, driving lattice QCD developers to write code that supports multiple APIs. Moreover, new developments in algorithms and physics research require frequent updates to existing software. These challenges have to be balanced against constantly changing personnel. At the same time, there is a wide range of applications for HISQ
-
Multipolynomial Monte Carlo for trace estimation in lattice QCD Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-08 Paul Lashomb, Ronald B. Morgan, Travis Whyte, Walter Wilcox
Estimating the trace of the inverse of a large matrix is an important problem in lattice quantum chromodynamics. A multilevel Monte Carlo method is proposed for this problem that uses different degree polynomials for the levels. The polynomials are developed from the GMRES algorithm for solving linear equations. To reduce orthogonalization expense, the highest degree polynomial is a composite or double
-
Olsson.wl & ROC2.wl: Mathematica packages for transformations of multivariable hypergeometric functions & regions of convergence for their series representations in the two variables case Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-08 B. Ananthanarayan, Souvik Bera, S. Friot, Tanay Pathak
-
ElasTool v3.0: Efficient computational and visualization toolkit for elastic and mechanical properties of materials Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-07 C.E. Ekuma, Z.-L. Liu
Efficient computation and visualization of elastic and mechanical properties are crucial in the selection of materials and the design of new materials. The toolkit marks a significant advancement in the computational analysis and visualization of elastic and mechanical properties of materials, essential in material selection and design. This enhanced version extends beyond standard calculations like
-
Band structure calculations of three-dimensional solid-fluid coupling phononic crystals using dual reciprocity boundary element method and wavelet compression method Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-07 Qi Wei, Jiawei Xiang, Weiping Zhu, Hongjiu Hu
The algorithm for calculating the band structures of two-dimensional phononic crystals (PCs) using the boundary element method (BEM) has been proposed for many years. However, it has not yet been extended to three-dimensional (3D) PCs because the fundamental solutions of 3D dynamics are complex and are related to angular frequency. In this study, the BEM is applied to calculate the band structures
-
PolyWeight: A free and open-source program for determination of molecular weight distribution of linear polymers Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-07 Atilio Minotto Neto, Otávio Bianchi, Leonardo Bresciani Canto, Janete Eunice Zorzi, Cláudio Antônio Perottoni
This paper introduces PolyWeight, a Python software featuring a user-friendly graphical user interface (GUI), which offers two distinct approaches for MWD determination: an analytical relation-based method and a parametric model-based method. By utilizing dynamic moduli, users can calculate MWD as well as molecular weight averages such as , , and . The functionality of PolyWeight is validated using
-
A class of arbitrarily high-order energy-preserving method for nonlinear Klein–Gordon–Schrödinger equations Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-04 Xuelong Gu, Yuezheng Gong, Wenjun Cai, Yushun Wang
In this paper, we develop a class of arbitrarily high-order energy-preserving time integrators for the nonlinear Klein–Gordon–Schrödinger equations. We employ Fourier pseudo-spectral method for spatial discretization, resulting in a semi-discrete system. Subsequently, we employ the Petrov-Galerkin method in time to obtain a fully-discrete system. We rigorously demonstrate that the proposed scheme preserves
-
An efficient solution of the multi-term multi-harmonic electron Boltzmann equation for use in global models Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-01 Joel E. Lynch, Travis R. Sippel, Shankar Subramaniam
Solving the electron Boltzmann equation is an essential but costly step in simulating low-temperature plasma kinetics. This work addresses the problem by introducing a solution of the general multi-term multi-harmonic Boltzmann equation (MTMH-BE) optimized for electrons in time-dependent non-equilibrium gases and electric fields. This is accomplished by configuring the numerical Jacobian as the product
-
A high-efficiency adaptive TENO scheme with optimal accuracy order for compressible flow simulation Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-03-01 Shujiang Tang
In this paper, we propose a new adaptive cut-off function and develop a fifth-order targeted ENO scheme that achieves optimal accuracy at any order of critical points, which performs excellently in conventional compressible gas dynamics. The cut-off value is crucial for controlling the dissipation in TENO schemes. Adjusting can improve shock capture and increase dissipation for small-scale features
-
Neutrinos from muon-rich ultra high energy electromagnetic cascades: The MUNHECA code Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-29 AmirFarzan Esmaeili, Arman Esmaili, Pasquale Dario Serpico
-
Verification of the Fourier-enhanced 3D finite element Poisson solver of the gyrokinetic full-f code PICLS Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-29 A. Stier, A. Bottino, M. Boesl, M. Campos Pinto, T. Hayward-Schneider, D. Coster, A. Bergmann, M. Murugappan, S. Brunner, L. Villard, F. Jenko
We introduce and derive the Fourier-enhanced 3D electrostatic field solver of the gyrokinetic full-f PIC code PICLS. The solver makes use of a Fourier representation in one periodic direction of the domain to make the solving of the system easily parallelizable and thus save run time. The presented solver is then verified using two different approaches of manufactured solutions. The test setup used
-
Peridynamic modeling for multiscale heat transport of phonon Boltzmann transport equation Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-29 Weier Liu, Yangde Feng, Ruilin Li, Chenhan Bai, Beifang Niu
Phonons are the main carriers in semiconductor materials, and the Boltzmann transport equation (BTE) can describe the phonon heat transport well. Numerically solving the phonon BTE is a challenging task due to its high dimensionality and nonlinearity. In this work, we develop a Peridynamic model for steady-state phonon heat transport of phonon Boltzmann transport equation based on the Peridynamic differential
-
The polarimeter vector for τ → 3πντ decays Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-28 Vladimir Cherepanov, Christian Veelken
The polarimeter vector of the represents an optimal observable for the measurement of the spin. In this paper we present an algorithm for the computation of the polarimeter vector for the decay channels and . The algorithm is based on a model for the hadronic current in these decay channels, which was fitted to data recorded by the CLEO experiment . PolarimetricVectorTau2a1, version 1.0.1
-
BIMBAMBUM: A potential flow solver for single cavitation bubble dynamics Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-28 Armand Baptiste Sieber, Henri Hugo Sieber, Davide Bernardo Preso, Mohamed Farhat
In the absence of analytical solutions for the dynamics of non-spherical cavitation bubbles, we have implemented a numerical simulation solver based on the boundary integral method (BIM) that models the behavior of a single bubble near an interface between two fluids. The density ratio between the two media can be adjusted to represent different types of boundaries, such as a rigid boundary or a free
-
Comparison of effective and stable Langevin dynamics integrators Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-23 Bogdan Tanygin, Simone Melchionna
Langevin and Brownian simulations play a prominent role in computational research, and state of the art integration algorithms provide trajectories with different stability ranges and accuracy in reproducing statistical averages. The practical usability of integrators is an important aspect to allow choosing large time steps while ensuring numerical stability and overall computational efficiency. In
-
ERCS24: An updated version of the ERCS08 program for calculations of the cross sections for atomic electron removal based on the ECPSSR theory and its variants Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-23 Vladimir Horvat
ERCS24, an updated version of the ERCS08 program, calculates the atomic lectron emoval ross ections. It is written in FORTRAN in order to make it more portable and easier to customize by a large community of physicists, but it also comes with a separate windows graphics user interface control application ERCS24w that makes it easy to quickly prepare the input file, run the program, as well as view
-
Massively parallel implementation of iterative eigensolvers in large-scale plane-wave density functional theory Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-22 Junwei Feng, Lingyun Wan, Jielan Li, Shizhe Jiao, Xinhui Cui, Wei Hu, Jinlong Yang
The Kohn-sham density functional theory (DFT) is a powerful method to describe the electronic structures of molecules and solids in condensed matter physics, computational chemistry and materials science. However, large and accurate DFT calculations within plane waves process a cubic-scaling computational complexity, which is usually limited by expensive computation and communication costs. The rapid
-
Quasi-optimal domain decomposition method for neural network-based computation of the time-dependent Schrödinger equation Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-21 Emmanuel Lorin, Xu Yang
In this paper, we derive and analyze the performance of optimal/quasi-optimal Schwarz Waveform Relaxation (SWR) domain decomposition methods (DDM) for the time-dependent Schrödinger equation when implement with neural network-based Partial Differential Equations (PDE) solvers. Optimal SWR methods, which are based on Dirichlet-to-Neumann operators, are known to have a higher convergence rate than classical
-
A binary filter inspired from the PIC sparse grid technique – Illustration on the XTOR-K code Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-21 T. Nicolas, V. Dubois, Q. Fang, H. Lütjens
It is known that the sparse grid method for Particle-In-Cell (PIC) solvers acts as a filter to reduce the PIC noise. In this paper, a simple rule to discard or keep modes in Fourier space (a binary filter with values either 0 or 1) is derived using the sparse grid combination formula. Its relation to the standard sparse grid filter, which is characterized quantitatively, is explained. The relations
-
ePDFpy: A Python-based interactive GUI tool for electron pair distribution function analysis of amorphous materials Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-19 Minhyo Kim, Pilsung Kim, Riccardo Bassiri, Kiran Prasai, Martin M. Fejer, Kyung-ha Lee
ePDFpy is an interactive analysis program with a graphical user interface (GUI), designed to process the electron Pair Distribution Function (PDF) analysis of diffraction patterns from Transmission Electron Microscope (TEM), to identify the local atomic structure of amorphous materials. The program offers a user-friendly Python-based interface, providing a straightforward and adaptable workflow for
-
Attosecond Chemistry Special Issue Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-16 Jimena D. Gorfinkiel, Fernando Martín
-
Development of a tokamak magnetohydrodynamic code with the discontinuous Galerkin and Weighted Essentially Non-Oscillatory methods Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-15 J. Ma, W. Guo, Y. Xie
In this study, we present the development of a new initial-value magnetohydrodynamic (MHD) code for toroidal geometry using discontinuous Galerkin (DG) and Weighted Essentially Non-Oscillatory (WENO) methods. The code utilizes a triangular mesh based on the flux of the fixed boundary equilibrium in the poloidal plane, which is uniformly divided in the toroidal direction. By solving the conservative
-
micrOMEGAs 6.0: N-component dark matter Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-15 G. Alguero, G. Bélanger, F. Boudjema, S. Chakraborti, A. Goudelis, S. Kraml, A. Mjallal, A. Pukhov
is a numerical code to compute dark matter (DM) observables in generic extensions of the Standard Model (SM) of particle physics. We present a new version of that includes a generalization of the Boltzmann equations governing the DM cosmic abundance evolution which can be solved to compute the relic density of N-component DM. The direct and indirect detection rates in such scenarios take into account
-
ERSN-OpenMC-Py: A python-based open-source software for OpenMC Monte Carlo code Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-15 M. Lahdour, T. El Bardouni, O. El Hajjaji, J. EL Bakkali, J. Al-Zain, S. Oulad-Belayachi, H. Ziani, Abdelghani Idrissi, S. El Maliki El Hlaibi
The graphical user interface is a key element in facilitating the use of complex simulation software. This project describes the development of a graphical user interface called “ERSN-OpenMC-Py” for an existing neutron simulation code, OpenMC. The main goal is to make simulation more accessible to a wider audience by providing a user-friendly and intuitive user interface. The process of developing
-
Trans-Net: A transferable pretrained neural networks based on temporal domain decomposition for solving partial differential equations Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-15 Dinglei Zhang, Ying Li, Shihui Ying
Physics-Informed Neural Networks (PINNs) has provided a novel direction for solving partial differential equations (PDEs) and has achieved significant advancements in the field of scientific computing. PINNs effectively incorporate the physical constraints of equations into the loss function, enabling neural networks to learn and approximate the behavior of physical systems by optimizing the loss function
-
A well-balanced all-Mach scheme for compressible two-phase flow Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-15 Sandro Malusà, Alessandro Alaia
We present an implicit-explicit finite volume scheme for the compressible two-phase model in all-Mach number regimes. In order to solve model equations efficiently and accurately in the low Mach regime, the convective term is split in a stiff part associated to fast acoustic waves, and a non-stiff part corresponding to mean flow advection. A Implicit-Explicit Runge-Kutta (IMEX-RK) method is used to
-
Evaluation of classical correlation functions from 2/3D images on CPU and GPU architectures: Introducing CorrelationFunctions.jl Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-15 Vasily Postnicov, Aleksei Samarin, Marina V. Karsanina, Mathieu Gravey, Aleksey Khlyupin, Kirill M. Gerke
Correlation functions are becoming one of the major tools for quantification of structural information that is usually represented as 2D or 3D images. In this paper we introduce ▪ open-source package developed in Julia and capable of computing all classical correlation functions based on imaging input data. Images include both binary and multi-phase representations. Our code is capable of evaluating
-
Handling shape optimization of superconducting cavities with DNMOGA Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-15 Peilin Wang, Kuangkuang Ye, Xuerui Hao, Jike Wang
Radiofrequency (RF) cavities hold immense importance in various accelerator applications, but their optimization poses significant challenges due to complex situations involved. In this study, a recently proposed multi-objective optimization algorithm is utilized to optimize the 325 MHz double spoke cavity, which is characterized by 38 geometric parameters and is one of the most complex cavities commonly
-
PolyHoop: Soft particle and tissue dynamics with topological transitions Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-15 Roman Vetter, Steve V.M. Runser, Dagmar Iber
We present PolyHoop, a lightweight standalone C++ implementation of a mechanical model to simulate the dynamics of soft particles and cellular tissues in two dimensions. With only few geometrical and physical parameters, PolyHoop is capable of simulating a wide range of particulate soft matter systems: from biological cells and tissues to vesicles, bubbles, foams, emulsions, and other amorphous materials
-
A SPIRED code for the reconstruction of spin distribution Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-15 Simon Buchwald, Gabriele Ciaramella, Julien Salomon, Dominique Sugny
-
Identification of self-interstitial atoms and vacancies in crystalline materials in atomistic simulation Comput. Phys. Commun. (IF 6.3) Pub Date : 2024-02-14 Jun Wang, Tao Li, Ziwen Fu, Baoqin Fu, Chengjun Gou
As the most important intrinsic point defects, self-interstitial atoms and vacancies play crucial roles in governing the microstructural evolution in crystalline materials. In atomistic simulation, the self-interstitial atoms and vacancies are identified by various location methods, such as the Wigner-Seitz cell method. Here we reveal that the mobile defects will give rise to the random drift of the