Coupling of External Electric Circuits with Computational Domains

Coupling of electrical circuits with 2D and 3D computational domains is very important for practical applications. To this aim, the notions of “electrical current” and “voltage” must be defined precisely and linked with local quantities (i.e., fields and potentials) in the computational domain. Apart from the static case, the definition of voltage is more complex than it may appear at a first glance, and it is usually tainted by unspoken and/or not justified assumptions. The purpose of this work is twofold: on one hand, to shed light on the definition and on the physical meaning of voltage in the case of time varying quasi-static fields and, on the other hand, to show how to establish coupling equations between lumped parameters circuit model and 2D/3D computational domains. It is demonstrated that a precise physical significance can be given to the voltage in terms of power balance only (the notion of potential is unnecessary). A couple of original operators which allow to express voltages and currents are introduced. Based on a critical analysis of the research literature, it is shown that existing coupling formulas can all be rewritten as particular cases of these two operators. The developed analysis is independent from any computational method and can be used to devise new coupling formulas.


Introduction
Coupling of external electrical circuits with Finite Element (FE) computational domains is very important for practical applications, and it has been investigated for a very long time (see [1] for a historical perspective in 1993). Coupling external circuits with 2D/3D computational domain is necessary to model electrotechnical appliances such as simple coils [2,3], non linear inductors [4], contactors [5,6], power transformers [7,8], electrical motors [9][10][11] and many others [12][13][14][15]. Maxwell's equations are solved in the computational domain by using various formulations, the unknown of which are local quantities such as potentials, or (less frequently) the magnetic fields directly [16][17][18][19]. Currents and voltages are global quantities that are applied to the computational domain through "ports", that is, interfaces between internal regions, or surfaces on the external boundary of the domain.
To this aim, the notions of "current" and "voltage" must be precisely defined and linked with the electric field e and current density j in the computational domain. The definition of the electrical current I which flows across a given surface Σ is easily expressed and it depends uniquely on the current density: where n is the unit normal vector, with outward orientation in respect of the computational domain. Conversely, the understanding of the physical significance of voltage is not as trivial as it may appear [20], and it is worth to be clarified. In the case of static fields J 2021, 4 (electrostatics, continuous currents), the electric field is conservative and the voltage U between two equipotential surfaces Σ a and Σ b can be uniquely defined as: where γ is any path which goes from Σ a to Σ b and v is an electric scalar potential such that e = − grad v. However, this definition of voltage cannot be truly satisfying in that it does not apply and cannot be easily extended to the case of time varying fields. In fact, in this case, the electric field is no more conservative and thus the choice of the path γ matters.
In the case of time varying fields, one may define the voltage as the circulation of the "electrostatic component" (− grad v) of the electric field: However, in the case of bounded domains, there is no such thing as the electrostatic component: at most, one may speak of an electrostatic component. In fact, Helmholtz's theorem [21] ensures that such a decomposition (3) exists, but in the case of bounded domains it is not unique [22] unless appropriate additional boundary conditions are imposed. Even assuming that such a definition is well posed from the mathematical point of view-that is, if one assumes that the value of the voltage hereby defined is independent on the gauge of (a, v), which is indeed the case-the question of the physical significance arises.
One observes that the definition (1) is not tainted by any of these issues: the current I is defined based on the current density only, and the definition is unambiguous in the case of time varying fields, and whatever the surface Σ.
The question of the well-posedness of the definition of voltage, in the case of time varying fields, is seldom addressed in undergraduate courses in physics and electrical engineering (because in practice "things work") but also in most of the scientific literature. In the author's opinion, any effort to shed light on this topic is worth to be explored.
The case of stranded coils [2,3,7,8] is considerably simpler. In fact, isopotential surfaces Σ a and Σ b reduce to points, so that in practice the integration path γ is unique, and it is imposed by the geometry of the electrical wires. In this case, the induced electromotive force can be expressed through the time derivative of the magnetic flux.
The case of solid conductors is more complex. In some cases [10], the voltage between two "isopotential" surfaces is mathematically defined as the difference of an electric scalar potential v, but the question of its physical significance is eluded. As already observed, this scalar potential is not unique and it is not clear at all why the voltage has the same value whatever the gauge of (a, v). The very same notion of "isopotential" surfaces is disturbing, because in the case of bounded domains it depends on a particular choice of the gauge, whereas the voltage ought to depend uniquely on the electric field.
The same considerations apply to [4,11,14], where a 2D modelling is taken on. In this case, Coulomb's gauge div a = 0 is arbitrarily fixed by imposing that the vector potential a is perpendicular to the xy plane, and thus a voltage between the front and the rear part of the domain can be uniquely defined as: where is the depth of the domain and n is the unit vector perpendicular to the plane. In [6,23], the following relationship is derived for a two-ports conductor in the electrokinetic case (also known as continuous current): where j 0 is any current density corresponding to a net current of 1 A. This equation is then generalised in various "flavours" [6,13,23], to the case of time varying fields, by rewriting grad v = −∂ t a − e: where j 0 = curl t 0 (see the aforementioned works for details). However, by doing so, one forgets that the starting Equation (5) holds under the hypothesis that fields do not vary with time. In spite of the fact that the derivation of (6)-(8) is somehow simplistic, all these equations are found to be correct. A completely different approach is to postulate that the current flows out of special thin "generator regions" Ω em f ,i , where a source of electromotive force exists [17,[24][25][26]. Inside the generator regions, the electric field is conservative and hence (2) is well posed, hence the voltage between the electrodes of a generator writes: where γ i is a path that joins the electrodes and e i is the electric field inside the generator.
In practice, the integral in (9) is never computed explicitly. Generators are removed from the computational domain so that their surfaces become new boundaries of the domain, over which appropriate boundary conditions must be imposed. Depending on the formulation, the voltage can be imposed either strongly (a formulation) or weakly (h formulation)-see details in [24,26]. This approach removes the practical difficulty of numerically defining voltages and has a correct power balance (all the electrical power which exits from generators enters in the domain), but it is not completely satisfying from the point of view of physics significance: voltage ought to be defined independently on how it is applied-that is, basing it exclusively on the electric field e inside the computational domain.
Finally, one observes that numerical modelling of the coupling of electrical circuits with computational domains is addressed in specific ways for each formulation. This gives rise to a large number of different formulas which express the link between global (currents and voltages) and local quantities (fields and potentials). This is a source of complication, including from the pedagogical point of view, because many different approaches ought to be introduced to students.
The purpose of this work is twofold: First, to clarify the notion of "voltage" in the case of time varying fields, and to expressing it by exclusively using the electric field inside the domain (without resorting to any potential), as it should be! Second, to devise a general method to write coupling equations between external electrical circuits and 2D/3D computational domains.
The article is organised as follows: an operator U ij [e] which expresses the voltage between two ports i, j for a given electric field e will be defined. It will be demonstrated that, by using this definition, the power balance is intrinsically respected: this gives a precise, physical significance to the notion of "voltage", even in the case of time-varying fields (Section 2). As a side product, a similar operator I n [h] for the current is obtained. It will be observed that the classical definition of current (1) can be rewritten as a particular case of usage of the new operator (Section 3). Then, the existing scientific literature is analysed critically and the most used Finite Element formulations for coupling external circuits and computational domains are reviewed in Section 4. It will be demonstrated that all of them can be retrieved by using the newly introduced operators U ij [e] and I n [h].
The obtained results are discussed in Section 5. Finally, the conclusions can be found in Section 6.

