Our ability to numerically model and understand the complex flow behavior of
solid-bearing suspensions has increased significantly over the last couple of
years, partly due to direct numerical simulations that compute flow around
individual interfaces and hence resolve unprecedented detail. While most
previous studies focus on spherical particles, we develop a direct numerical
approach to capture rectangular particles. Our approach uses distributed
Lagrange multipliers to enforce rigid-body motion in the solid domain in
combination with an immersed boundary method to correctly enforce the no-slip
constraint on the solid-fluid interfaces. An important component of our model
is an efficient particle collision scheme that prevents overlap between
particles of different shapes and allows for the transfer of both translational
and angular momentum during particle collision. We verify and validate our
numerical method through several benchmark cases. Applied to suspension flow,
we test the hypothesize that particle rotations alter the...

more |
pdf
| html
PhysicsPaper:
Rotation matters - Direct numerical simulations of rectangular particles in suspensions at low to intermediate solid fraction. https://t.co/eglEtDuxOg

None.

None.

Sample Sizes : None.

Authors: 3

Total Words: 13995

Unqiue Words: 3483

We present a generalized algorithm based on a spherical harmonics expansion
method for efficient computation of the three-dimensional gravitational
potential on a multi-patch grid in spherical geometry. Instead of solving for
the gravitational potential by superposition of separate contributions from the
mass density distribution on individual grid patch our new algorithm computes
directly the gravitational potential due to contributions from all grid patches
in one computation step, thereby reducing the computational cost of the gravity
solver. This is possible by considering a set of angular weights which are
derived from rotations of spherical harmonics functions defined in a global
coordinate system that is common for all grid patches. Additionally, our
algorithm minimizes data communication between parallel compute tasks by
eliminating its proportionality to the number of subdomains in the grid
configuration, making it suitable for parallelized computation on a multi-patch
grid configuration with any number of subdomains....

more |
pdf
| html
None.

None.

Sample Sizes : None.

Authors: 1

Total Words: 9180

Unqiue Words: 2241

Boundary integral equation methods for analyzing electromagnetic scattering
phenomena typically suffer from several of the following problems: (i)
ill-conditioning when the frequency is low; (ii) ill-conditioning when the
discretization density is high; (iii) ill-conditioning when the structure
contains global loops (which are computationally expensive to detect); (iv)
incorrect solution at low frequencies due to current cancellations; (v)
presence of spurious resonances. In this paper, quasi-Helmholtz projectors are
leveraged to obtain a magnetic field integral equation (MFIE) formulation that
is immune to drawbacks (i)-(iv). Moreover, when this new MFIE is combined with
a regularized electric field integral equation, a new quasi-Helmholtz projector
combined field integral equation is obtained that also is immune to (v).
Numerical results corroborate the theory and show the practical impact of the
newly proposed formulations.

more |
pdf
| html
PhysicsPaper:
Magnetic and Combined Field Integral Equations Based on the Quasi-Helmholtz Projectors. https://t.co/agTIv5gSdw

None.

None.

Sample Sizes : None.

Authors: 5

Total Words: 7864

Unqiue Words: 2078

Kinetic Monte-Carlo (KMC) simulations are a well-established numerical tool
to investigate the time-dependent surface morphology in molecular beam epitaxy
(MBE) experiments. One drawback of the KMC method, however, is the high
computational cost at growth conditions, where particle diffusion dominates
over particle adsorption. As an alternative, we introduce in this study a
limited mobility (LM) model to realistically simulate the low-temperature
regime (where detachment processes can be neglected). The model is based on the
stochastic transition rules of the so-called Das Sarma-Tamborena (DT) model,
but differs from the latter via a variable diffusion length. A first guess for
this length can be extracted from the mean-squared displacement calculated from
short KMC simulations. Comparing the resulting surface morphologies in the
submonolayer growth regime to those obtained from KMC simulations, we find
deviations which can be cured by adding fluctuations to the diffusion length.
This mimics the stochastic nature of particle...

more |
pdf
| html
None.

None.

Sample Sizes : None.

Authors: 2

Total Words: 10112

Unqiue Words: 2327

Surrogate models for hotspot ignition and growth rates were presented in Part
I, where the hotspots were formed by the collapse of single cylindrical voids.
Such isolated cylindrical voids are idealizations of the void morphology in
real meso-structures. This paper therefore investigates the effect of
non-cylindrical void shapes and void-void interactions on hotspot ignition and
growth. Surrogate models capturing these effects are constructed using a
Bayesian Kriging approach. The training data for machine learning the
surrogates are derived from reactive void collapse simulations spanning the
parameter space of void aspect ratio (AR), void orientation (\Theta), and void
fraction (\phi). The resulting surrogate models portray strong dependence of
the ignition and growth rates on void aspect ratio and orientation,
particularly when they are oriented at acute angles with respect to the imposed
shock. The surrogate models for void interaction effects show significant
changes in hotspot ignition and growth rates as the void fraction...

more |
pdf
| html
None.

None.

Sample Sizes : None.

Authors: 6

Total Words: 11868

Unqiue Words: 2604

Convection is an important physical process in astrophysics well-studied
using numerical simulations under the Boussinesq and/or anelastic
approximations. However these approaches reach their limits when compressible
effects are important in the high Mach flow regime, e.g. in stellar atmospheres
or in the presence of accretion shocks. In order to tackle these issues, we
propose a new high performance and portable code, called "ARK" with a numerical
solver well-suited for the stratified compressible Navier-Stokes equations. We
take a finite volume approach with machine precision conservation of mass,
transverse momentum and total energy. Based on previous works in applied
mathematics we propose the use of a low Mach correction to achieve a good
precision in both low and high Mach regimes. The gravity source term is
discretized using a well-balanced scheme in order to reach machine precision
hydrostatic balance. This new solver is implemented using the Kokkos library in
order to achieve high performance computing and portability...

more |
pdf
| html
None.

None.

Sample Sizes : None.

Authors: 5

Total Words: 9716

Unqiue Words: 2488

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 99,586 papers.*

Sort results based on if they are interesting or reproducible.

Interesting

Reproducible