9,116 research outputs found
Affordable, Entropy Conserving and Entropy Stable Flux Functions for the Ideal MHD Equations
In this work, we design an entropy stable, finite volume approximation for
the ideal magnetohydrodynamics (MHD) equations. The method is novel as we
design an affordable analytical expression of the numerical interface flux
function that discretely preserves the entropy of the system. To guarantee the
discrete conservation of entropy requires the addition of a particular source
term to the ideal MHD system. Exact entropy conserving schemes cannot dissipate
energy at shocks, thus to compute accurate solutions to problems that may
develop shocks, we determine a dissipation term to guarantee entropy stability
for the numerical scheme. Numerical tests are performed to demonstrate the
theoretical findings of entropy conservation and robustness.Comment: arXiv admin note: substantial text overlap with arXiv:1509.06902;
text overlap with arXiv:1007.2606 by other author
Entropy Stable Finite Volume Approximations for Ideal Magnetohydrodynamics
This article serves as a summary outlining the mathematical entropy analysis
of the ideal magnetohydrodynamic (MHD) equations. We select the ideal MHD
equations as they are particularly useful for mathematically modeling a wide
variety of magnetized fluids. In order to be self-contained we first motivate
the physical properties of a magnetic fluid and how it should behave under the
laws of thermodynamics. Next, we introduce a mathematical model built from
hyperbolic partial differential equations (PDEs) that translate physical laws
into mathematical equations. After an overview of the continuous analysis, we
thoroughly describe the derivation of a numerical approximation of the ideal
MHD system that remains consistent to the continuous thermodynamic principles.
The derivation of the method and the theorems contained within serve as the
bulk of the review article. We demonstrate that the derived numerical
approximation retains the correct entropic properties of the continuous model
and show its applicability to a variety of standard numerical test cases for
MHD schemes. We close with our conclusions and a brief discussion on future
work in the area of entropy consistent numerical methods and the modeling of
plasmas
An entropy stable discontinuous Galerkin method for the shallow water equations on curvilinear meshes with wet/dry fronts accelerated by GPUs
We extend the entropy stable high order nodal discontinuous Galerkin spectral
element approximation for the non-linear two dimensional shallow water
equations presented by Wintermeyer et al. [N. Wintermeyer, A. R. Winters, G. J.
Gassner, and D. A. Kopriva. An entropy stable nodal discontinuous Galerkin
method for the two dimensional shallow water equations on unstructured
curvilinear meshes with discontinuous bathymetry. Journal of Computational
Physics, 340:200-242, 2017] with a shock capturing technique and a positivity
preservation capability to handle dry areas. The scheme preserves the entropy
inequality, is well-balanced and works on unstructured, possibly curved,
quadrilateral meshes. For the shock capturing, we introduce an artificial
viscosity to the equations and prove that the numerical scheme remains entropy
stable. We add a positivity preserving limiter to guarantee non-negative water
heights as long as the mean water height is non-negative. We prove that
non-negative mean water heights are guaranteed under a certain additional time
step restriction for the entropy stable numerical interface flux. We implement
the method on GPU architectures using the abstract language OCCA, a unified
approach to multi-threading languages. We show that the entropy stable scheme
is well suited to GPUs as the necessary extra calculations do not negatively
impact the runtime up to reasonably high polynomial degrees (around ). We
provide numerical examples that challenge the shock capturing and positivity
properties of our scheme to verify our theoretical findings
A boundary integral formalism for stochastic ray tracing in billiards
Determining the flow of rays or non-interacting particles driven by a force or velocity field is fundamental to modelling many physical processes. These include particle flows arising in fluid mechanics and ray flows arising in the geometrical optics limit of linear wave equations. In many practical applications, the driving field is not known exactly and the dynamics are determined only up to a degree of uncertainty. This paper presents a boundary integral framework for propagating flows including uncertainties, which is shown to systematically interpolate between a deterministic and a completely random description of the trajectory propagation. A simple but efficient discretisation approach is applied to model uncertain billiard dynamics in an integrable rectangular domain
An Entropy Stable Nodal Discontinuous Galerkin Method for the Two Dimensional Shallow Water Equations on Unstructured Curvilinear Meshes with Discontinuous Bathymetry
We design an arbitrary high-order accurate nodal discontinuous Galerkin
spectral element approximation for the nonlinear two dimensional shallow water
equations with non-constant, possibly discontinuous, bathymetry on
unstructured, possibly curved, quadrilateral meshes. The scheme is derived from
an equivalent flux differencing formulation of the split form of the equations.
We prove that this discretisation exactly preserves the local mass and
momentum. Furthermore, combined with a special numerical interface flux
function, the method exactly preserves the mathematical entropy, which is the
total energy for the shallow water equations. By adding a specific form of
interface dissipation to the baseline entropy conserving scheme we create a
provably entropy stable scheme. That is, the numerical scheme discretely
satisfies the second law of thermodynamics. Finally, with a particular
discretisation of the bathymetry source term we prove that the numerical
approximation is well-balanced. We provide numerical examples that verify the
theoretical findings and furthermore provide an application of the scheme for a
partial break of a curved dam test problem
The effects of short-lived radionuclides and porosity on the early thermo-mechanical evolution of planetesimals
The thermal history and internal structure of chondritic planetesimals,
assembled before the giant impact phase of chaotic growth, potentially yield
important implications for the final composition and evolution of terrestrial
planets. These parameters critically depend on the internal balance of heating
versus cooling, which is mostly determined by the presence of short-lived
radionuclides (SLRs), such as aluminum-26 and iron-60, as well as the heat
conductivity of the material. The heating by SLRs depends on their initial
abundances, the formation time of the planetesimal and its size. It has been
argued that the cooling history is determined by the porosity of the granular
material, which undergoes dramatic changes via compaction processes and tends
to decrease with time. In this study we assess the influence of these
parameters on the thermo-mechanical evolution of young planetesimals with both
2D and 3D simulations. Using the code family I2ELVIS/I3ELVIS we have run
numerous 2D and 3D numerical finite-difference fluid dynamic models with
varying planetesimal radius, formation time and initial porosity. Our results
indicate that powdery materials lowered the threshold for melting and
convection in planetesimals, depending on the amount of SLRs present. A subset
of planetesimals retained a powdery surface layer which lowered the thermal
conductivity and hindered cooling. The effect of initial porosity was small,
however, compared to those of planetesimal size and formation time, which
dominated the thermo-mechanical evolution and were the primary factors for the
onset of melting and differentiation. We comment on the implications of this
work concerning the structure and evolution of these planetesimals, as well as
their behavior as possible building blocks of terrestrial planets.Comment: 19 pages, 11 figures, 5 tables; accepted for publication in Icarus;
for associated video files, see http://timlichtenberg.net/2015_porosity.html
or http://dx.doi.org/10.1016/j.icarus.2016.03.00
A Provably Stable Discontinuous Galerkin Spectral Element Approximation for Moving Hexahedral Meshes
We design a novel provably stable discontinuous Galerkin spectral element
(DGSEM) approximation to solve systems of conservation laws on moving domains.
To incorporate the motion of the domain, we use an arbitrary
Lagrangian-Eulerian formulation to map the governing equations to a fixed
reference domain. The approximation is made stable by a discretization of a
skew-symmetric formulation of the problem. We prove that the discrete
approximation is stable, conservative and, for constant coefficient problems,
maintains the free-stream preservation property. We also provide details on how
to add the new skew-symmetric ALE approximation to an existing discontinuous
Galerkin spectral element code. Lastly, we provide numerical support of the
theoretical results
Time-frequency analysis of ship wave patterns in shallow water: modelling and experiments
A spectrogram of a ship wake is a heat map that visualises the time-dependent
frequency spectrum of surface height measurements taken at a single point as
the ship travels by. Spectrograms are easy to compute and, if properly
interpreted, have the potential to provide crucial information about various
properties of the ship in question. Here we use geometrical arguments and
analysis of an idealised mathematical model to identify features of
spectrograms, concentrating on the effects of a finite-depth channel. Our
results depend heavily on whether the flow regime is subcritical or
supercritical. To support our theoretical predictions, we compare with data
taken from experiments we conducted in a model test basin using a variety of
realistic ship hulls. Finally, we note that vessels with a high aspect ratio
appear to produce spectrogram data that contains periodic patterns. We can
reproduce this behaviour in our mathematical model by using a so-called
two-point wavemaker. These results highlight the role of wave interference
effects in spectrograms of ship wakes.Comment: 14 pages, 7 figure
Dipolar spin correlations in classical pyrochlore magnets
We study spin correlations for the highly frustrated classical pyrochlore
lattice antiferromagnets with O(N) symmetry in the limit T->0. We conjecture
that a local constraint obeyed by the extensively degenerate ground states
dictates a dipolar form for the asymptotic spin correlations, at all N 2
for which the system is paramagnetic down to T=0. We verify this conjecture in
the cases N=1 and N=3 by simulations and to all orders in the 1/N expansion
about the solvable N=infinity limit. Remarkably, the N=infinity formulae are an
excellent fit, at all distances, to the correlators at N=3 and even at N=1.
Thus we obtain a simple analytical expression also for the correlations of the
equivalent models of spin ice and cubic water ice, I_h.Comment: 4 pages revtex
- …