Space-time Formulation Research Articles

Existence results for the A – φ – B magnetodynamic formulation of the Maxwell system with skin and proximity effects

ABSTRACT The magnetodynamic Maxwell system given in its potential and space-time formulation is a popular model considered in the engineering community. It allows to model some phenomena such as eddy current losses in multiple turn winding. Indeed, in some cases, they can significantly alter the performance of the devices, and consequently can no more be neglegted. It turns out that this model is not yet analytically studied, therefore we here consider its well-posedness. First, the existence of strong solutions with the help of the theory of Showalter on degenerated parabolic problems is established. Second, using energy estimates, existence and uniqueness of the weak solution of the is deduced.

A Novel Staggered Semi-implicit Space-Time Discontinuous Galerkin Method for the Incompressible Navier-Stokes Equations

A new high order accurate staggered semi-implicit space-time discontinuous Galerkin (DG) method is presented for the simulation of viscous incompressible flows on unstructured triangular grids in two space dimensions. The staggered DG scheme defines the discrete pressure on the primal triangular mesh, while the discrete velocity is defined on a staggered edge-based dual quadrilateral mesh. In this paper, a new pair of equal-order-interpolation velocity-pressure finite elements is proposed. On the primary triangular mesh (the pressure elements) the basis functions are piecewise polynomials of degree $N$ and are allowed to jump on the boundaries of each triangle. On the dual mesh instead (the velocity elements), the basis functions consist in the union of piecewise polynomials of degree $N$ on the two subtriangles that compose each quadrilateral and are allowed to jump only on the dual element boundaries, while they are continuous inside. In other words, the basis functions on the dual mesh are built by continuous finite elements on the subtriangles. This choice allows the construction of an efficient, quadrature-free and memory saving algorithm. In our coupled space-time pressure correction formulation for the incompressible Navier-Stokes equations, arbitrary high order of accuracy in time is achieved through the use of time-dependent test and basis functions, in combination with simple and efficient Picard iterations. Several numerical tests on classical benchmarks confirm that the proposed method outperforms existing staggered semi-implicit space-time DG schemes, not only from a computer memory point of view, but also concerning the computational time.

Massless Infinite Spin Representations

A new twistorial formulation of a massless particle with infinite spin is presented. The twistor field of infinite spin is found; the twistor field transformation for obtaining the space-time field formulation is found. The $$N = 1$$ infinite spin supermultiplet formed by the fields of the infinite integer and half-integer spins is presented.

A space-time discontinuous Galerkin method for the elastic wave equation

In this work we present a new high order space-time discretization method based on a discontinuous Galerkin paradigm for the second order visco-elastodynamics equation. After introducing the method, we show that the resulting space-time discontinuous Galerkin formulation is well-posed, stable and retains optimal rate of convergence with respect to the discretization parameters, namely the mesh size and the polynomial approximation degree. A set of two and three-dimensional numerical experiments confirms the theoretical bounds.

Newton-Cartan submanifolds and fluid membranes.

We develop the geometric description of submanifolds in Newton-Cartan spacetime. This provides the necessary starting point for a covariant spacetime formulation of Galilean-invariant hydrodynamics on curved surfaces. We argue that this is the natural geometrical framework to study fluid membranes in thermal equilibrium and their dynamics out of equilibrium. A simple model of fluid membranes that only depends on the surface tension is presented and, extracting the resulting stresses, we show that perturbations away from equilibrium yield the standard result for the dispersion of elastic waves. We also find a generalization of the Canham-Helfrich bending energy for lipid vesicles that takes into account the requirements of thermal equilibrium.

POD–Galerkin Model Order Reduction for Parametrized Time Dependent Linear Quadratic Optimal Control Problems in Saddle Point Formulation

