Using the Schur Complement to Reduce Runtime in KULL's Magnetic Diffusion Package
Recently a Resistive Magnetohydrodynamics (MHD) package has been added to the KULL code. In order to be compatible with the underlying hydrodynamics algorithm, a new sub-zonal magnetics discretization was developed that supports arbitrary polygonal and polyhedral zones. This flexibility comes at the cost of many more unknowns per zone - approximately ten times more for a hexahedral mesh. We can eliminate some (or all, depending on the dimensionality) of the extra unknowns from the global matrix during assembly by using a Schur complement approach. This trades expensive global work for cache-friendly local work, while still allowing solution for the full system. Significant improvements in the solution time are observed for several test problems.
- Research Organization:
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- W-7405-ENG-48
- OSTI ID:
- 1018784
- Report Number(s):
- LLNL-PROC-464197; TRN: US1103462
- Resource Relation:
- Conference: Presented at: Nuclear Explosives Code Developers' Conference 2010, Los Alamos, NM, United States, Oct 18 - Oct 22, 2010
- Country of Publication:
- United States
- Language:
- English
Similar Records
An analysis of options available for developing a common laser ray tracing package for Ares and Kull code frameworks
Using the Schur Complement to Reduce Runtime in KULL's Magnetic Diffusion Package