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
StatsPapers:
Leveraging Bayesian Analysis To Improve Accuracy of Approximate Models. https://t.co/b7D9krXD7N

None.

None.

Sample Sizes : None.

Authors: 3

Total Words: 12410

Unqiue Words: 3221

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.

CondensedPapers:
Large fluctuations in non-ideal coarse-grained systems. https://t.co/slHAX8aKmE

None.

None.

Sample Sizes : None.

Authors: 2

Total Words: 3596

Unqiue Words: 1159

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.

Sample Sizes : None.

Authors: 5

Total Words: 6491

Unqiue Words: 2217

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.

Sample Sizes : None.

Authors: 3

Total Words: 9219

Unqiue Words: 2381

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.

Sample Sizes : None.

Authors: 4

Total Words: 14377

Unqiue Words: 3511

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.

Sample Sizes : None.

Authors: 9

Total Words: 0

Unqiue Words: 0

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.

Sample Sizes : None.

Authors: 4

Total Words: 0

Unqiue Words: 0

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.

Sample Sizes : None.

Authors: 2

Total Words: 0

Unqiue Words: 0

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.

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.*

Sort results based on if they are interesting or reproducible.

Interesting

Reproducible