In this work we deal with parametrized time dependent optimal control problems governed by partial differential equations. We aim at extending the standard saddle point framework of steady constraints to time dependent cases. We provide an analysis of the well-posedness of this formulation both for parametrized scalar parabolic constraint and Stokes governing equations and we propose reduced order methods as an effective strategy to solve them. Indeed, on one hand, parametrized time dependent optimal control is a very powerful mathematical model which is able to describe several physical phenomena, on the other, it requires a huge computational effort. Reduced order methods are a suitable approach to have rapid and accurate simulations. We rely on POD–Galerkin reduction over the physical and geometrical parameters of the optimality system in a space-time formulation. Our theoretical results and our methodology are tested on two examples: a boundary time dependent optimal control for a Graetz flow and a distributed optimal control governed by time dependent Stokes equations. With these two test cases the convenience of the reduced order modelling is further extended to the field of time dependent optimal control.

High Order ADER Schemes for Continuum Mechanics

In this paper we first review the development of high order ADER finite volume and ADER discontinuous Galerkin schemes on fixed and moving meshes, since their introduction in 1999 by Toro et al. We show the modern variant of ADER based on a space-time predictor-corrector formulation in the context of ADER discontinuous Galerkin schemes with a posteriori subcell finite volume limiter on fixed and moving grids, as well as on space-time adaptive Cartesian AMR meshes. We then present and discuss the unified symmetric hyperbolic and thermodynamically compatible (SHTC) formulation of continuum mechanics developed by Godunov, Peshkov and Romenski (GPR model), which allows to describe fluid and solid mechanics in one single and unified first order hyperbolic system. In order to deal with free surface and moving boundary problems, a simple diffuse interface approach is employed, which is compatible with Eulerian schemes on fixed grids as well as direct Arbitrary-Lagrangian-Eulerian methods on moving meshes. We show some examples of moving boundary problems in fluid and solid mechanics.

Quasistatic analysis of elastoplastic structures by the proper generalized decomposition in a space-time approach

In this paper, we present the space-time formulation using the Proper Generalized Decomposition (PGD) to solve an elastoplastic structure under quasistatic loading. Here, the solution, dependent on space and time, is built by iteratively summing up enrichments until convergence is obtained. Thereby, each enrichment functional is represented by the dyadic product of spatial and temporal functions, which we call spatial and temporal modes. The proposed strategy is included into a finite element framework resulting in two coupled equations, one in space and one in time. Those two equations are solved simultaneously using the fixed-point algorithm. We visualize and mechanically interpret the evolution of the spatial and temporal modes within the proposed iterative scheme. The test structure is represented by the discretization of an elastoplastic frame geometry subjected to an external cyclic force history. The results are compared with the conventional iteration method. It is shown that inelastic deformation can be directly located by the evolution of the temporal functions. However, the corresponding spatial modes show a high correlation containing global elastic and plastic deformation patterns that can be interpreted by an orthonormal basis representation using the singular value decomposition.

A Unified Framework for the Solution of Hyperbolic PDE Systems Using High Order Direct Arbitrary-Lagrangian–Eulerian Schemes on Moving Unstructured Meshes with Topology Change

In this work, we review the family of direct Arbitrary-Lagrangian–Eulerian (ALE) finite vlume (FV) and discontinuous Galerkin (DG) schemes on moving meshes that at each time step are rearranged by explicitly allowing topology changes, in order to guarantee a robust mesh evolution even for high shear flow and very long evolution times. Two different techniques are presented: a local nonconforming approach for dealing with sliding lines, and a global regeneration of Voronoi tessellations for treating general unpredicted movements. Corresponding elements at consecutive times are connected in space-time to construct closed space-time control volumes, whose bottom and top faces may be polygons with a different number of nodes, with different neighbors, and even degenerate space-time sliver elements. Our final ALE FV-DG scheme is obtained by integrating, over these arbitrary shaped space-time control volumes, the space-time conservation formulation of the governing hyperbolic PDE system: so, we directly evolve the solution in time avoiding any remapping stage, being conservative and satisfying the GCL by construction. Arbitrary high order of accuracy in space and time is achieved through a fully discrete one-step predictor–corrector ADER approach, also integrated with well balancing techniques to further improve the accuracy and to maintain exactly even at discrete level many physical invariants of the studied system. A large set of different numerical tests has been carried out in order to check the accuracy and the robustness of our methods for both smooth and discontinuous problems, in particular in the case of vortical flows.

