A New Stabilization of Adaptive Step Trapezoid Rule Based on Finite Difference Interrupts

The adaptive step trapezoid rule (TR) is a generally effective numerical integrator, but it is prone to ringing instability and solution stall. We introduce a new stabilization algorithm based on finite difference interrupts (FDI) that reduces ringing and prevents stall. Unlike previously reported stabilization schemes, our algorithm achieves stability and second order accuracy without incurring significant computational cost or spurious diffusion. TR with FDI is at least as stable and accurate as existing methods when solving the prototypical scalar problem. We demonstrate that it has better stability, accuracy, and cost when solving the tightly coupled system describing free surface flows of viscoelastic liquids. Though we demonstrate TR with FDI in the context of a finite element method-of-lines, it is applicable to any TR-based algorithm.

[1]  Charlie H. Cooke.  A characterization of stiffly stable linear multistep methods , 1973 .

[2]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947, Mathematical Proceedings of the Cambridge Philosophical Society.

[3]  Rolf Rannacher,et al.  On the smoothing property of the crank-nicolson scheme , 1982 .

[4]  Scott D. Phillips,et al.  Computational and experimental analysis of dynamics of drop formation , 1999 .

[5]  Ole Østerby.  Five Ways of Reducing the Crank–Nicolson Oscillations , 2003 .

[6]  M. Pasquali,et al.  Free surface flows of polymer solutions with models based on the conformation tensor , 2002 .

[7]  J. C. Simo,et al.  Unconditional stability and long-term behavior of transient algorithms for the incompressible Navier-Stokes and Euler equations , 1994 .

[8]  Wulf G. Dettmer,et al.  An analysis of the time integration algorithms for the finite element solutions of incompressible Navier-Stokes equations based on a stabilised formulation , 2003 .

[9]  P. Kloucek,et al.  Stability of the fractional step T-scheme for the nonstationary Navier-Stokes equations , 1994 .

[10]  Gareth H. McKinley,et al.  Formation of beads-on-a-string structures during break-up of viscoelastic filaments , 2010 .

[11]  G. Dahlquist.  A special stability problem for linear multistep methods , 1963 .

[12]  P. A. Sackinger,et al.  A finite element method for free surface flows of incompressible fluids in three dimensions. Part I. Boundary fitted mesh motion , 2000 .

[13]  Miroslav Grmela,et al.  Conformation tensor rheological models , 1987 .

[14]  Ole Østerby,et al.  Some numerical investigations of the stability of electrochemical digital simulation, particularly as affected by first-order homogeneous reactions , 1994 .

[15]  J. Rothstein,et al.  Computational study of viscoelastic effects on liquid transfer during gravure printing , 2013 .

[16]  J. Szmelter.  Incompressible flow and the finite element method , 2001 .

[17]  Sandra M. C. Malta,et al.  Avoiding spurious modes of time discretized operators in transport problems , 2010 .

[18]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[19]  Randolph E. Bank,et al.  Transient simulation of silicon devices and circuits , 1985 .

[20]  Catherine E. Powell,et al.  A framework for the development of implicit solvers for incompressible flow problems , 2012 .

[21]  Bamin Khomami,et al.  The influence of polymer concentration and chain architecture on free surface displacement flows of polymeric fluids , 2005 .

[22]  Lawrence F. Shampine,et al.  Analysis and implementation of TR-BDF2 , 1996 .

[23]  J. M. Watt.  Numerical Initial Value Problems in Ordinary Differential Equations , 1972 .

[24]  S. Turek.  a Comparative Study of Time-Stepping Techniques for the Incompressible Navier-Stokes Equations: from Fully Implicit Non-Linear Schemes to Semi-Implicit Projection Methods , 1996 .

[25]  Bamin Khomami,et al.  A computational study of the influence of viscoelasticity on the interfacial dynamics of dip coating , 2011 .

[26]  Suhas V. Patankar,et al.  A NEW FINITE-DIFFERENCE SCHEME FOR PARABOLIC DIFFERENTIAL EQUATIONS , 1978 .

