hyper.deal: An Efficient, Matrix-free Finite-element Library for High-dimensional Partial Differential Equations

This work presents the efficient, matrix-free finite-element library hyper.deal for solving partial differential equations in two up to six dimensions with high-order discontinuous Galerkin methods. It builds upon the low-dimensional finite-element library deal.II to create complex low-dimensional meshes and to operate on them individually. These meshes are combined via a tensor product on the fly, and the library provides new special-purpose highly optimized matrix-free functions exploiting domain decomposition as well as shared memory via MPI-3.0 features. Both node-level performance analyses and strong/weak-scaling studies on up to 147,456 CPU cores confirm the efficiency of the implementation. Results obtained with the library hyper.deal are reported for high-dimensional advection problems and for the solution of the Vlasov–Poisson equation in up to six-dimensional phase space.

Complexity of COVID-19 Dynamics

With population explosion and globalization, the spread of infectious diseases has been a major concern. In 2019, a newly identified type of Coronavirus caused an outbreak of respiratory illness, popularly known as COVID-19, and became a pandemic. Although enormous efforts have been made to understand the spread of COVID-19, our knowledge of the COVID-19 dynamics still remains limited. The present study employs the concepts of chaos theory to examine the temporal dynamic complexity of COVID-19 around the world. The false nearest neighbor (FNN) method is applied to determine the dimensionality and, hence, the complexity of the COVID-19 dynamics. The methodology involves: (1) reconstruction of a single-variable COVID-19 time series in a multi-dimensional phase space to represent the underlying dynamics; and (2) identification of “false” neighbors in the reconstructed phase space and estimation of the dimension of the COVID-19 series. For implementation, COVID-19 data from 40 countries/regions around the world are studied. Two types of COVID-19 data are analyzed: (1) daily COVID-19 cases; and (2) daily COVID-19 deaths. The results for the 40 countries/regions indicate that: (1) the dynamics of COVID-19 cases exhibit low- to medium-level complexity, with dimensionality in the range 3 to 7; and (2) the dynamics of COVID-19 deaths exhibit complexity anywhere from low to high, with dimensionality ranging from 3 to 13. The results also suggest that the complexity of the dynamics of COVID-19 deaths is greater than or at least equal to that of the dynamics of COVID-19 cases for most (three-fourths) of the countries/regions. These results have important implications for modeling and predicting the spread of COVID-19 (and other infectious diseases), especially in the identification of the appropriate complexity of models.

Spray dispersion regimes following atomization in a turbulent co-axial gas jet

A canonical co-axial round-jet two-fluid atomizer where atomization occurs over a wide range of momentum ratios: $M=1.9 - 376.4$ is studied. The near field of the spray, where the droplet formation process takes place, is characterized and linked to droplet dispersion in the far field of the jet. Counterintuitively, our results indicate that in the low-momentum regime, increasing the momentum in the gas phase leads to less droplet dispersion. A critical momentum ratio of the order of $M_c=50$ , that separates this regime from a high-momentum one with less dispersion, is found in both the near and far fields. A phenomenological model is proposed that determines the susceptibility of droplets to disperse beyond the nominal extent of the gas phase based on a critical Stokes number, $St=\tau _p/T_E=1.9$ , formulated based on the local Eulerian large scale eddy turnover time, $T_E$ , and the droplets’ response time, $\tau _p$ . A two-dimensional phase space summarizes the extent of these different regimes in the context of spray characteristics found in the literature.

Constant-roll inflation in the generalized SU(2) Proca theory

The generalized SU(2) Proca theory (GSU2P for short) is a variant of the well known generalized Proca theory (GP for short) where the vector field belongs to the Lie algebra of the SU(2) group of global transformations under which the action is made invariant. New interesting possibilities arise in this framework because of the existence of new interactions of purely non-Abelian character and new configurations of the vector field that result in spatial spherical symmetry and the cosmological dynamics being driven by the propagating degrees of freedom. We study the two-dimensional phase space of the system that results when the cosmic triad configuration is employed in the Friedmann-Lemaitre-Robertson-Walker background and find an attractor curve whose attraction basin covers almost all the allowed region. Such an attractor curve corresponds to a primordial inflationary solution that has the following characteristic properties: 1.) it is a de Sitter solution whose Hubble parameter is regulated by a generalized version of the SU(2) group coupling constant, 2.) it is constant-roll including, as opposite limiting cases, the slow-roll and ultra slow-roll varieties, 3.) a number of e-folds $N > 60$ is easily reached, 4.) it has a graceful exit into a radiation dominated period powered by the canonical kinetic term of the vector field and the Einstein-Hilbert term. The free parameters of the action are chosen such that the tensor sector of the theory is the same as that of general relativity at least up to second-order perturbations, thereby avoiding the presence of ghost and Laplacian instabilities in the tensor sector as well as making the gravity waves propagate at light speed. This is a proof of concept of the interesting properties we could find in this scenario when the coupling constants be replaced by general coupling functions.

