### Top 9 Arxiv Papers Today in Computational Physics

##### #1. Leveraging Bayesian Analysis To Improve Accuracy of Approximate Models
###### Balasubramanya T. Nadiga, Chiyu Jiang, Daniel Livescu
We focus on improving the accuracy of an approximate model of a multiscale dynamical system that uses a set of parameter-dependent terms to account for the effects of unresolved or neglected dynamics on resolved scales. We start by considering various methods of calibrating and analyzing such a model given a few well-resolved simulations. After presenting results for various point estimates and discussing some of their shortcomings, we demonstrate (a) the potential of hierarchical Bayesian analysis to uncover previously unanticipated physical dependencies in the approximate model, and (b) how such insights can then be used to improve the model. In effect parametric dependencies found from the Bayesian analysis are used to improve structural aspects of the model. While we choose to illustrate the procedure in the context of a closure model for buoyancy-driven, variable-density turbulence, the statistical nature of the approach makes it more generally applicable. Towards addressing issues of increased computational cost associated...
more | pdf | html
###### Tweets
StatsPapers: Leveraging Bayesian Analysis To Improve Accuracy of Approximate Models. https://t.co/b7D9krXD7N
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 3
Total Words: 12410
Unqiue Words: 3221

##### #2. Large fluctuations in non-ideal coarse-grained systems
###### M. Reza Parsa, Alexander J. Wagner
Using the recently introduced Molecular Dynamics Lattice Gas (MDLG) approach, we test fluctuations of coarse-grained quantities. We show that as soon as the system can no longer be considered an ideal gas fluctuations fail to diminish upon coarse-graining as is usually expected. These results suggest that current approaches to simulating fluctuating hydrodynamics may have to be augmented to achieve quantitative results for systems with a non-ideal equation of state. The MDLG method gives a guidance to the exact nature of the fluctuation in such systems.
more | pdf | html
None.
###### Tweets
CondensedPapers: Large fluctuations in non-ideal coarse-grained systems. https://t.co/slHAX8aKmE
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 3596
Unqiue Words: 1159

##### #3. First fully kinetic three-dimensional simulation of the AWAKE baseline scenario
###### N. Moschüring, K. V. Lotov, K. Bamberg, F. Deutschmann, H. Ruhl
The "Advanced Proton Driven Plasma Wakefield Acceleration Experiment" (AWAKE) aims to accelerate leptons via proton-beam-driven wakefield acceleration. It comprises extensive numerical studies as well as experiments at the CERN laboratory. The baseline scenario incorporates a plasma volume of approximately $62\,\mathrm{cm}^3$. The plasma wavelength is about $1.25\,\mathrm{mm}$ and needs to be adequately resolved, using a minimum of $130$ points per plasma wavelength, in order to accurately reproduce the physics. The baseline scenario incorporates the proton beam micro-bunching, the concurrent non-linear wakefield growth as well as the off-axis electron beam injection, trapping and acceleration. We present results for the first three-dimensional simulation of this baseline scenario with a full model, using a sufficient resolution. The simulation consumed about $22\,\mathrm{Mch}$ of computer resources and scaled up to 32768 cores, thanks to a multitude of adaptions, improvements and optimization of the simulation code PSC. Through...
more | pdf | html
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 5
Total Words: 6491
Unqiue Words: 2217

##### #4. Machine Learning of Algebraic Stress Models using Deterministic Symbolic Regression
###### Martin Schmelzer, Richard P. Dwight, Paola Cinnella
A novel deterministic symbolic regression method SpaRTA (Sparse Regression of Turbulent Stress Anisotropy) is introduced to infer algebraic stress models for the closure of RANS equations directly from high-fidelity LES or DNS data. The models are written as tensor polynomials and are built from a library of candidate functions. The machine-learning method is based on elastic net regularisation which promotes sparsity of the inferred models. By being data-driven the method relaxes assumptions commonly made in the process of model development. Model-discovery and cross-validation is performed for three cases of separating flows, i.e. periodic hills ($Re$=10595), converging-diverging channel ($Re$=12600) and curved backward-facing step ($Re$=13700). The predictions of the discovered models are significantly improved over the $k$-$\omega$ SST also for a true prediction of the flow over periodic hills at $Re$=37000. This study shows a systematic assessment of SpaRTA for rapid machine-learning of robust corrections for standard...
more | pdf | html
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 3
Total Words: 9219
Unqiue Words: 2381

##### #5. Building nonparametric $n$-body force fields using Gaussian process regression
###### Aldo Glielmo, Claudio Zeni, Ádám Fekete, Alessandro De Vita
Constructing a classical potential suited to simulate a given atomic system is a remarkably difficult task. This chapter presents a framework under which this problem can be tackled, based on the Bayesian construction of nonparametric force fields of a given order using Gaussian process (GP) priors. The formalism of GP regression is first reviewed, particularly in relation to its application in learning local atomic energies and forces. For accurate regression it is fundamental to incorporate prior knowledge into the GP kernel function. To this end, this chapter details how properties of smoothness, invariance and interaction order of a force field can be encoded into corresponding kernel properties. A range of kernels is then proposed, possessing all the required properties and an adjustable parameter $n$ governing the interaction order modelled. The order $n$ best suited to describe a given system can be found automatically within the Bayesian framework by maximisation of the marginal likelihood. The procedure is first tested on...
more | pdf | html
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 4
Total Words: 14377
Unqiue Words: 3511