A quantum cellular automaton for one-dimensional QED

We propose a discrete spacetime formulation of quantum electrodynamics in one dimension (a.k.a the Schwinger model) in terms of quantum cellular automata, i.e. translationally invariant circuits of local quantum gates. These have exact gauge covariance and a maximum speed of information propagation. In this picture, the interacting quantum field theory is defined as a “convergent” sequence of quantum cellular automata, parameterized by the spacetime lattice spacing—encompassing the notions of continuum limit and renormalization, and at the same time providing a quantum simulation algorithm for the dynamics.

Space-time topology optimization for additive manufacturing

The design of optimal structures and the planning of (additive manufacturing) fabrication sequences have been considered typically as two separate tasks that are performed consecutively. In the light of recent advances in robot-assisted (wire-arc) additive manufacturing which enable addition of material along curved surfaces, we present a novel topology optimization formulation which concurrently optimizes the structure and the fabrication sequence. For this, two sets of design variables, i.e., a density field for defining the structural layout, and a time field which determines the fabrication process order, are simultaneously optimized. These two fields allow to generate a sequence of intermediate structures, upon which manufacturing constraints (e.g., fabrication continuity and speed) are imposed. The proposed space-time formulation is general, and is demonstrated on three fabrication settings, considering self-weight of the intermediate structures, process-dependent critical loads, and time-dependent material properties.

Bayesian modeling of spatiotemporal patterns of TB-HIV co-infection risk in Kenya