Position-dependent mass quantum systems and ADM formalism

The classical Einstein-Hilbert (EH) action for general relativity (GR) is shown to be formally analogous to the classical system with position-dependent mass (PDM) models. The analogy is developed and used to build the covariant classical Hamiltonian as well as defining an alternative phase portrait for GR. The set of associated Hamilton’s equations in the phase space is presented as a first-order system dual to the Einstein field equations. Following the principles of quantum mechanics, I build a canonical theory for the classical general. A fully consistent quantum Hamiltonian for GR is constructed based on adopting a high dimensional phase space. It is observed that the functional wave equation is timeless. As a direct application, I present an alternative wave equation for quantum cosmology. In comparison to the standard Arnowitt-Deser-Misner(ADM) decomposition and quantum gravity proposals, I extended my analysis beyond the covariant regime when the metric is decomposed into the 3+13+1 dimensional ADM decomposition. I showed that an equal dimensional phase space can be obtained if one applies ADM decomposed metric.

Confinement in 4D: An Attempt at Classical Understanding

In this review, we revisit our approach to constructing an effective theory for Abelian and Non-Abelian gauge theories in 4D. Our goal is to have an effective theory that provides a simple classical picture of the main qualitatively important features of these theories. We set out to ensure the presence of the massless photons—Goldstone bosons in Abelian theory and their disappearance in the Non-Abelian case—accompanied by the formation of confining strings between charged states. Our formulation avoids using vector fields and instead operates with the basic degrees of freedom that are the scalar fields of a nonlinear σ-model. The Mark 1 model we study turns out to have a large global symmetry group-the 2D diffeomorphism invariance in the Abelian limit, which is isomorphic to the group of all canonical transformations in the classical two dimensional phase space. This symmetry is not present in QED, and we eliminate it by “gauging” this infinite dimensional global group. Introducing additional modifications to the model (Mark 2), we are able to prove that the “Abelian” version is equivalent to the theory of a free photon. Achieving the desired property in the “Non-Abelian” regime turns out to be tricky. We are able to introduce a perturbation that leads to the formation of confining strings in our Mark 1 model. These strings have somewhat unusual properties, in that their profile does not decay exponentially away from the center of the string. In addition, the perturbation explicitly breaks the diffeomorphism invariance. Preserving this invariance in the gauged model as well as achieving confining strings in Mark 2 model remains an open question.

CONTROL SYSTEM DEPENDING ON A PARAMETER

A nonlinear control system depending on a parameter is considered in a finite-dimensional Euclidean space and on a finite time interval. The dependence on the parameter of the reachable sets and integral funnels of the corresponding differential inclusion system is studied. Under certain conditions on the control system, the degree of this dependence on the parameter is estimated. Problems of targeting integral funnels to a target set in the presence of an obstacle in strict and soft settings are considered. An algorithm for the numerical solution of this problem in the soft setting has been developed. An estimate of the error of the developed algorithm is obtained. An example of solving a specific problem for a control system in a two-dimensional phase space is given.

Geometric particle-in-cell methods for the Vlasov–Maxwell equations with spin effects

We propose a numerical scheme to solve the semiclassical Vlasov–Maxwell equations for electrons with spin. The electron gas is described by a distribution function $f(t,{\boldsymbol x},{{{\boldsymbol p}}}, {\boldsymbol s})$ that evolves in an extended 9-dimensional phase space $({\boldsymbol x},{{{\boldsymbol p}}}, {\boldsymbol s})$ , where $\boldsymbol s$ represents the spin vector. Using suitable approximations and symmetries, the extended phase space can be reduced to five dimensions: $(x,{{p_x}}, {\boldsymbol s})$ . It can be shown that the spin Vlasov–Maxwell equations enjoy a Hamiltonian structure that motivates the use of the recently developed geometric particle-in-cell (PIC) methods. Here, the geometric PIC approach is generalized to the case of electrons with spin. Total energy conservation is very well satisfied, with a relative error below $0.05\,\%$ . As a relevant example, we study the stimulated Raman scattering of an electromagnetic wave interacting with an underdense plasma, where the electrons are partially or fully spin polarized. It is shown that the Raman instability is very effective in destroying the electron polarization.