##### #6. Exploiting Parallelism on Shared Memory in the QED Particle-in-Cell Code PICADOR with Greedy Load Balancing
###### Iosif Meyerov, Sergei Bastrakov, Aleksei Bashinov, Evgeny Efimenko, Alexander Panov, Elena Panova, Igor Surmin, Valentin Volokitin, Arkady Gonoskov
State-of-the-art numerical simulations of laser plasma by means of the Particle-in-Cell method are often extremely computationally intensive. Therefore there is a growing need for development of approaches for efficient utilization of resources of modern supercomputers. In this paper, we address the problem of a substantially non-uniform and dynamically varying distribution of macroparticles in a computational area in simulating quantum electrodynamic (QED) cascades. We propose and evaluate a load balancing scheme for shared memory systems, which allows subdividing individual cells of the computational domain into work portions with subsequent dynamic distribution of these portions between OpenMP threads. Computational experiments on 1D, 2D, and 3D QED simulations show that the proposed scheme outperforms the previously developed standard and custom schemes in the PICADOR code by 2.1 to 10 times when employing several Intel Cascade Lake CPUs.
more | pdf | html
None.
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 9
Total Words: 0
Unqiue Words: 0

##### #7. A new Fragile Points Method (FPM) in computational mechanics, based on the concepts of Point Stiffnesses and Numerical Flux Corrections
###### Leiting Dong, Tian Yang, Kailei Wang, Satya N. Atluri
In this paper, a new method, named the Fragile Points Method (FPM), is developed for computer modeling in engineering and sciences. In the FPM, simple, local, polynomial, discontinuous and Point-based trial and test functions are proposed based on randomly scattered points in the problem domain. The local discontinuous polynomial trial and test functions are postulated by using the Generalized Finite Difference method. These functions are only piece-wise continuous over the global domain. By implementing the Point-based trial and test functions into the Galerkin weak form, we define the concept of Point Stiffnesses as the contribution of each Point in the problem domain to the global stiffness matrix. However, due to the discontinuity of trial and test functions in the domain, directly using the Galerkin weak form leads to inconsistency. To resolve this, Numerical Flux Corrections, which are frequently used in Discontinuous Galerkin methods are further employed in the FPM. The resulting global stiffness matrix is symmetric and...
more | pdf | html
None.
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 4
Total Words: 0
Unqiue Words: 0

##### #8. Performance Enhancement for High-order Gas-kinetic Scheme Based on WENO-adaptive-order Reconstruction
###### Xing Ji, Kun Xu
High-order gas-kinetic scheme (HGKS) has been well-developed in the past years. Abundant numerical tests including hypersonic flow, turbulence, and aeroacoustic problems, have been used to validate its accuracy, efficiency, and robustness. However, there are still rooms for its further improvement. Firstly, the reconstruction in the previous scheme mainly achieves a third-order accuracy for the initial non-equilibrium states. At the same time, the equilibrium state in space and time in HGKS has to be reconstructed separately. Secondly, it is complicated to get reconstructed data at Gaussian points from the WENO-type method in high dimensions. For HGKS, besides the point-wise values at the Gaussian points it also requires the slopes in both normal and tangential directions of a cell interface. Thirdly, there exists visible spurious overshoot/undershoot at weak discontinuities from the previous HGKS with the standard WENO reconstruction. In order to overcome these difficulties, in this paper we use an improved reconstruction for...
more | pdf | html
None.
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

##### #9. Solving Elliptic Interface Problems with Jump Conditions on Cartesians Grids
###### Daniil Bochkov, Frederic Gibou
We present a simple numerical algorithm for solving elliptic equations where the diffusion coefficient, the source term, the solution and its flux are discontinuous across an irregular interface. The algorithm produces second-order accurate solutions and first-order accurate gradients in the $L^\infty$-norm on Cartesian grids. The condition number is bounded, regardless of the ratio of the diffusion constant and scales like that of the standard 5-point stencil approximation on a rectangular grid with no interface. Numerical examples are given in two and three spatial dimensions.
more | pdf | html
None.
None.
None.
###### Other stats
Sample Sizes : None.
Authors: 2
Total Words: 0
Unqiue Words: 0

Assert is a website where the best academic papers on arXiv (computer science, math, physics), bioRxiv (biology), BITSS (reproducibility), EarthArXiv (earth science), engrXiv (engineering), LawArXiv (law), PsyArXiv (psychology), SocArXiv (social science), and SportRxiv (sport research) bubble to the top each day.

Papers are scored (in real-time) based on how verifiable they are (as determined by their Github repos) and how interesting they are (based on Twitter).

To see top papers, follow us on twitter @assertpub_ (arXiv), @assert_pub (bioRxiv), and @assertpub_dev (everything else).

To see beautiful figures extracted from papers, follow us on Instagram.

Tracking 129,961 papers.

###### Search
Sort results based on if they are interesting or reproducible.
Interesting
Reproducible
Online
###### Stats
Tracking 129,961 papers.