BackgroundTuberculosis (TB) and Human Immunodeficiency Virus (HIV) diseases are globally acknowledged as a public health challenge that exhibits adverse bidirectional relations due to the co-epidemic overlap. To understand the co-infection burden we used the case notification data to generate spatiotemporal maps that described the distribution and exposure hypotheses for further epidemiologic investigations in areas with unusual case notification levels.MethodsWe analyzed the TB and TB-HIV case notification data from the Kenya national TB control program aggregated for forty-seven counties over a seven-year period (2012–2018). Using spatiotemporal poisson regression models within the Integrated Nested Laplace Approach (INLA) paradygm, we modeled the risk of TB-HIV co-infection. Six competing models with varying space-time formulations were compared to determine the best fit model. We then assessed the geographic patterns and temporal trends of coinfection risk by mapping the posterior marginal from the best fit model.ResultsOf the total 608,312 TB case notifications, 194,129 were HIV co-infected. The proportion of TB-HIV co-infection was higher in females (39.7%) than in males (27.0%). A significant share of the co-infection was among adults aged 35 to 44 years (46.7%) and 45 to 54 years (42.1%). Based on the Bayesian Defiance Information (DIC) and the effective number of parameters (pD) comparisons, the spatiotemporal model allowing space-time interaction was the best in explaining the geographical variations in TB-HIV coinfection. The model results suggested that the risk of TB-HIV coinfection was influenced by infrastructure index (Relative risk (RR) = 5.75, Credible Interval (Cr.I) = (1.65, 19.89)) and gender ratio (RR = 5.81e−04, Cr. I = (1.06e−04, 3.18e−03). The lowest and highest temporal relative risks were in the years 2016 at 0.9 and 2012 at 1.07 respectively. The spatial pattern presented an increased co-infection risk in a number of counties. For the spatiotemporal interaction, only a few counties had a relative risk greater than 1 that varied in different years.ConclusionsWe identified elevated risk areas for TB/HIV co-infection and fluctuating temporal trends which could be because of improved TB case detection or surveillance bias caused by spatial heterogeneity in the co-infection dynamics. Focused interventions and continuous TB-HIV surveillance will ensure adequate resource allocation and significant reduction of HIV burden amongst TB patients.

Petrov–Galerkin space-time hp-approximation of parabolic equations in H1/2

Abstract We analyse a class of variational space-time discretizations for a broad class of initial boundary value problems for linear, parabolic evolution equations. The space-time variational formulation is based on fractional Sobolev spaces of order $1/2$ and the Riemann–Liouville derivative of order $1/2$ with respect to the temporal variable. It accommodates general, conforming space discretizations and naturally accommodates discretization of infinite horizon evolution problems. We prove an inf-sup condition for $hp$-time semidiscretizations with an explicit expression of stable test functions given in terms of Hilbert transforms of the corresponding trial functions; inf-sup constants are independent of temporal order and the time-step sequences, allowing quasi-optimal, high-order discretizations on graded time-step sequences, and also $hp$-time discretizations. For solutions exhibiting Gevrey regularity in time and taking values in certain weighted Bochner spaces, we establish novel exponential convergence estimates in terms of $N_t$, the number of (elliptic) spatial problems to be solved. The space-time variational setting allows general space discretizations and, in particular, for spatial $hp$-FEM discretizations. We report numerical tests of the method for model problems in one space dimension with typical singular solutions in the spatial and temporal variable. $hp$-discretizations in both spatial and temporal variables are used without any loss of stability, resulting in overall exponential convergence of the space-time discretization.

Space–time variational methods for Maxwell's equations

AbstractThe efficient and accurate numerical solution of the time–dependent Maxwell equations is one of the most challenging tasks, see, e.g., [1]. Besides semi–discretization methods such as the method of lines and Laplace transformation approaches, space–time variational formulations became popular in recent years. Here the variational principle is applied simultaneously in space and time, which later requires the solution of the global linear system of algebraic equations. But this can be done in parallel, and the space–time formulation allows for an adaptive resolution of the solution in space and time simultaneously. Following previous work [3,5] on the acoustic wave equation we present two variational formulations for the solution of the electromagnetic wave equation.

Noncommutative gauge theories on D-branes in non-geometric backgrounds

We investigate the noncommutative gauge theories arising on the worldvolumes of D-branes in non-geometric backgrounds obtained by T-duality from twisted tori. We revisit the low-energy effective description of D-branes on three-dimensional T-folds, examining both cases of parabolic and elliptic twists in detail. We give a detailed description of the decoupling limits and explore various physical consequences of the open string non-geometry. The T-duality monodromies of the non-geometric backgrounds lead to Morita duality monodromies of the noncommutative Yang-Mills theories induced on the D-branes. While the parabolic twists recover the well-known examples of noncommutative principal torus bundles from topological T-duality, the elliptic twists give new examples of noncommutative fibrations with non-geometric torus fibres. We extend these considerations to D-branes in backgrounds with R-flux, using the doubled geometry formulation, finding that both the non-geometric background and the D-brane gauge theory necessarily have explicit dependence on the dual coordinates, and so have no conventional formulation in spacetime.

High order direct Arbitrary-Lagrangian-Eulerian (ALE) PNPM schemes with WENO Adaptive-Order reconstruction on unstructured meshes

In this work we present a conservative WENO Adaptive Order (AO) reconstruction operator applied to an explicit one-step Arbitrary-Lagrangian-Eulerian (ALE) discontinuous Galerkin (DG) method. The spatial order of accuracy is improved by reconstructing higher order piecewise polynomials of degree M>N, starting from the underlying polynomial solution of degree N provided by the DG scheme. High order of accuracy in time is achieved by the ADER approach, making use of an element-local space-time Galerkin finite element predictor that arises from a one-step time integration procedure. As a result, space-time polynomials of order M+1 are obtained and used to perform the time evolution of the numerical solution adopting a fully explicit DG scheme.To maintain algorithm simplicity, the mesh motion is restricted to be carried out using straight lines, hence the old mesh configuration at time tn is connected with the new one at time tn+1 via space-time segments, which result in space-time control volumes on which the governing equations have to be integrated in order to obtain the time evolution of the discrete solution. Our algorithm falls into the category of direct Arbitrary-Lagrangian-Eulerian (ALE) schemes, where the governing PDE system is directly discretized relying on a space-time conservation formulation and which already takes into account the new grid geometry directly during the computation of the numerical fluxes. A local rezoning strategy might be used in order to locally optimize the mesh quality and avoiding the generation of invalid elements with negative determinant. The proposed approach reduces to direct ALE finite volume schemes if N=0, while explicit direct ALE DG schemes are recovered in the case of N=M.In order to stabilize the DG solution, an a priori WENO based limiting technique is employed, that makes use of the numerical solution inside the element under consideration and its neighbor cells to find a less oscillatory polynomial approximation. By using a modal basis in a reference element, the evaluation of the oscillation indicators is very easily and efficiently carried out, hence allowing higher order modes to be properly limited, while leaving the zero-th order mode untouched for ensuring conservation.Numerical convergence rates for 2≤N,M≤4 are presented as well as a wide set of benchmark test problems for hydrodynamics on moving and fixed unstructured meshes.

Variational formulations for explicit Runge-Kutta Methods

Variational space-time formulations for partial differential equations have been of great interest in the last decades, among other things, because they allow to develop mesh-adaptive algorithms. Since it is known that implicit time marching schemes have variational structure, they are often employed for adaptivity. Previously, Galerkin formulations of explicit methods were introduced for ordinary differential equations employing specific inexact quadrature rules. In this work, we prove that the explicit Runge-Kutta methods can be expressed as discontinuous-in-time Petrov-Galerkin methods for the linear diffusion equation. We systematically build trial and test functions that, after exact integration in time, lead to one, two, and general stage explicit Runge-Kutta methods. This approach enables us to reproduce the existing time-domain (goal-oriented) adaptive algorithms using explicit methods in time.

Thermal Quantum Spacetime

The intersection of thermodynamics, quantum theory and gravity has revealed many profound insights, all the while posing new puzzles. In this article, we discuss an extension of equilibrium statistical mechanics and thermodynamics potentially compatible with a key feature of general relativity, background independence; and we subsequently use it in a candidate quantum gravity system, thus providing a preliminary formulation of a thermal quantum spacetime. Specifically, we emphasise an information-theoretic characterisation of generalised Gibbs equilibrium that is shown to be particularly suited to background independent settings, and in which the status of entropy is elevated to being more fundamental than energy. We also shed light on its intimate connections with the thermal time hypothesis. Based on this, we outline a framework for statistical mechanics of quantum gravity degrees of freedom of combinatorial and algebraic type, and apply it in several examples. In particular, we provide a quantum statistical basis for the origin of covariant group field theories, shown to arise as effective statistical field theories of the underlying quanta of space in a certain class of generalised Gibbs states.

Adaptive asynchronous time-stepping, stopping criteria, and a posteriori error estimates for fixed-stress iterative schemes for coupled poromechanics problems

In this paper we develop adaptive iterative coupling schemes for the Biot system modeling coupled poromechanics problems. We particularly consider the space–time formulation of the fixed-stress iterative scheme, in which we first solve the problem of flow over the whole space–time interval, then exploiting the space–time information for solving the mechanics. Two common discretizations of this algorithm are then introduced based on two coupled mixed finite element methods in-space and the backward Euler scheme in-time. Therefrom, adaptive fixed-stress algorithms are build on conforming reconstructions of the pressure and displacement together with equilibrated flux and stresses reconstructions. These ingredients are used to derive a posteriori error estimates for the fixed-stress algorithms, distinguishing the different error components, namely the spatial discretization, the temporal discretization, and the fixed-stress iteration components. Precisely, at the iteration k≥1 of the adaptive algorithm, we prove that our estimate gives a guaranteed and fully computable upper bound on the energy-type error measuring the difference between the exact and approximate pressure and displacement. These error components are efficiently used to design adaptive asynchronous time-stepping and adaptive stopping criteria for the fixed-stress algorithms. Numerical experiments illustrate the efficiency of our estimates and the performance of the adaptive iterative coupling algorithms.

An Introduction to Space–Time Exterior Calculus

The basic concepts of exterior calculus for space–time multivectors are presented: Interior and exterior products, interior and exterior derivatives, oriented integrals over hypersurfaces, circulation and flux of multivector fields. Two Stokes theorems relating the exterior and interior derivatives with circulation and flux, respectively, are derived. As an application, it is shown how the exterior-calculus space–time formulation of the electromagnetic Maxwell equations and Lorentz force recovers the standard vector-calculus formulations, in both differential and integral forms.