Definition of Voltage
In [27], Hiptmair and Sterz pointed out the difficulties to rigorously define what voltage is and suggested it could be defined through the notion of electric power. Based on this idea, it will be shown that it is possible to give a precise meaning to the notion of "voltage" in the case of time varying fields and for an arbitrary number of ports. The proposed definition of voltage respects the power balance and it is coherent with all previous works. Moreover, the proposed definition of voltage relies exclusively on the electric field; therefore, it is independent of any formulation or numerical method. To go further, the following theorem has to be demonstrated: Theorem 1 (fake power). Let e 1 and h 2 be an electric and a magnetic field, defined in a domain Ω c and Ω, where Ω c ⊂ Ω is the conductive part of Ω. Let b 1 and j 2 the magnetic flux density and the current density associated respectively with e 1 and h 2 , and assume that displacement currents are neglected: Assume that the boundary ∂Ω = Γ h ∪ Γ e can be parted in two surfaces Γ e and Γ h where the following boundary conditions hold: and that there exists a couple of potentials (a 1 , v 1 ) such that: Under this hypothesis, one has: 1.
The scalar potential v 1 is equipotential on each of the N connected components Γ e,n of Moreover, for a given electric field e 1 , the set of potentials {V 1,n } is unique up to an arbitrary constant value.

