A Comparison of Several Preconditioners in the Solution of Parabolic Systems in Three Space Dimensions Using DASPK and a High-Order Finite Element Method
نویسندگان
چکیده
We describe an algorithm for solving three-dimensional vector systems of parabolic partial diierential equations. Solutions are calculated using Galerkin's method with a high-order, piecewise polynomial hierarchical basis in space and the DAE code DASPK in time. A variant of DASPK that computes Jacobian vector products with a stored analytic Jacobian is compared with the default matrix-free approximation of the Jacobian vector product on three test problems. Several preconditioners are also compared including incomplete Cholesky, banded LU, and element-by-element factorizations.
منابع مشابه
A new positive definite semi-discrete mixed finite element solution for parabolic equations
In this paper, a positive definite semi-discrete mixed finite element method was presented for two-dimensional parabolic equations. In the new positive definite systems, the gradient equation and flux equations were separated from their scalar unknown equations. Also, the existence and uniqueness of the semi-discrete mixed finite element solutions were proven. Error estimates were also obtaine...
متن کاملMultilevel Schwarz and Multigrid preconditioners for the Bidomain system
Two parallel and scalable multilevel preconditioners for the Bidomain system in computational electrocardiology are introduced and studied. The Bidomain system, consisting of two degenerate parabolic reaction-diffusion equations coupled with a stiff system of several ordinary differential equations, generates very ill-conditioned discrete systems when discretized with semi-implicit methods in t...
متن کاملVARIATIONAL DISCRETIZATION AND MIXED METHODS FOR SEMILINEAR PARABOLIC OPTIMAL CONTROL PROBLEMS WITH INTEGRAL CONSTRAINT
The aim of this work is to investigate the variational discretization and mixed finite element methods for optimal control problem governed by semi linear parabolic equations with integral constraint. The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control is not discreted. Optimal error estimates in L2 are established for the state...
متن کاملکاربرد روش معادله سهموی در تحلیل مسائل انتشار امواج داخل ساختمان
With the rapid growth of indoor wireless communication systems, the need to accurately model radio wave propagation inside the building environments has increased. Many site-specific methods have been proposed for modeling indoor radio channels. Among these methods, the ray tracing algorithm and the finite-difference time domain (FDTD) method are the most popular ones. The ray tracing approach ...
متن کاملFree Vibration Analysis of a Sloping-frame: Closed-form Solution versus Finite Element Solution and Modification of the Characteristic Matrices (TECHNICAL NOTE)
This article deals with the free vibration analysis and determination of the seismic parameters of a sloping-frame which consists of three members; a horizontal, a vertical, and an inclined member. The both ends of the frame are clamped, and the members are rigidly connected at joint points. The individual members of the frame are assumed to be governed by the transverse vibration theory of an ...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2007