Abstract
Systems of partial differential equations (PDEs) comprising a combination of constraints and evolution equations are ubiquitous in physics. For both theoretical and practical reasons, such as numerical integration, it is desirable to have a systematic understanding of the well-posedness of the Cauchy problem for these systems. In this article, we first review the use of hyperbolic reductions, where the evolution equations are singled out for consideration. We then examine in greater detail the extensions, namely, systems in which constraints are evolved as auxiliary variables alongside the original variables, resulting in evolution systems with no constraints. Assuming a particular structure of the original system, we provide sufficient conditions for the strong hyperbolicity of an extension. Finally, this theory is applied to the examples of electromagnetism and a toy model of magnetohydrodynamics.
1 Introduction
In this work, we continue [1–4] the study of first-order systems of equations in which there are more equations than unknowns, but with a structure that permits, in principle, splitting suitable linear combinations of them into “evolution” and “constraint” equations. We restrict to the case of consistent systems, in which the number of equations is equal to the number of constraints plus the number of independent variables, and furthermore to the special case in which the number of independent variables matches the number of evolution equations. The latter means that we do not consider systems with gauge freedom remaining, which would imply the existence of variables with unspecified equations of motion. In this case, one can attempt a solution by carefully restricting the initial data and then directly solving the evolution equations. For an introductory review, see Hilditch [5]. One must then check that the constraint equations are satisfied in the time development. For this, integrability identities among the whole system of equations must be satisfied. These conditions will be assumed and spelled out in detail below. This “free evolution approach” requires us to establish the well-posedness of the Cauchy problem Gustafsson et al. [6]; Kreiss and Ortiz [7] (for a review of well-posedness applied to general relativity, see Sarbach and Tiglio [8]). We restrict ourselves to the concepts arising from the theory of strongly hyperbolic systems, in which well-posedness is determined by algebraic properties of the principal symbol of the equation system. For first-order systems, the principal symbol is simply the set of matrices multiplying the derivatives of the variables. The algebraic properties leading to well-posedness have several equivalent characterizations summarized in the Kreiss matrix theorem Kreiss [9]. To assert well-posedness for the systems under consideration, we need to find a suitable square system, that is, a system where the number of variables equals the number of equations. This can be achieved by taking a subset of the equation system, called a reduction, resulting in a pure evolution system. The use of reductions is customary, but another possibility, which is often employed in numerical schemes, consists of making an extension, that is, extending the system by adding more variables. These extensions are commonly referred to as divergence cleaning [10]; Munz et al. [11, 12], from their use in magnetohydrodynamics, or as [13] or Z-systems [14] from their use in general relativity.
A paradigmatic example is given by the Maxwell equations,where the unknowns are the components of the Faraday tensor , an anti-symmetric tensor (so there are a total of six independent variables). , the current vector, is a given vector fixed in space-time, which has vanishing divergence. This is necessary due to the integrability identity . We work here in four-dimensional space-time with the Levi–Civita derivative associated with . There are thus a total of equations for , so six of them should be evolution equations, and the remaining two should be constraints. Introducing a time-like covector , one finds that contraction with that vector on both equations gives constraint, that is, equations which have derivatives only in directions perpendicular to ; while projection on the space perpendicular to gives equations that have derivatives along for each of the independent components of . Thus, in the terminology introduced above, a reduction is obtained by taking only these projections as the evolution equations. The integrability identity and divergence property of together imply that constraints are satisfied in the time development if they are at an initial surface.
On the other hand, an extension is given by adding two auxiliary constraint variables , one for each Maxwell constraint, and making a choice for their equations of motion. To accomplish this in a covariant fashion, we need to define two tensor fields . The proposed extended system is
It turns out that if the symmetric parts of are Lorentzian metrics whose cones have non-zero intersections among each other and with the cone of , then the extended system is well-posed. (We use the mathematical notion of a cone; when needed, we use the term light cone to refer to their boundaries). The equations that were constraints are now evolution equations for , and the others acquire spatial derivatives of these fields. As mentioned above, such extensions have been employed with enormous success in numerical relativity [15–20] and computational astrophysics, with works introducing this approach for magnetohydrodynamics [11, 12]; Dedner et al. [10] is particularly influential. Here, we investigate the space of possible extensions that lead to well-posed Cauchy problems and how to construct them in a natural, covariant fashion.
The article is organized as follows. In Section 2, we define the type of systems to be considered, including the necessary conditions they must satisfy in order to have a well-posed Cauchy problem. In Section 3, we introduce the Kronecker decomposition of matrix pencils and explain its implications to the study of strongly hyperbolic systems. In Section 4, we formalize the framework for extensions. Given the considerable freedom in choosing them, we use the Kronecker decomposition as a guide for making these choices. In Section 5, we demonstrate how this framework applies to two concrete examples: Maxwell’s electrodynamics and a toy model of magnetohydrodynamics (MHD). Finally, in Section 6, we conclude with discussions and provide comments on how this line of research is being further developed.
2 Preliminaries and notation
To fix notation, we specify the systems we consider, following the notation of Geroch [1]; Abalos and Reula [3]; Abalos [4]. We consider a manifold of dimension , and the following system of the quasi-linear first-order partial differential equations on the fields ,where the indices , , are abstract, grouping several tensorial indices into one and merely indicating where the contractions are. We use lower-case Latin indices to denote single vector indices, lower-case Greek indices to indicate variable fields, and upper-case Latin to label the equations space. The function on indices indicates their total dimension.
We impose the following conditions on :
Condition 1: the generalized Kreiss condition.
We assume that the matrix is smooth in all arguments and that there exists a hypersurface orthogonal covector such that for all values of , not proportional to , the matrix pencilhas a kernel only for a finite set of real values of (the term matrix pencil refers here to the uni-parametric combination , where and are matrices that do not depend on the parameter ).
In addition, the corresponding singular values of approach zero in a linear way, that is, , with in a neighborhood of . We recall that the singular values are the square roots of the eigenvalues of . Because this is an matrix, there are singular values (see Abalos [2] for more details and for a more general definition).
These conditions imply two things: i) the rank of is maximal. Therefore, by defining any vector transversal to the surface flat defined by (i.e., ), we can obtain all field derivatives along from their values and their derivatives at that surface. This means that we have enough evolution equations for each field . Observe that once we have a choice of satisfying Condition 1, then there is an open set of covectors satisfying the same condition. Thus, we can always form hypersurfaces in a neighborhood of any point, leading to a local initial value problem; ii) In the case that the number of equations equals the number of variables, these conditions imply there is a well-posed Cauchy problem, in the usual sense for strongly hyperbolic systems, off of the mentioned surface. This is the classic Kreiss condition.
In case there are more equations than variables, we need to make sure that there are no more linearly independent equations having derivatives along the transversal vector ; otherwise, we would have an inconsistency because two equations could give different values for the same transversal derivative. To accomplish that, we impose:
Condition 2: the Geroch constraint condition.
If the number of equations is larger than the number of variables , then we assume there exists a set of matrices , which are labeled by upper-case Greek indices, withand that . This condition ensures that the rest of the equations do not have derivatives off of the surface defined by , so that the system is consistent. Indeed, the following linear combination of equations, called constraints,have only derivatives on the flat defined by .
There is a further consistency condition that would guarantee that if the initial data are such that constraint quantities vanish at the initial surface, then they would also vanish along evolution [4]. We require the following:
Condition 3: integrability.
In other words, there is a particular on-shell identity among derivatives of our equation system. In most cases of physical interest, this identity is a consequence of gauge or diffeomorphism invariance.
3 Kronecker decomposition
When studying the well-posedness of the Cauchy problem, the relevant aspect is the behavior of the system in the limit of high frequencies. We can thus restrict our attention to a neighborhood of each point and work in the frequency domain, employing the Fourier–Laplace transform in space and time, respectively. Explicitly, we consider a time function and a foliation given by its level surfaces. We define and take a vector transversal to the foliation and adjust it such that . We choose covectors such that and define . We perform Fourier in , and Laplace in . Thus, we replace space derivatives by and time derivatives by . Furthermore, in what follows, once any particular is chosen, we take a coordinate base so that , and , and so . Finally, in the high frequency limit, we obtain .
The Kronecker decomposition of a matrix pencil is a canonical transformation that generalizes the Jordan decomposition of a square matrix pencil. Considering the (square or non-square) pencil , the Kronecker decomposition is achieved by multiplying this pencil by specific matrices and , which are independent of (as in the square Jordan decomposition case). This transformation results in a new pencil that has a block structure with particular canonical blocks (see Gantmakher [21, 22], for a detailed description and Equation 3 for an example).
It turns out that the Kronecker decomposition can be used naturally in the analysis of systems with constraints or gauge freedom. With the first two conditions assumed above, the Kronecker decomposition of the pencil is given by
Ultimately, this represents a change of basis of both the variable and equation spaces, which depends on but not on . The first block is a diagonal block, this diagonal represents the true degrees of freedom of the entire system. It contains as many elements as the “zeros” of the singular value decomposition, counting their multiplicity. The blocks, called in the literature, are due to the constraints; there are a total of blocks. Because each block occupies two rows, we see that the number of zero rows is . The zero rows are present in many systems; they represent differential constraints among the constraints themselves. The numbers defined above also satisfy:
With this decomposition at hand, it is easy to see how to choose among them linear combinations that give evolution equations for all . Observe that the equations (rows) with a are certain to contain derivatives transversal to the flats. So, we must include them, but we can add any combination of the other rows to them. It turns out that, by simply adding to each of these rows the immediate row below, multiplied by any number , and discarding all the remaining rows, we obtain the evolution equations.Thus, we have constructed a map from the equation space to the variable space, which we refer to as a reduction and denote by . Thus, is a map from the variable space into itself consisting of a set of diagonal matrices satisfying the classic Kreiss conditions (see point ii. within Condition 1). Notice that we can choose the extra roots of (i.e., the ) as we please. They are the propagation speed of extra constraint modes. This simple observation is the principle behind the results in Reula [23]; Abalos [4].
Thus, there is a reduction (a linear combination of the equations) such that the Cauchy problem of the system is well-posed. Furthermore, Condition 3 asserts that if the initial data satisfy all equations (including the vanishing of the constraints), then all the equations are satisfied for all times as long as the solution exists. See Abalos [4] for details.
4 Extensions
A generic extension would imply the addition of an extra matrix, (and extra variables ), to obtain a square systemHere, is a term we can also freely choose that does not include derivatives of or and that goes to 0 when Z goes to 0. In general, represents damping terms [13]; [10]; [24], which are important in numerical applications. For simplicity in our discussion, however, we omit it.
Because we are interested in solving Equation 2 for , our extension proposal only makes sense if we can show that for suitable initial data (for ), the solution of Equation 4 has throughout the development, thereby ensuring that is a solution of Equation 2.
As we explained before, if we assume Conditions 1, 2, and 3 and take any initial data for satisfying the constraints, we know that the initial value problem for Equation 2 is “well-posed” and has a unique solution . (Here, by well-posed, we mean that the map from Cauchy data to solutions is continuous. To establish this, one finds a hyperbolic reduction from which we may assert that the reduced system is well-posed for arbitrary initial data. Then, one shows that if the initial data satisfy the constraints, then the solutions of the reduced system also satisfy them. Thus, they are solutions to the whole system, and we call the whole system well-posed). Therefore, if we choose such that the extended system, Equation 4, is well-posed, then for any initial data, there will be a unique solution. If we choose as initial data , then will be a solution, and by uniqueness is the solution. Therefore, we only need to show that system Equation 4 satisfies Kreiss’s condition.
4.1 Strong hyperbolicity of the extensions
A particularly interesting set of extensions is obtained by noticing the symmetry between the Kronecker decomposition of and . So, we start by computing it:
Recalling that the matrices can be thought of as a basis, labeled by , for the kernel of , it is easy to understand its structure. Here, the rows with zeros are in number. This is so because the diagonal part of cannot contribute to the kernel. We then have blocks , observing that they have a minus sign on them. This is because they are kernels for the corresponding blocks of . Finally, there is a block that is a kernel of the zero rows of . This part is completely undetermined, so we have simply added a diagonal matrix.
To make more apparent the extension we proposed, we reorganize the rows of and such thatHere, all the matrices are blocks matrices where of size , of size , and is the identity matrix of size . The zero rows of are of size , and the zero rows of are of .
From this reorganization, it is apparent that a natural choice of is given bywhere now must be chosen to render the system diagonalizable. This is, of course, not the most general extension but is a natural and fully covariant proposal for . The principal symbol of Equation 4 becomes thena square matrix.
We now propose a particular expression for , namely,where is of size , and is the identity matrix of size .
Using expressions Equations 5, 6, we concludeIt is easy to verify that this matrix is pencil-similar to the following diagonal matrix:and so it satisfies Kreiss’s condition. The extra eigenvalues , introduced by , come in pairs, which means that there are new null cones as characteristic. We shall see this in the examples below, where Lorentzian metrics are used to realize these null cones.
5 Examples
In this section, we present two implementation examples of our proposal, showing that they produce well-posed systems while largely preserving the covariance of the original theories. In all cases, extra Lorentzian metrics are introduced to avoid light cone intersections.
5.1 Maxwell’s equations
We start with the example given in the introduction Equation 1. For them, we have a space of variables (anti-symmetric tensors), which is dimensional in a four-dimensional space-time of metric . The space of equations is , namely, two space-time vectors. We have (see Geroch [1])
Given a time-like , we have
So, it is the map , which is of the maximal rank. This system satisfies Condition 1; see Abalos and Reula [3] for more details.
The tensor is also of maximal rank for any b1. Since the dimension of the image is 2-dimensional, we have , and the system is consistent, satisfying Condition 2.
We also haveand so Condition 3 is also satisfied.
A suitable reduction isThis renders the evolution equations symmetric hyperbolic. As we saw above, a simple extension is obtained introducing two tensors and defining
If we take their symmetric parts to be any two Lorentzian metrics, each one of them sharing a common time-like covector with , but not touching their null cones (for brevity, we do not consider here such degenerate cases), then the system is strongly hyperbolic and so has a well-posed Cauchy problem. To check this, we compute the characteristics of the system and the corresponding eigenvectors and see when we get a complete set, that is, a total of eight eigenvectors.
The characteristic equations arewhere we need to solve these equations for with and given and for the eigenvectors and . The solutions split into three cases: first, when is null with respect to (physical case), then when it is null with respect to or (extended cases), as we explain below.
We already know four of the eigenvectors, namely, the physical ones arising from the original system. To recover these, we set and search for the value of . The second equation then implies that for some vector , while the first implies that where indices are raised with the space-time metric. Because cannot be proportional to (otherwise would vanish), both terms must vanish and so we concludewhich admits two real solutions for . Hence, is orthogonal to , which leaves two options remaining for for each of the two values of .
Now, we want to find the rest of the eigenvectors. For that, we first choose , . Contracting the first equation with , and using the anti-symmetry of , we get a condition for ,which again admits two real values of . Repeating the argument above, the first equation becomes
Because the null cones of and are by assumption not touching, we have . It follows that satisfies Equation 8 provided that Equation 7 holds. Observe furthermore that satisfies the same equations and results in the same Faraday tensor for any . Thus, Equation 7 gives two additional eigenvectors.
If we drop the assumption that the null cones of and are non-touching and assume that they touch at , then to have a solution, we need that must be proportional to .
The final case is similar to the second. We choose , and obtainand the same equations for the dual of , so we need not discuss it separately.
In summary, we have obtained the eight eigenvectors we require to satisfy the Kreiss condition and conclude that the system is strongly hyperbolic.
5.2 Toy MHD
Here we look at the evolution of a magnetic field driven by a given velocity field in a space-time . The system isHere, we take to be time-like and of norm one, . We also take . This last is a gauge condition to make the solutions unique for the whole system because otherwise, if is a solution, then also is a solution, with an arbitrary function.
We observe that there are four equations for three variables. Three of them are evolution equations for the three components of . We shall see below that the other is a constraint. Thus, Condition 2 is also satisfied.
The principal part of system Equation 9 is
It is easy to check that Condition 1 is satisfied. The Geroch matrices are also easy to obtain as . They form a basis of the left kernel of and, as we explained before, this means that when Equation 9 is contracted with , a constraint is generated; this iswhere . We notice that this is the spatial divergence of in disguise.
On the other hand, the following integrability condition holds trivially; thus, the system satisfies Condition 3.
The extended system consists of adding a term to Equation 9, with as in the previous example and with the extra variable . Its principal part is , with . The characteristic equation iswhere we need to solve this equation for with given, and for the eigenvectors and (with ).
Without loss of generality, we choose such that , and we rewrite the characteristic equations projecting on to and perpendicular to it (with the projector ). We obtain
The physical solution comes from choosing , and the eigenvectors and orthogonal to . Because has two possible directions, we obtain two eigenvectors.
The remaining eigenvectors come from choosing such thatand , . This expression satisfies the second characteristic equation trivially, and it is easy to verify that the first one reduces to
Because, as before, there are two solutions for from Equation 11, we obtain two more eigenvectors. In summary, we have obtained the four eigenvectors we require to satisfy the Kreiss condition and conclude that the extended system is strongly hyperbolic. Finally, we notice that Equation 11 can also be rederived from the integrability condition, i.e., by multiplying Equation 10 by .
6 Conclusion
Similar extensions to those proposed here were previously known, starting with the divergence cleaning used in magnetohydrodynamics and later generalized as -systems for generic symmetric hyperbolic systems. To implement them, it was necessary to break the covariance of the system in the usual sense of performing a decomposition. For symmetric hyperbolic systems, such extensions can be obtained in our framework by committing to a frame and a reduction and adding an extra term that annihilates the time component of the constraint basis. This results in an extended symmetric hyperbolic system.
In this article, we have presented an extension scheme for first-order PDEs. With appropriate adaptation, however, these results can be applied to systems of two or even more orders. We will show in future articles how to apply these ideas to gravity theories to extend the system and to fix the gauge, allowing us to reinterpret and generalize known results such as those of Bona et al. [25]; Hilditch and Richter [26]; Kovács and Reall [27].
Although the existence of a strongly hyperbolic extension is performed in Fourier space and results in a system of pseudodifferential equations, our examples show that in cases of physical interest, one may obtain differential extensions. These extensions furthermore retain covariance of the theory in the sense that, contrary to earlier -system extensions, at least in the principal part, they do not rely on a preferred time direction but instead the addition of other Lorentzian metric tensors. Further details and a complete proof will be provided in a longer version of this work.
In our analysis, we resorted to previous work to argue that the constraints, if initially satisfied, are satisfied at later times. This helped us conclude that remains zero throughout the evolution. There are, however, more elegant ways to show this when the constraints do not have any kernel from the left, that is, no set of zero rows in their Kronecker decomposition (see Equation 3). In such cases, it can be shown that the fields satisfy a second-order evolution system that is decoupled from and has a well-posed initial value problem. Choosing these fields to vanish at the initial surface and the fields satisfy the original constraints of the system, all derivatives of vanish on the initial surface, in particular any transversal derivative, so the unique solution to the second-order system is 0, and the constraints are satisfied for all times. Unfortunately, the presence of zeros may prevent the second-order system from being well-posed, so more care is needed. This will be further considered in the aforementioned longer article.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
FA: conceptualization, formal analysis, investigation, methodology, project administration, validation, writing–original draft, and writing–review and editing. OR: conceptualization, formal analysis, methodology, validation, writing–original draft, and writing–review and editing. DH: investigation, validation, writing–review and editing, and formal analysis.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Grant PID2022-138963NB-I00 funded by MCIN/AEI/10.13039/501100011033/FEDER, UE, by CONICET, SeCyT-UNC, and MinCyT-Argentina and by FCT Project No. UIDB/00099/2020.
Acknowledgments
We thank Carlos Palenzuela for several helpful discussions associated with the article. We also thank the “Programa de projectes de recerca amb investigadors convidats de prestigi reconegut” of the Universitat de les Illes Balears through which Oscar Reula visited the UIB, making it possible to finalize this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
1.^Here the target space is two copies of , and the image is 1-dimensional on each one of them.
References
1.
GerochR. Partial differential equations of physics. Scotland: General Relativity, Aberdeen (1996). p. 19–60.
2.
AbalosF. A necessary condition ensuring the strong hyperbolicity of first-order systems. J Hyperbolic Differential Equations (2019) 16:193–221. 10.1142/s0219891619500073
3.
AbalosFReulaO. On necessary and sufficient conditions for strong hyperbolicity in systems with constraints. Class Quant Grav (2020) 37:185012. 10.1088/1361-6382/ab954c
4.
AbalosJF. On constraint preservation and strong hyperbolicity. Classical Quan Gravity (2022) 39:215004. 10.1088/1361-6382/ac88af
5.
HilditchD. An introduction to well-posedness and free-evolution. Int J Mod Phys A (2013) 28:1340015. 10.1142/S0217751X13400150
6.
GustafssonBKreissH-OOligerJ. Time dependent problems and difference methods, 24. John Wiley and Sons (1995). 10.1002/9781118548448
7.
KreissH-OOrtizOE. Introduction to numerical methods for time dependent differential equations. John Wiley and Sons (2014).
8.
SarbachOTiglioM. Continuum and discrete initial-boundary value problems and einstein’s field equations. Living Rev relativity (2012) 15:9. 10.12942/lrr-2012-9
9.
KreissH-O. Über die stabilitätsdefinition für differenzengleichungen die partielle differentialgleichungen approximieren. BIT Numer Mathematics (1962) 2:153–81. 10.1007/BF01957330
10.
DednerAKemmFKrönerDMunzC-DSchnitzerTWesenbergM. Hyperbolic divergence cleaning for the mhd equations. J Comput Phys (2002) 175:645–73. 10.1006/jcph.2001.6961
11.
MunzC-DOmnesPSchneiderRSonnendrückerEVoßU. Divergence correction techniques for maxwell solvers based on a hyperbolic model. J Comput Phys (2000) 161:484–511. 10.1006/jcph.2000.6507
12.
MunzC-DOmmesPSchneiderR. A three-dimensional finite-volume solver for the maxwell equations with divergence cleaning on unstructured meshes. Computer Phys Commun (2000) 130:83–117. 10.1016/s0010-4655(00)00045-x
13.
BrodbeckOFrittelliSHübnerPReulaOA. Einstein’s equations with asymptotically stable constraint propagation. J Math Phys (1999) 40:909–23. 10.1063/1.532694
14.
BonaCLedvinkaTPalenzuelaCZacekM. General covariant constraint free evolution system for numerical relativity. arXiv (2002). 10.48550/arXiv.gr-qc/0209082
15.
BaumgarteTWShapiroSL. Numerical integration of Einstein’s field equations. Phys Rev D (1998) 59:024007. 10.1103/PhysRevD.59.024007
16.
ShibataMNakamuraT. Evolution of three-dimensional gravitational waves: harmonic slicing case. Phys Rev D (1995) 52:5428–44. 10.1103/PhysRevD.52.5428
17.
NakamuraTOoharaKKojimaY. General relativistic collapse to black holes and gravitational waves from black holes. Prog Theor Phys Suppl (1987) 90:1–218. 10.1143/ptps.90.1
18.
PretoriusF. Numerical relativity using a generalized harmonic decomposition. Class Quant Grav (2005) 22:425–51. 10.1088/0264-9381/22/2/014
19.
BernuzziSHilditchD. Constraint violation in free evolution schemes: comparing the BSSNOK formulation with a conformal decomposition of the Z4 formulation. Phys Rev D (2010) 81:084003. 10.1103/PhysRevD.81.084003
20.
AlicDBona-CasasCBonaCRezzollaLPalenzuelaC. Conformal and covariant formulation of the z4 system with constraint-violation damping. Phys Rev D (2012) 85:064040. 10.1103/physrevd.85.064040
21.
GantmakherFR. The theory of matrices, 1. New York: Chelsea Publishing Company (1998).
22.
GantmakherFR. The theory of matrices, 2. Providence, RI: American Mathematical Soc. (1998). p. 131.
23.
ReulaOA. Strongly hyperbolic systems in general relativity. J Hyperbolic Differential Equations (2004) 1:251–69. 10.1142/s0219891604000111
24.
GundlachCCalabreseGHinderIMartín-GarcíaJM. Constraint damping in the z4 formulation and harmonic gauge. Classical Quan Gravity (2005) 22:3767–73. 10.1088/0264-9381/22/17/025
25.
BonaCLedvinkaTPalenzuela-LuqueCZacekM. Constraint-preserving boundary conditions in the Z4 numerical relativity formalism. Class Quant Grav (2005) 22:2615–33. 10.1088/0264-9381/22/13/007
26.
HilditchDRichterR. Hyperbolicity of physical theories with application to general relativity. Phys Rev D (2016) 94:044028. 10.1103/PhysRevD.94.044028
27.
KovácsADReallHS. Well-posed formulation of lovelock and horndeski theories. Phys Rev D (2020) 101:124003. 10.1103/PhysRevD.101.124003
Summary
Keywords
well-posed initial value problem, constraint equations, evolution equations, extensions, singular value decomposition (SVD), Kronecker decomposition, electromagnetism, magnetohydrodynamics
Citation
Abalos F, Reula O and Hilditch D (2025) Hyperbolic extensions of constrained PDEs. Front. Phys. 12:1517192. doi: 10.3389/fphy.2024.1517192
Received
25 October 2024
Accepted
30 December 2024
Published
17 February 2025
Volume
12 - 2024
Edited by
Jose Luis Jaramillo, Université de Bourgogne, France
Reviewed by
Guangqing Feng, Henan Polytechnic University, China
Philipp OJ Scherer, Technical University of Munich, Germany
Updates
Copyright
© 2025 Abalos, Reula and Hilditch.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fernando Abalos, j.abalos@uib.es
ORCID: Fernando Abalos, https://orcid.org/0000-0001-7863-3711; Oscar Reula, https://orcid.org/0000-0003-2517-7454; David Hilditch, https://orcid.org/0000-0001-9960-5293
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.