2.
The following equality holds: 3. Let I 2,n be the total current which enters in Ω c through Γ e,n : where n is the unitary vector normal to Γ e,n oriented outwardly. Then, the following equality holds: Before giving the demonstration, some remarks are mandatory: • It has to be highlighted that the hypothesis of the theorem does not require that constitutive laws are enforced (they could eventually be, but this is not necessary): In this case, the E side and the H side of Tonti's diagram are unlinked, thus e 1 and h 2 cannot be both the true electric field e and magnetic fields h: Hence, the quantity P [e 1 , h 2 ] defined by (17), which is homogeneous to a power, has no physical significance. For this reason, in the general case, P [e 1 , h 2 ] will be hereafter called "fake" power. Conversely, when constitutive laws are enforced, then e 1 = e and h 2 = h are the true electric and magnetic fields, and P [e, h] is the true, physically significant, instantaneous power which is provided to the domain by the generators. • The boundary Γ e writes: The surfaces Γ e,n are the ports which connect the computational domain with the external electric circuit. As such, they are also part of the boundary of ∂Ω c , where the scalar potential v 1 is defined hence. The surface Γ b is the remaining part of Γ e , where the scalar potential are not defined ( Figure 1). The boundary of Ω c writes: Γ j is the interface between Ω c and the remaining (insulator) part of Ω, where n · j = 0.
Proof of Theorem 1. By writing grad v 1 = −e 1 − a 1 and by using (15) on each port Γ e,n one has: Hence, v 1 is constant on each port Γ e,n . By using (14), the first integral of (17) writes: By integrating by parts, the second integral of (17) writes: where the boundary integral vanishes due to the boundary conditions (12) and (13). By using together (26) and (27), one obtains (18). Finally, by using integration by parts, The first integral on the right-hand side vanishes because div j 2 = div curl h 2 = 0. The boundary of Ω c can be parted as: where Γ j is the interface between the conductive and insulator part of the domain Ω i = Ω\Ω c . The current cannot flow in or out from Γ j , that is: hence the boundary integral in (28) writes: Finally, the uniqueness of the set {V 1,n } can be demonstrated by observing that these values depend uniquely on the electric field e 1 . Hence, a set of N − 1 independent equations can be obtained from (31) by selecting an arbitrary port Γ e,re f (namely Γ e,re f = Γ e,N ) and by considering N − 1 magnetic fields h 2,n such that I 2,re f = −I 2,n = 1 A.
In order to define the voltage U ij between two ports Γ e,i and Γ e,j , the set H ij of "test" magnetic field is defined: where I n is the net current due to j = curl h which enters in Ω through the n th port Γ e,n .
The following theorem provides a well-posed definition of the voltage U ij : Theorem 2 (definition of voltage). Let e be the electric field in Ω c , and Γ e,i and Γ e,j two ports.
There is one and only one value U ij [e] such that: for any h ∈ H ij .
Proof of Theorem 2. Theorem 1 holds with e 1 = e and h 2 = h , ∀h ∈ H ij , and thus the fake power writes: In order to demonstrate that this value is unique, it must be shown that it does not depend on h nor on the gauge of (a, v). The left-hand side term of (35) does not depend on the gauge because of (17), and the right-hand side of (35) does not depend on h . Hence, P [e, h ] is independent of both h and the gauge of (a, v). Therefore, P [e, h ] depends exclusively on the electric field e, which proves its uniqueness. Hence, This definition of voltage is well posed because U ij depends exclusively on the electric field e, and it holds as well for static and time varying fields. One observes that when V n = v| Γ e,n are the voltages of generators that feed the device (Figure 1), and e 1 = e and h 2 = h are the true electric and magnetic fields, the power balance is verified in that (20) writes: This result is important because it provides a precise physical significance to "voltages", that is the unique set of values U ij = V i − V j such that the power balance (36) is respected. In particular, in the case of a simple two-ports domain, the voltage can be defined like the ratio between the instantaneous power injected into the domain and the instantaneous current. One observes that the electric vector potential t 0 defined beforehand belongs to H ij : this shows that indeed (5) and (7) are particular cases of (34).

