A Fast Matrix Compression Method for Large Scale Numerical Modelling of Rotationally Symmetric 3D Passive Structures in Fusion Devices
Abstract
:1. Introduction
- The extension of the method when both the structure and the sources exhibit the same symmetries. This is quite important because symmetry is typical in devices for fusion applications.
- A new approach (namely the DOF-based method) is introduced for the QR-recursive method, which compared to the old one (ELEMENT-based method) is numerically more effective.
- An efficient numerical approach is used in solving the system in order to handle the “electrodes” case. This case is very usual for practical fusion device application.
- We tackle the problem of the “small boxes” in QR-recursive, which could degrade the performances of the overall method.
2. Mathematical Formulation
2.1. Integral Formulation
2.2. Numerical Model
- is a column vector discretizing the unknown density current (i.e., )
- is a column vector representing the unknown voltages at the electrodes (i.e., )
- is the contribution to the current flowing in the electrode i due to the DOF j
- is zero for those DOF j, which do not belong to the boundary.
- is the time step
- and are the unknowns at integration instant
- are the prescribed currents flowing in the electrodes at integration instant
- and are source terms at integration instant
2.3. Symmetric Periodic Boundary Condition
- (i)
- The assembly of has a cost of .
- (ii)
- The inversion (usually factorization) of the matrix , using a direct method, as well know, requires operations.
3. QR-Recursive Method
3.1. Summary of QR-Recursive Approach
- (1)
- Boxes generation
- (2)
- Definition of DOFs associated to the boxes
- (3)
- Definition of interaction matrix between two boxes
- (4)
- Setting Lnear, Lfar
- (5)
- Overall compressed operator optimization and practical considerations
3.1.1. Boxes Generation
- (a)
- The cell is able to tessellate the 3D space.
- (b)
- The cell can be exactly subdivided in eight smaller replicas.
- (c)
- The cell should present an aspect ratio of the order of unity. In other words, the cell should not be either flattened or elongated.
3.1.2. Definition of DOFs Associated to the Boxes
3.1.3. Definition of the Interaction Matrix between Two Boxes
- for each matrix entry of for which i or j is a boundary DOF.
- for each matrix entry of for which both i and j are internal DOFs.
3.1.4. Setting
- is the set of sources boxes c, such that the box is in far zone of the box c.
- is the set of sources boxes c, such that the box is in the near of the box c.
- When the two boxes are classified as near, hence, they have a large rank and cannot be efficiently compressed.
- The interactions arising from non-local DOFs. These DOFs could not be compressed because they are related to nonlocal shape functions.
- (i)
- assembling of
- (ii)
- product evaluation of (i.e., )
3.2. The Preconditioner and the Initial Guess Estimation
- The computational cost (memory and time) of the factorization is linear versus the number of unknowns. Actually, it uses Cholesky decomposition (tailored for sparse matrix).
- The back/forward substitution (involved in the Cholesky method) is also very cheap.
3.3. Extension of Compression to Symmetries
- Reducing the solution domain only to an elementary part of the whole structure. In the following, this part of the domain is called the main sector.
- Assuming that basis functions automatically verify the symmetry conditions, by two suitable operators: reflection and rotation.
- (i)
- The number of DOFs is reduced by a factor of . This gives a huge gain in matrix storing, factorizing, and inverting.
- (ii)
- The integration (43) can be seen as limiting the outer integral defined in (12) only to the main sector. Clearly, this reduces the matrix assembly time.
- (1)
- Compressing each single interaction matrix appearing in (47).
- (2)
- Summing the interactions matrices and after compressing the resulting matrix .
3.4. DOF-Based Approach
- (1)
- In the Boxes Generation Algorithm, we assume the object to be individual DOFs.
- (2)
- In Equation (24) we change the definition of the set as follows:
- time required to assembly the compressed operator L,
- the compression versus accuracy, signature.
3.5. Handling Small Size Boxes
3.6. Handling Electrodes
- 1.
- Apply an iterative method to solve the augmeted system
- 2.
- Elimining by substituing the first equation of the system (22) into the second one. This subsitution yields the resulting linear system
4. Results
4.1. Testcase #1
4.2. Testcase #2
5. Discussion and Conclusions
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Acknowledgments
Conflicts of Interest
Appendix A
Algorithm A1. Boxes Generation |
Input |
Choose the “objects” (elements or DOF) |
smin1 is the minimum number of objects in a box |
Output |
BoxSet, FatherSet, ChildrenSet |
is the set of all not empty boxes |
is the set of all box having a number of objects number of objects less equal to smin1 |
returns the father of the box b |
returns the eight children of the box b |
Create first box b0 // the smallest cube containing the mesh |
level = 0 // current level starting from 0 |
// temporary |
// temporary |
While // Loop util there are still fathers |
level = level + 1 // current level |
//current level children list |
For // loop over previous level father box F |
// set the number of objects in box F |
If |
// Subdivide father F in eight children C |
For = 1,8 |
C = // get i-th children C |
// get the number of objects in box B |
If |
// Update Sets |
// C is the son of F |
// F is the father of C |
// C will be next level fathers |
End if |
End for |
Else |
// father F is childless |
End if |
End while |
// Generation Update: current sons will be next level fathers |
End Algorithm |
Appendix B
- is the set of boxes made by and all childless boxes adjacent to .
- is the set of boxes not adjacent to at the same level of and are well separated by b.
- , for a childless box b, is the set of descendent of b’s colleagues that are not adjacent to b, but whose parent boxes are adjacent to b.
- is the set of c such that .
Appendix C
Algorithm A2. Small Boxes Fusion |
Input |
is the minimum number of objects |
is a parameter controlling maximum distance |
Output |
is the set of objects in the box |
// is the set of the boxes to be retained |
//is the set of the boxes to be fused |
For // loop over all childless boxes |
// get number of objects in the box b |
If smin2 |
// put b in list of boxes to be retained |
Else |
// put b into list of boxes to be fused |
End if |
End for |
For // loop over the boxes to be fused |
//compute the distance from b |
For |
d(x) = distance(b,x) // distance between b and x |
End |
// Compute the minimum distance |
// set the maximum allowed distance |
// Create the set of all boxes whose distance from b is less than max_d |
// get the number of objects for this set |
For |
End for |
// Choose in the set Nearest, the one having the minimum number of objects |
// fuse all objects of the box b into the z box |
// the objects of the box b must be deleted |
End for |
End of the algorithm |
References
- Villone, F.; Albanese, R.; Liu, Y.Q.; Portone, A.; Rubinacci, G. 3D effects of conducting structures on RWMs control in ITER. In Proceedings of the 34th EPS Conference, Warsaw, Poland, 2–6 July 2007. [Google Scholar]
- Bossavit, A. Computational Electromagnetism: Variational Formulations, Complementarity, Edge Elements; Academic Press: Cambridge, MA, USA, 1998. [Google Scholar]
- Biro, O.; Preis, K. An edge finite element eddy current formulation using a reduced magnetic and a current vector potential. IEEE Trans. Magn. 2000, 36, 3128–3130. [Google Scholar] [CrossRef]
- Bossavit, A.; Verite, J. The “TRIFOU” Code: Solving the 3-D eddy-currents problem by using H as state variable. IEEE Trans. Magn. 1983, 19, 2465–2470. [Google Scholar] [CrossRef]
- Kurz, S.; Fetzer, J.; Lehner, G.; Rucker, W.M. A novel formulation for 3D eddy current problems with moving bodies using a Lagrangian description and BEM-FEM coupling. IEEE Trans. Magn. 1998, 34, 3068–3073. [Google Scholar] [CrossRef]
- Kameari, A. Transient eddy current analysis on thin conductors with arbitrary connections and shapes. J. Comput. Phys. 1981, 42, 124–140. [Google Scholar] [CrossRef]
- Bossavit, A. On the numerical analysis of eddy-current problems. Comput. Methods Appl. Mech. Eng. 1981, 27, 303–318. [Google Scholar] [CrossRef]
- Blum, J.; Leloup, C.; Dupas, L.; Thooris, B. Eddy current calculations for the Tore Supra Tokamak. IEEE Trans. Magn. 1983, 19, 2461–2464. [Google Scholar] [CrossRef]
- Albanese, R.; Rubinacci, G. Integral formulation for 3D eddy-current computation using edge elements. IEE Proc. A Phys. Sci. Meas. Instrum. Manag. Educ. Rev. 1988, 135, 457–462. [Google Scholar] [CrossRef]
- Albanese, R.; Rubinacci, G. Finite Element Methods for the Solution of 3D Eddy Current Problems. Adv. Imaging Electron. Phys. 1998, 102, 1–86. [Google Scholar] [CrossRef]
- Roccella, R.; Boccaccini, L.; Meyder, R.; Raff, S. Assessment of EM loads on the EU HCPB TBM during plasma disruption and normal operating scenario including the ferromagnetic effect. Fusion Eng. Des. 2008, 83, 1212–1216. [Google Scholar] [CrossRef]
- Testoni, P.; Albanese, R.; Lucca, F.; Roccella, M.; Portone, A.; Rubinacci, G.; Ventre, S.; Villone, F. Status of the EU domestic agency electromagnetic analyses of ITER vacuum vessel and blanket modules. Fusion Eng. Des. 2013, 88, 1934–1937. [Google Scholar] [CrossRef]
- Zhu, J.; Bykov, V.; Pagani, I.; Lucca, F.; Wegener, L.; Bosch, H.-S. Specific features of eddy current analysis with ANSYS® for fast plasma current decay event in W7-X. Fusion Eng. Des. 2021, 163, 112136. [Google Scholar] [CrossRef]
- Kwon, S.; Hong, S.-H. Estimation of electromagnetic loads including eddy and halo current on the major disruption for the K-DEMO divertor. Fusion Eng. Des. 2021, 168, 112592. [Google Scholar] [CrossRef]
- Albanese, R.; Rubinacci, G.; Tamburrino, A.; Ventre, S.; Villone, F. A fast 3D eddy current integral formulation. COMPEL Int. J. Comput. Math. Electr. Electron. Eng. 2001, 20, 317–331. [Google Scholar] [CrossRef]
- Phillips, J.R.; White, J.K. A precorrected-FFT method for electrostatic analysis of complicated 3-D structures. IEEE Trans. Comput.-Aided Des. Integr. Circuits Syst. 1997, 16, 1059–1072. [Google Scholar] [CrossRef] [Green Version]
- Tamburrino, A.; Ventre, S.; Rubinacci, G. An FFT integral formulation using edge-elements for Eddy Current Testing. Int. J. Appl. Electromagn. Mech. 2000, 11, 141–162. [Google Scholar] [CrossRef]
- Greengard, L.; Rokhlin, V. A fast algorithm for particle simulations. J. Comput. Phys. 1987, 73, 325–348. [Google Scholar] [CrossRef]
- Cheng, H.; Greengard, L.; Rokhlin, V. A Fast Adaptive Multipole Algorithm in Three Dimensions. J. Comput. Phys. 1999, 155, 468–498. [Google Scholar] [CrossRef]
- Rubinacci, G.; Ventre, S.; Villone, F.; Liu, Y. A fast technique applied to the analysis of Resistive Wall Modes with 3D conducting structures. J. Comput. Phys. 2008, 228, 1562–1572. [Google Scholar] [CrossRef]
- Maffucci, A.; Rubinacci, G.; Tamburrino, A.; Ventre, S.; Villone, F. Fast Low-Frequency Impedance Extraction using a Volumetric Three-Dimensional Integral Formulation. In Proceedings of the 23rd Annual Review of Progress in Applied Computational Electromagnetics (ACES 2007), Verona, Italy, 19–23 March 2007. [Google Scholar]
- Van Loan, C.F.; Golub, G.H. Matrix Computations, 3rd ed.; Johns Hopkins University Press: London, UK, 1996. [Google Scholar]
- EFDA, European Fusion Development. The ITER Project. Available online: www.iter.org (accessed on 15 September 2019).
- Rubinacci, G.; Tamburrino, A.; Villone, F. Circuits/fields coupling and multiply connected domains in integral formulations. IEEE Trans. Magn. 2002, 38, 581–584. [Google Scholar] [CrossRef]
- Börm, S.; Grasedyck, L.; Hackbusch, W. Hierarchical Matrices. Institut für Mathematik in den Naturwissenschaften. Available online: https://www.researchgate.net/publication/277293203_Hierarchical_Matrices (accessed on 24 March 2022).
- Ventre, S.; Cau, F.; Chiariello, A.; Giovinco, G.; Maffucci, A.; Villone, F. Fast and Accurate Solution of Integral Formulations of Large MQS Problems Based on Hybrid OpenMP–MPI Parallelization. Appl. Sci. 2022, 12, 627. [Google Scholar] [CrossRef]
- Fresa, R.; Rubinacci, G.; Ventre, S. An eddy current integral formulation on parallel computer systems. Int. J. Numer. Methods Eng. 2005, 62, 1127–1147. [Google Scholar] [CrossRef]
- Scalable Linear Algebra PACKage. Available online: www.scalapack.org (accessed on 12 November 2019).
- Gruber, O.; Lackner, K.; Pautasso, G.; Seidel, U.; Streibl, B. Vertical displacement events and halo currents. Plasma Phys. Control. Fusion 1993, 35, B191–B204. [Google Scholar] [CrossRef]
- Kapur, S.; Long, D.E. N-body problems: IES3: Efficient electrostatic and electromagnetic simulation. IEEE Comput. Sci. Eng. 1998, 5, 60–67. [Google Scholar] [CrossRef]
- Voltolina, D.; Torchio, R.; Bettini, P.; Specogna, R.; Alotto, P. Optimized cycle basis in volume integral formulations for large scale eddy-current problems. Comput. Phys. Commun. 2021, 265, 108004. [Google Scholar] [CrossRef]
- Rubinacci, G.; Tamburrino, A. Automatic Treatment of Multiply Connected Regions in Integral Formulations. IEEE Trans. Magn. 2010, 46, 2791–2794. [Google Scholar] [CrossRef]
- Burkholder, R.J.; Lee, J.-F. Fast Dual-MGS Block-Factorization Algorithm for Dense MoM Matrices. IEEE Trans. Antennas Propag. 2004, 52, 1693–1699. [Google Scholar] [CrossRef]
- The MPI Forum. MPI: A Message Passing Interface. In Proceedings of the 1993 ACM/IEEE Conference on Supercomputing. Supercomputing ‘93, Portland, OR, USA, 15–19 November 1993.
- Barney, B. Introduction to Parallel Computing; Livermore National Laboratory: Livermore, CA, USA, 2015.
- Saad, Y.; Schultz, M.H. GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems. SIAM J. Sci. Stat. Comput. 1986, 7, 856–869. [Google Scholar] [CrossRef] [Green Version]
- The HLS Mathematical Software Library, Symmetric Sparse Matrix: Compute elimination orderings. Available online: https://www.hsl.rl.ac.uk/specs/hsl_mc68.pdf (accessed on 24 March 2022).
- Rubinacci, G.; Tamburrino, A.; Ventre, S.; Villone, F. Application of n-Fold Rotational Symmetries to Eddy Currents Integral Model in the Time Domain. IEEE Trans. Magn. 2019, 56, 1–4. [Google Scholar] [CrossRef]
- Cau, F.; Bessette, D.; Bauer, P.; Gauthier, F.; Portone, A.; Testoni, P.; Ventre, S. Update of Joule Losses Calculation in the ITER Cold Structures during Fast Plasma Transients. IEEE Trans. Appl. Supercond. 2020, 30, 1–4. [Google Scholar] [CrossRef]
- Cau, F.B.; Bessette, D.; D’Amico, G.; Portone, A.; Rubinacci, G.; Testoni, P.; Ventre, S.; Villone, F. Joule Losses in the ITER Cold Structures during Plasma Transients. IEEE Trans. Appl. Supercond. 2016, 26, 1–5. [Google Scholar] [CrossRef]
- Foussat, A.; Libeyre, P.; Mitchell, N.; Gribov, Y.; Jong, C.T.J.; Bessette, D.; Gallix, R.; Bauer, P.; Sahu, A. Overview of the ITER Correction Coils Design. IEEE Trans. Appl. Supercond. 2010, 20, 402–406. [Google Scholar] [CrossRef]
- Bachmann, C.; Sugihara, M.; Roccella, R.; Sannazzaro, G.; Gribov, Y.; Riccardo, V.; Hender, T.; Gerasimov, S.; Pautasso, G.; Belov, A.V.; et al. Specification of asymmetric VDE loads of the ITER tokamak. Fusion Eng. Des. 2011, 86, 1915–1919. [Google Scholar] [CrossRef]
- Börm, S.; Grasedyck, L.; Hackbusch, W. Introduction to hierarchical matrices with applications. Eng. Anal. Bound. Elem. 2003, 27, 405–422. [Google Scholar] [CrossRef] [Green Version]
- Hackbusch, W. A sparse matrix arithmetic based on H-matrices. Part I: Introduction to H-matrices. Computing 1999, 62, 89–108. [Google Scholar] [CrossRef]
- Hackbusch, W.; Khoromskij, B.N. A Sparse -Matrix Arithmetic, Part II: Application to Multidimensional Problems. Computing 2000, 64, 21–47. [Google Scholar] [CrossRef]
- H-Lib Pro Documentation. Available online: https://www.hlibpro.com/ (accessed on 21 March 2022).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Cau, F.; Chiariello, A.G.; Rubinacci, G.; Scalera, V.; Tamburrino, A.; Ventre, S.; Villone, F. A Fast Matrix Compression Method for Large Scale Numerical Modelling of Rotationally Symmetric 3D Passive Structures in Fusion Devices. Energies 2022, 15, 3214. https://doi.org/10.3390/en15093214
Cau F, Chiariello AG, Rubinacci G, Scalera V, Tamburrino A, Ventre S, Villone F. A Fast Matrix Compression Method for Large Scale Numerical Modelling of Rotationally Symmetric 3D Passive Structures in Fusion Devices. Energies. 2022; 15(9):3214. https://doi.org/10.3390/en15093214
Chicago/Turabian StyleCau, Francesca, Andrea Gaetano Chiariello, Guglielmo Rubinacci, Valentino Scalera, Antonello Tamburrino, Salvatore Ventre, and Fabio Villone. 2022. "A Fast Matrix Compression Method for Large Scale Numerical Modelling of Rotationally Symmetric 3D Passive Structures in Fusion Devices" Energies 15, no. 9: 3214. https://doi.org/10.3390/en15093214
APA StyleCau, F., Chiariello, A. G., Rubinacci, G., Scalera, V., Tamburrino, A., Ventre, S., & Villone, F. (2022). A Fast Matrix Compression Method for Large Scale Numerical Modelling of Rotationally Symmetric 3D Passive Structures in Fusion Devices. Energies, 15(9), 3214. https://doi.org/10.3390/en15093214