[27]  C. W. Gear,et al.  The automatic integration of stiff ordinary differential equations , 1968, IFIP Congress.

[28]  G. Strang,et al.  Optimal stability for trapezoidal–backward difference split-steps , 2010 .

[29]  Howard C. Elman,et al.  Fast iterative solvers for buoyancy driven flow problems , 2011, J. Comput. Phys..

[30]  Carl E. Pearson,et al.  Impulsive end condition for diffusion equation , 1965 .

[31]  Jeff Cash.  Split linear multistep methods for the numerical integration of stiff differential systems , 1983 .

[32]  R. Glowinski,et al.  Viscous Flow Simulation by Finite Element Methods and Related Numerical Techniques , 1985 .

[33]  L. E. Scriven,et al.  Flows in Forward Deformable Roll Coating Gaps , 1997 .

[34]  N. J. Salamon,et al.  A self-optimizing single-step algorithm for the numerical integration of parabolic equations , 1995 .

[35]  Uri M. Ascher,et al.  Computer methods for ordinary differential equations and differential-algebraic equations , 1998 .

[36]  Bamin Khomami,et al.  Time-dependent simulations of non-axisymmetric patterns in Taylor-Couette flow of dilute polymer solutions , 2006 .

[37]  P. J. Dotson,et al.  Polymer solution rheology based on a finitely extensible bead—spring chain model , 1980 .

[38]  L. E. Scriven,et al.  Theoretical modeling of microstructured liquids: a simple thermodynamic approach , 2004 .

[39]  David F. Griffiths,et al.  Adaptive Time-Stepping for Incompressible Flow Part I: Scalar Advection-Diffusion , 2008, SIAM J. Sci. Comput..

[40]  K. Morgan,et al.  Recent advances in numerical methods in fluids , 1980 .

[41]  B. Lindberg.  On smoothing and extrapolation for the trapezoidal rule , 1971 .

[42]  Roland W. Lewis,et al.  A comparison of time marching schemes for the transient heat conduction equation , 1975 .

[43]  Ole Østerby,et al.  Damping of Crank-Nicolson error oscillations , 2003, Comput. Biol. Chem..

[44]  L. E. Scriven,et al.  Discretization of free surface flows and other moving boundary problems , 1992 .

[45]  David F. Griffiths,et al.  Adaptive Time-Stepping for Incompressible Flow Part II: Navier--Stokes Equations , 2010, SIAM J. Sci. Comput..

[46]  Jaroslav Hron,et al.  Numerical Study of a Modified Time-Stepping θ-Scheme for Incompressible Flow Simulations , 2006, J. Sci. Comput..

[47]  Barry Lee,et al.  Finite elements and fast iterative solvers: with applications in incompressible fluid dynamics , 2006, Math. Comput..

[48]  Matteo Pasquali,et al.  Beads-on-string formation during filament pinch-off: Dynamics with the PTT model for non-affine motion , 2009 .

[49]  Pradeep P. Bhat,et al.  Dynamics of viscoelastic liquid filaments: low capillary number flows , 2008 .

[50]  Oldrich Joel Romero,et al.  Response of slot coating flows to periodic disturbances , 2008 .

[51]  Oreste S. Bursi,et al.  The analysis of the Generalized-α method for non-linear dynamic problems , 2002 .

[52]  Gregory M. Hulbert,et al.  The unimportance of the spurious root of time integration algorithms for structural dynamics , 1994 .

[53]  Joonas Sorvari,et al.  Time integration in linear viscoelasticity—a comparative study , 2010 .

[54]  L. Shampine,et al.  Computer solution of ordinary differential equations : the initial value problem , 1975 .

[55]  Milan D. Mihajlovic,et al.  Unstructured finite element method for the solution of the Boussinesq problem in three dimensions , 2013 .

[56]  C F Curtiss,et al.  Integration of Stiff Equations. , 1952, Proceedings of the National Academy of Sciences of the United States of America.

[57]  L. E. Scriven,et al.  Slot coating of mildly viscoelastic liquids , 2006 .

[58]  Jintai Chung,et al.  A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method , 1993 .