Expression of the Current
Now that a precise definition has been given to the term "voltage", one observes that the roles of e 1 and h 2 can be exchanged, so as to provide a useful expression of the net current which enters in the domain through electrodes. To this aim, let us define the set E n of "test" electric fields: with v ≡ 1 V on Γ e,n , v ≡ 0 on Γ e, =n } The following theorem provides a useful formula to compute the electric current which enters in the domain through a port: Theorem 3 (computation of current). Let h be the magnetic field in Ω, and Γ e,n a port. The electric current I n , which enters in the domain through Γ e,n , is equal to: for any e ∈ E n .
Proof of Theorem 3. Theorem 1 holds with e 1 = e and h 2 = h, ∀e ∈ E n , thus the fake power writes: In this case, there is no need to prove the uniqueness of I n because currents are correctly defined. Nevertheless, an argument similar to the one which has been used to prove Theorem 2 applies as well (that is, P [e , h] does not depend on the gauge, I n does not depend on e ).
One observes that even if (38) expresses the current based on h, only its curl matters: Finally, one observes that even the classical definition of current (1) can be rewritten as a particular case of (38). In fact, consider as test electric field the gradient of a scalar potential e = − grad v , the support of which is a thin layer of thickness w which lays on the surface Σ. Assume that v ≡ 1 V on Σ. In the limit of w → 0, the electric field writes e = (−n)δ Σ , where δ Σ is the Dirac distribution associated with the surface Σ. Hence, by using this particular test electric field with (38), the current writes:

Formulations
Hereafter, several Finite Element formulations are reviewed with the purpose of showing that, for all of them, the coupling equations with electrical circuits can be easily derived by using (34) and (38). In order to somehow simplify the notation hereafter, the divisions by 1 A or 1 V will be dropped and the expressions of the voltage and current will write simply: For each formulation, particular test fields h ∈ H ij or e ∈ E n will be chosen in order to simplify the writing of formulations at the discrete level. The purpose of this section is to show that known expressions to compute voltages and currents can be viewed as particular cases of (34) and (38), respectively.

Electrokinetics
The scalar potential electrokinetic formulation writes: that is div j = 0, with e = − grad v, with σ = electrical conductivity. The weak form writes: where v and v belong to appropriate function spaces. At the discrete level, the domain Ω = Ω c is meshed and the set of nodes N is parted as: where N i,n is the set of nodes on the n th port Γ e,n , and N v is the set of the remaining nodes. Classical nodal shape functions s n (x) are associated with each node of N v , whereas special isopotential shape functions f n (x) are associated with each port Γ e,n [24][25][26] (see Figure 4 of [25]): The support of f n is bound to the layer of elements which touch Γ e,n . One observes that f n ≡ 1 on Γ e,n and f n ≡ 0 on any other port, thus e n = − grad f n ∈ E n . By using (38) the net current I n which enters in Ω c through the n th writes: where the last equality comes from (45). One observes that the latter equation writes where h is the magnetic field (which is not computed by using this formulation) and j = curl h. The discrete approximation of v writes: With this formulation, the potential V n on the port Γ e,n can be imposed strongly, as it belongs directly to the set of unknowns. When Galerkin's method is used, the same shape functions are used both to expand the unknowns and as test functions; hence, the coupling equation which allows to impose the current I n appears naturally for v = f n : In order to illustrate this formulation with a practical example, consider a conductor (σ = 10 S/m) with a hole and a low conductive region (σ = 10 −6 S/m), where a continuous voltage of 5 V is imposed between the two electrodes on the left and on the right (Figure 2a). A test electric field e is computed by solving the electrokinetic formulation with a random conductivity (Figure 2b). A voltage of 1 V is imposed between the electrodes. A realisation of the test electric field e is depicted in Figure 2c, while the true current density j is depicted in Figure 2d. The total current is computed by using the classical expressions (1) and (49). It is observed that the results obtained by the two expressions converge to the same value as the mesh is refined ( Table 1). The corresponding program is written in the FreeFem++ environment [28], and it is freely available at https://github.com/scorretti/examples_ circuit_coupling (accessed on 23 November 2021).  Figure 2. The domain composed of a conductor with a poorly conductive region (green) and a hole (a). A random conductivity (b) is used to compute the test electric field e (c) that, on its turn, is used to compute the total current I corresponding to the real current density (d).

