A Parallel Direct Method for Finite Element Electromagnetic Computations Based on Domain Decomposition

High performance parallel computing and direct (factorization-based) solution methods have been the two main trends in electromagnetic computations in recent years. When time-harmonic (frequency-domain) Maxwell\u27s equation are directly discretized with the Finite Element Method (FEM) or other Partial Differential Equation (PDE) methods, the resulting linear system of equations is sparse and indefinite, thus harder to efficiently factorize serially or in parallel than alternative methods e.g. integral equation solutions, that result in dense linear systems. State-of-the-art sparse matrix direct solvers such as MUMPS and PARDISO don\u27t scale favorably, have low parallel efficiency and high memory footprint. This work introduces a new class of sparse direct solvers based on domain decomposition method, termed Direct Domain Decomposition Method (D3M), which is reliable, memory efficient, and offers very good parallel scalability for arbitrary 3D FEM problems. Unlike recent trends in approximate/low-rank solvers, this method focuses on `numerically exact\u27 solution methods as they are more reliable for complex `real-life\u27 models. The proposed method leverages physical insights at every stage of the development through a new symmetric domain decomposition method (DDM) with one set of Lagrange multipliers. Applying a special regularization scheme at the interfaces, either artificial loss or gain is introduced to each domain to eliminate non-physical internal resonances. A block-wise recursive algorithm based on Takahashi relationship is proposed for the efficient computation of discrete Dirichlet-to-Neumann (DtN) map to reduce the volumetric problem from all domains into an auxiliary surfacial problem defined on the domain interfaces only. Numerical results show up to 50% run-time saving in DtN map computation using the proposed block-wise recursive algorithm compared to alternative approaches. The auxiliary unknowns on the domain interfaces form a considerably (approximately an order of magnitude) smaller block-wise sparse matrix, which is efficiently factorized using a customized block LDL$^T$ factorization with restricted pivoting to ensure stability. The parallelization of the proposed D3M is realized based on Directed Acyclic Graph (DAG). Recent advances in parallel dense direct solvers, have shifted toward parallel implementation that rely on DAG scheduling to achieve highly efficient asynchronous parallel execution. However, adaptation of such schemes to sparse matrices is harder and often impractical. In D3M, computation of each domain\u27s discrete DtN map ``embarrassingly parallel\u27\u27, whereas the customized block LDLT is suitable for a block directed acyclic graph (B-DAG) task scheduling, similar to that used in dense matrix parallel direct solvers. In this approach, computations are represented as a sequence of small tasks that operate on domains of DDM or dense matrix blocks of the reduced matrix. These tasks can be statically scheduled for parallel execution using their DAG dependencies and weights that depend on estimates of computation and communication costs. Comparisons with state-of-the-art exact direct solvers on electrically large problems suggest up to 20% better parallel efficiency, 30% - 3X less memory and slightly faster in runtime, while maintaining the same accuracy

Correlations in Nuclear Matter

We analyze the nuclear matter correlation properties in terms of the pair correlation function. To this aim we systematically compare the results for the variational method in the Lowest Order Constrained Variational (LOCV) approximation and for the Bruekner-Hartree-Fock (BHF) scheme. A formal link between the Jastrow correlation factor of LOCV and the Defect Function (DF) of BHF is established and it is shown under which conditions and approximations the two approaches are equivalent. From the numerical comparison it turns out that the two correlation functions are quite close, which indicates in particular that the DF is approximately local and momentum independent. The Equations of State (EOS) of Nuclear Matter in the two approaches are also compared. It is found that once the three-body forces (TBF) are introduced the two EOS are fairly close, while the agreement between the correlation functions holds with or without TBF.Comment: 11 figure

In This Edition

Over the past two decades, considerable attention has been given to the subject of climate mitigation, of the development of laws and policies to reduce the amount of greenhouse gas emissions that are contributing to global warming. More recently, in addition to climate mitigation, attention has turned to the question of climate adaptation, of the development of law and policies that respond to the environmental consequences of climate change. In this Pacific Region Edition of the Golden Gate University Environmental Law Journal, titled Climate Resiliency – California Prepares for an Altered Environment, we develop this theme of climate adaptation. Our edition features six articles, three from professionals in the legal field and three from students at Golden Gate University

In This Edition

Over the past two decades, considerable attention has been given to the subject of climate mitigation, of the development of laws and policies to reduce the amount of greenhouse gas emissions that are contributing to global warming. More recently, in addition to climate mitigation, attention has turned to the question of climate adaptation, of the development of law and policies that respond to the environmental consequences of climate change. In this Pacific Region Edition of the Golden Gate University Environmental Law Journal, titled Climate Resiliency – California Prepares for an Altered Environment, we develop this theme of climate adaptation. Our edition features six articles, three from professionals in the legal field and three from students at Golden Gate University

Hybrid stars within the framework of the Sigma-Omega-Rho model combined with the MIT and NJL models

In this paper, we investigate the structure of hybrid stars consisting of hadrons (neutrons, protons, sigmas, lambdas), leptons (electrons, muons), and quarks (up, down, strange). We use a relativistic mean-field (RMF) model namely the Sigma-omega-rho model for the hadronic phase and the MIT bag model as well as the NJL model for the quark phase. In addition, Maxwell and Gibbs conditions are employed to investigate the hadron-Quark phase transition. Finally, by obtaining the mass-radius relation, $M (M_{sun}) \leqslant 2.07$ is predicted for such hybrid stars.Comment: 23 pages, 10 figure