An adaptive moving grid method for one-dimensional systems of partial differential equations and its numerical solution
In this paper we examine a scheme for choosing a moving mesh based on minimizing the time rate of change of the solution in the moving coordinates. We show how to apply this method to systems where the time derivatives cannot be solved for explicitly, writing the moving mesh equations in an implicit form. We give a geometrical interpretation of the moving mesh equation which exposes some of its weaknesses, and suggest some modifications based on this interpretation which increase the efficiency of the scheme. Upon discretization of the original PDE system coupled with the moving mesh equations, a system of differential/algebraic equations (DAEs) is generated. We discuss the solution of this system using the code DASSL, and show how changing the error estimate for time step selection from the usual strategy can significantly reduce the computation time. We introduce a new DAE code, DASPK, which combines the time integration methods of DASSL with the preconditioned GMRES method for solving large sparse linear systems. Finally, we present the results of some numerical experiments on reaction-diffusion equations which illustrate how well the resulting methods work. 16 refs., 2 figs., 2 tabs.
- Research Organization:
- Lawrence Livermore National Lab., CA (USA)
- DOE Contract Number:
- W-7405-ENG-48
- OSTI ID:
- 6487746
- Report Number(s):
- UCRL-100289; CONF-8810307-1; ON: DE89006743
- Resource Relation:
- Conference: Workshop on adaptive methods for partial differential equations, Troy, NY, USA, 13 Oct 1988; Other Information: Portions of this document are illegible in microfiche products
- Country of Publication:
- United States
- Language:
- English
Similar Records
Using Krylov methods in the solution of large-scale differential-algebraic systems
Marker ReDistancing/Level Set Method for High-Fidelity Implicit Interface Tracking