Eddy Current Formulations
A very large number of formulations have been developed to model eddy currents. In general, either the current or the voltage can be imposed strongly, whereas the other quantity has to be imposed weakly. Hereafter, some of the most used formulations are reviewed, but similar considerations can be taken on for any other formulation.
With the t − t 0 − φ formulation, the equations on the E side of Tonti's diagram (curl e = −∂ t b, div b = 0) are imposed weakly: with ρ = 1/σ = electrical resistivity and ν = 1/µ = magnetic reluctivity. The magnetic field h and the current density j write, respectively [13,23,29]: where φ is a magnetic scalar potential and t 0,n are precomputed electric vector potentials defined on Ω 0 ⊃ Ω c . For each two-ports solid conductor, a potential t 0,n is precomputed in such a way that the net current associated with j 0,n = curl t 0,n = 1 A (see Figure 5.29, p. 171 of [29]), thus t 0,n ∈ H n (in order to simplify the notation, where two-ports cases are concerned, one writes n instead of i, j). Currents in coils I n belong to the set of unknowns and are, therefore, imposed strongly. By taking h = t 0,n , one finds with (34) the expression of the voltage U n associated with each coil: In the case of nonsimply connected conductors, additional potentials t 0,k and associated unknowns I k must be added to cope with the so-called connexity problem (that is, to be able to impose Ampère's theorem). The corresponding voltages U k are set to 0 by using (55): By the way, one observes that in the case of static fields (∂ t ≡ 0) one seamlessly obtains the dual electrokinetic formulation [30], with respect of the aforementioned v formulation, by imposing curl e = 0: Voltages can be expressed by (34) as:

a − v Formulation
With the a − v formulation, the equations on the H side of Tonti's diagram (curl h = j, div j = 0) are imposed weakly: The electric field e and the flux density b write, respectively where a and v are, respectively, a magnetic vector potential and an electric scalar potential. By using (34) together with (20), one observes that voltages is determined by imposing the scalar potential v only. In fact, for any h ∈ H ij one has: For instance, in [25], this is accomplished by discretising the scalar potential v as: where the function v 0,n is any potential which corresponds to a voltage of 1 V for the n th solid conductor, and hence − grad v 0,n ∈ E n (in [25], only two-ports conductors are considered). It is found that the current I n which flows through the n th conductor writes: where e = − grad v o,n . It seems that, in order to retrieve (38), the term´Ω h · ∂ t b is missing. However, in fact, notice that − grad v o,n is defined only on Ω c and it is curl free. It is possible to find a curl-free extension of − grad v o,n over all the domain: whereṽ is any scalar potential which enforce the continuity with v on ∂Ω c . The thereby defined test electric field e is curl free and it is defined over all the domain; thus: Hence, the integral´Ω h · ∂ t b vanishes and the expression (38) is found again. Notice that the real electric field e is not curl free, but test electric fields e are allowed to.

h − φ Formulation
With the h − φ formulation, the equations on the E side of Tonti's diagram (curl e = −∂ t b, div b = 0) are imposed weakly. The computational domain is parted as Ω = Ω c ∪ Ω C c and a number of "cuts" have to be conducted to make Ω C c (insulators) simply connected. Assume that only two-ports solid conductors are present, some of which are coils. The field is impressed by generators through thin sections in the coils. The weak form of the formulation writes [17]: The unknown of this formulation is, directly, the magnetic field h, which is discretised as: where E c is the set of edges internal to Ω c , N C c is the set of nodes on Ω C c and on its boundary, and C is the set of cuts required to make Ω C c simply connected (see [17,24,31] for details). The key point is that j n = curl c n is a current density which corresponds to a unit net current in the n th coil, and hence c n ∈ H n . The coefficients I n represent the currents which flow through the n th coil, which can therefore be imposed strongly. It is found [17] that the corresponding voltage U n writes: where the test magnetic field is h = c n (see Figure 2 of [26]).

Stranded Coils
The case of stranded coils, composed of a multitude of thin insulated wires where eddy currents in individual wires can be neglected, is considerably simpler because the direction of j is known in advance. These kinds of coils are handled through a wire density vector j 0 which corresponds to a unit current. The homogenised current density in stranded coils writes: where ||j 0 || = n c /S c , with n c = number of turns, and S c = section of the coil [7]. Regardless of the formulation used, it is natural to impose currents which flow in stranded coils because the variable I appears explicitly in (70). Voltage in stranded coils writes: where the resistance R is known. Depending on the used formulation, the linkage flux Ψ may be written in different ways [7,12,26]: where curl t 0 = curl h s = j 0 (see the aforementioned works for details). By observing that RI =´Ω c j 0 · e and by using the latter equation, (71) boils up to one of (6) and (7), which are particular cases of (34) with h = t 0 or h = h s (cf. end of Section 2).

