6,585 research outputs found
Finite Element Simulation of Dense Wire Packings
A finite element program is presented to simulate the process of packing and
coiling elastic wires in two- and three-dimensional confining cavities. The
wire is represented by third order beam elements and embedded into a
corotational formulation to capture the geometric nonlinearity resulting from
large rotations and deformations. The hyperbolic equations of motion are
integrated in time using two different integration methods from the Newmark
family: an implicit iterative Newton-Raphson line search solver, and an
explicit predictor-corrector scheme, both with adaptive time stepping. These
two approaches reveal fundamentally different suitability for the problem of
strongly self-interacting bodies found in densely packed cavities. Generalizing
the spherical confinement symmetry investigated in recent studies, the packing
of a wire in hard ellipsoidal cavities is simulated in the frictionless elastic
limit. Evidence is given that packings in oblate spheroids and scalene
ellipsoids are energetically preferred to spheres.Comment: 17 pages, 7 figures, 1 tabl
Method of lines transpose: High order L-stable O(N) schemes for parabolic equations using successive convolution
We present a new solver for nonlinear parabolic problems that is L-stable and
achieves high order accuracy in space and time. The solver is built by first
constructing a single-dimensional heat equation solver that uses fast O(N)
convolution. This fundamental solver has arbitrary order of accuracy in space,
and is based on the use of the Green's function to invert a modified Helmholtz
equation. Higher orders of accuracy in time are then constructed through a
novel technique known as successive convolution (or resolvent expansions).
These resolvent expansions facilitate our proofs of stability and convergence,
and permit us to construct schemes that have provable stiff decay. The
multi-dimensional solver is built by repeated application of dimensionally
split independent fundamental solvers. Finally, we solve nonlinear parabolic
problems by using the integrating factor method, where we apply the basic
scheme to invert linear terms (that look like a heat equation), and make use of
Hermite-Birkhoff interpolants to integrate the remaining nonlinear terms. Our
solver is applied to several linear and nonlinear equations including heat,
Allen-Cahn, and the Fitzhugh-Nagumo system of equations in one and two
dimensions
Scalable explicit implementation of anisotropic diffusion with Runge-Kutta-Legendre super-time-stepping
An important ingredient in numerical modelling of high temperature magnetised
astrophysical plasmas is the anisotropic transport of heat along magnetic field
lines from higher to lower temperatures.Magnetohydrodynamics (MHD) typically
involves solving the hyperbolic set of conservation equations along with the
induction equation. Incorporating anisotropic thermal conduction requires to
also treat parabolic terms arising from the diffusion operator. An explicit
treatment of parabolic terms will considerably reduce the simulation time step
due to its dependence on the square of the grid resolution () for
stability. Although an implicit scheme relaxes the constraint on stability, it
is difficult to distribute efficiently on a parallel architecture. Treating
parabolic terms with accelerated super-time stepping (STS) methods has been
discussed in literature but these methods suffer from poor accuracy (first
order in time) and also have difficult-to-choose tuneable stability parameters.
In this work we highlight a second order (in time) Runge Kutta Legendre (RKL)
scheme (first described by Meyer et. al. 2012) that is robust, fast and
accurate in treating parabolic terms alongside the hyperbolic conversation
laws. We demonstrate its superiority over the first order super time stepping
schemes with standard tests and astrophysical applications. We also show that
explicit conduction is particularly robust in handling saturated thermal
conduction. Parallel scaling of explicit conduction using RKL scheme is
demonstrated up to more than processors.Comment: 15 pages, 9 figures, incorporated comments from the referee. This
version is now accepted for publication in MNRA
Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms
Some numerical aspects of finite-difference algorithms for nonlinear multidimensional hyperbolic conservation laws with stiff nonhomogenous (source) terms are discussed. If the stiffness is entirely dominated by the source term, a semi-implicit shock-capturing method is proposed provided that the Jacobian of the soruce terms possesses certain properties. The proposed semi-implicit method can be viewed as a variant of the Bussing and Murman point-implicit scheme with a more appropriate numerical dissipation for the computation of strong shock waves. However, if the stiffness is not solely dominated by the source terms, a fully implicit method would be a better choice. The situation is complicated by problems that are higher than one dimension, and the presence of stiff source terms further complicates the solution procedures for alternating direction implicit (ADI) methods. Several alternatives are discussed. The primary motivation for constructing these schemes was to address thermally and chemically nonequilibrium flows in the hypersonic regime. Due to the unique structure of the eigenvalues and eigenvectors for fluid flows of this type, the computation can be simplified, thus providing a more efficient solution procedure than one might have anticipated
Adaptive multiresolution computations applied to detonations
A space-time adaptive method is presented for the reactive Euler equations
describing chemically reacting gas flow where a two species model is used for
the chemistry. The governing equations are discretized with a finite volume
method and dynamic space adaptivity is introduced using multiresolution
analysis. A time splitting method of Strang is applied to be able to consider
stiff problems while keeping the method explicit. For time adaptivity an
improved Runge--Kutta--Fehlberg scheme is used. Applications deal with
detonation problems in one and two space dimensions. A comparison of the
adaptive scheme with reference computations on a regular grid allow to assess
the accuracy and the computational efficiency, in terms of CPU time and memory
requirements.Comment: Zeitschrift f\"ur Physicalische Chemie, accepte
- …