6,585 research outputs found

    Finite Element Simulation of Dense Wire Packings

    Full text link
    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

    Get PDF
    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

    Full text link
    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 (Δx\Delta x) 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 10410^4 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

    Get PDF
    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

    Full text link
    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
    corecore