Discussion
In this work, an original couple of operators U ij [e] and I n [h] are introduced. These operators provide a rigorous definition of voltages and currents in terms of the electric and magnetic field only.
The notion of voltage is analysed and a general physical interpretation is given based on the electrical power balance of the computational domain. In the static case, the classical definition of voltage (2) is obtained as a particular case of (34), by taking as test magnetic field h γ , the Biôt-Savart field generated by a unit current which flows along the path γ: In the case of time varying fields, by taking the same test magnetic field h γ , the integral´Ω h γ · ∂ t b remains but the voltage can still be expressed as the difference of the scalar potential on the two (isopotential) ports, provided that the hypothesis of Theorem 1 is respected: Notice that (74) does not contradict the fact that, in the case of time varying fields, the circulation of e depends on the particular path γ.
It has to be remarked that not all couples (a 1 , v 1 ) satisfy the hypothesis of Theorem 1. A notable example of gauge where the hypothesis is not satisfied is the so-called temporal gauge [27] (also called Weyl gauge): This potential a * 1 is sometimes called modified vector potential [32]. Luckily, this is not a problem because, in order to demonstrate the well-posedness of the definition of voltage, it is enough that a single couple of potentials which satisfies the hypothesis of the theorem exists; the uniqueness of voltages is then demonstrated. Moreover, potentials are only intermediate actors in the demonstration and, in practice, it is not necessary to compute them; it is enough that at least a couple of potentials exists. Notice also that, in the case of unbounded domains, Helmholtz decomposition ensures, by itself, the uniqueness of the scalar potential [22], and thus of voltages.
Finally, notice that Theorem 1 requires that the electric field is normal to the ports which connect the computational domain to electric circuits. This is a limitative hypothesis which is required to define voltages but not for currents. Perhaps it is possible to weaken the hypothesis of the theorem so as to define voltages between arbitrary ports (that is, where the electric field is not necessarily normal) based on the power balance. Moreover, it can be conjectured that, in the case of physical electric fields, it is always possible to find at least a couple of potentials (a 1 , v 1 ) which satisfy the hypothesis of Theorem 1.
The findings on existing Finite Element formulations are resumed in

Conclusions
In this work, the existing scientific literature on methods for coupling external electric circuits with the Finite Element method is critically analysed. It is observed that, in the cases of time varying fields, the definition of "voltage" is at best unclear and usually tainted with unspoken and/or unjustified assumptions.
In order to overcome these problems, a couple of original operators are defined so as to provide workable expressions for voltages and currents. In particular, it has been possible to define voltages based uniquely on the electric field-that is, without relying on any potential. By doing this, a clear physical significance is given to voltages, based on the fact that power balance must be respected. This is a first contribution of this work.
This result restores the "symmetry" between current and voltage, in that current is defined based on the current density only (or, equivalently, on the magnetic field). However, this symmetry is not complete; currents can be uniquely defined through any surfaces (also known as ports), whereas in the case of voltages, the choice of surfaces is constrained by the fact that the electric field must be perpendicular. In the case of static fields, this constraint leads to the notion of equipotential surface (which is flawed for time varying fields); if it were possible to remove this constraint, it would be possible to define a "voltage" between two non-equipotential surfaces! Notice that this is not a problem from the practical point of view because the details of connections between external circuits and computational domains are never modelled accurately in the quasi-static regime.
By using the newly defined operators, it is possible to devise a general method to write down coupling equations between external circuits and computational domains.
In particular, it is demonstrated that all the coupling expressions found in the reviewed scientific literature for the Finite Element method can be retrieved as particular cases of these operators.
It is important to observe that, even if the analysis of existing formulations has been taken on with the Finite Element method only, these new operators do not rely on any particular numerical method. Hence, in principle, they can be used to devise the implementation of coupling between external electrical circuits and computational domains with any other numerical method (cell method, finite difference method, or whatever).
The major limitation of this work is that it is limited to the quasi-static case; in particular, displacement currents are completely neglected. It would be interesting to extend this work in order to take into account capacitive effects [33,34] or even wave propagation.
Funding: This research received no external funding.