Electrostatic Approximation of Vector Fields

. This paper provides expressions for the boundary potential that provides the best electrostatic potential approximation of a given L 2 vector ﬁeld on a nice bounded region in R N . The permittivity of the region is assumed to be known and the potential is required to be zero on the conducting part of the boundary. The boundary potential is found by solving the minimization conditions and using a special basis of the trace space for the space of allowable potentials. The trace space is identiﬁed by its representation with respect to a basis of Σ − Steklov eigenfunctions.


Introduction
Quite often in physical applications one wishes to produce electrostatic fields in a region of known permittivity (x) that approximate a prescribed (given) field F. Such fields are determined by their boundary values, so a natural question is what imposed potentials on the boundary provide good approximations, in an energy norm, to F? Very often the boundary includes surface patches that are conductors as well as patches where nonzero potentials may be imposed.This may be regarded as a problem of stationary control or approximation.
The difficulty with such problems has been how to work with the control space of allowable boundary conditions as it will be a trace space of allowable H 1 functions.The standard Lions-Magenes description of trace spaces is not amenable to nice constructions of solutions for problems of this type.This problem is treated here using methods based on the spectral characterization of trace spaces as described in Auchmuty [3] and [4] which provides constructive methods and explicit bases for the traces.The allowable trace space is characterized as being isomorphic to a class of weak solutions of a linear elliptic equation and a basis of Σ−Steklov eigenfunctions of this space is identified.
Here an explicit expression for the best approximation in terms of the data is found.Moreover boundary data that is close in a boundary norm to this optimal solution will provide good approximations to the field F in an energy norm on Ω.The analysis described here is described for quite general N-dimensional regions since the results are essentially independent of the dimension N ≥ 2 and the methods may be of interest for other approximation questions.

Definitions and Requirements.
To analyze this problem, standard definitions, terminology and assumptions will be used as in Attouch, Buttazzo and Michaille [1].All functions in this paper will take values in R := [−∞, ∞], derivatives should be taken in a weak sense and N ≥ 2 throughout.
A region is a non-empty, connected, open subset of R N .Its closure is denoted Ω and its boundary is ∂Ω := Ω\Ω.Let L p (Ω), H 1 (Ω) be the usual real Lebesgue and Sobolev spaces of functions on Ω.The norm on L p (Ω) is denoted .p and the inner product on L 2 (Ω) by ., . .The basic requirement on Ω is (B1): Ω is a bounded region in R N whose boundary ∂Ω is the union of a finite number of disjoint closed Lipschitz surfaces; each surface having finite surface area.
The region Ω is said to satisfy Rellich's theorem provided the imbedding of The trace map is the linear extension of the map restricting Lipschitz continuous functions on Ω to ∂Ω.When (B1) holds, this map has an extension to W 1,1 (Ω) and then the trace of u on ∂Ω will be Lebesgue integrable with respect to σ, see [5], Section 4.2 for details.The region Ω is said to satisfy the compact trace theorem provided the trace mapping γ : on H 1 (Ω) and the associated norm is denoted u ∂ .This is an equivalent inner product to the usual inner product when Ω obeys (B1) -see [2] for a proof.Here ∇u := (D 1 u, . . ., D n u) is the gradient of the function u Our interest is in a problem that arises in electrostatics where part of the boundary Σ is a conductor and a potential can be imposed on the complementary part of the boundary Σ := ∂Ω \ Σ. Mathematically our requirements are (B2): Σ is an nonempty open subset of ∂Ω, Σ and Σ have strictly positive surface measure and σ(∂Σ 3

The Boundary Control Problem
The problem to be studied here is given a vector field F on Ω to find the potential ϕ that provides the best L 2 −approximation when the region Ω has known permittivity tensor (.) and part of the boundary Σ is a conductor held at zero potential.That is we want to find the function φ that minimizes with |.| is the Euclidean norm on R N .Since F is known this reduces to minimizing the functional where A(x) := (x) T (x) is real symmetric, G(x) := (x) T F(x) on Ω and the superscript T denotes the vector transpose.The following will be assumed.(A1): A(x) := (a jk (x)) is a real symmetric matrix whose components are bounded Lebesgue-measurable functions on Ω and there exist constants c 0 , c 1 such that Existence uniqueness and extremality conditions for this problem may be obtained using standard methods.The problem is a convex quadratic minimization problem on a Hilbert space so the existence may be stated as follows.
The functional E also is G-differentiable so the minimizers satisfy the following.
Theorem 32 Assume that (A1), (B1), (B2) hold and F ∈ L 2 (Ω : R N ) is given, then the minimizer φ of E on H 1 Σ0 (Ω) satisfies the equation To obtain further results about this problem a decomposition of the space H 1 Σ0 (Ω) will be used.Consider the bilinear form a : This bilinear form defines the a-inner product on H 1 Σ0 (Ω) and is equivalent to the ∂−norm (2.1) when A satisfies (A1).
Observe that a function w ∈ H 1 Σ0 (Ω) is a-orthogonal to H 1 0 (Ω) if and only if That is γ(w) is zero on Σ and L A w(x) := div (A∇w)(x) = 0 in a weak sense on Ω.The space of all such functions will be denoted N (L A , Σ) so we have the orthogonal decomposition where ⊕ a indicates that the a-inner product is used.
In light of this result, the minimizer φ has a decomposition of the form φ = φ 0 + φ b where φ 0 is the minimizer of E on H 1 0 (Ω) and (3.10)The fact that φ 0 is a solution of the extremality condition on H 1 0 (Ω) implies that div ( A∇ϕ 0 − G ) = 0 on Ω in a weak sense so this last equation may be written The solution of the problem for φ 0 is a standard Dirichlet boundary value problem and it is worth noting that the value of E(φ 0 ) = 0 if and only if div G = 0 on Ω in a weak sense.In this case the general problem reduces to that of solving (3.10) or (3.11) alone.
Our interest is in the problem of finding an expression for the boundary trace of φ or φ b on Σ.That is what boundary data gives the best approximating potential for the given field F on Ω?

Bases and Representations of N (L A , Σ)
To find the boundary data that provides the best L 2 −approximation to the field F, an orthogonal basis of the space consisting of certain Steklov-type eigenfunctions of L A is constructed and used.This will yield a spectral representation of φ b as described in theorem 5.xx below.
An a-orthonormal basis of N (L A , Σ) may be found using the algorithm described in Auchmuty [4].In the notation of that paper take V = N (L A , Σ), a as above and m(u, v) := Σ γ(u) γ(v) dσ. (4.13) This problem has the form of equation 2.1 of [4] and the bilinear forms a, m satisfy conditions (A1)-(A4) of that paper.Moreover condition (A5) there holds with H = L 2 ( Σ, dσ).
Define A(χ) = a(χ, χ), M (χ) := m(χ, χ) and C 1 to be the closed unit ball in H 1 Σ0 (Ω) with respect to the a-norm.Consider the variational problem of maximizing M on C 1 .This problem has maximizers ±χ 1 that have a-norm 1 and are solutions of (4.14) associated with an eigenvalue λ 1 > 0. Moreover one has the coercivity inequality Using the construction of section 4 of [4], a countably infinite a-orthonormal basis B := {χ j : j ≥ 1} of N (L A , Σ) may be constructed using a sequence of constrained maximization problems for M .

The best approximating potential
We are now in a position to specify the boundary data for the potential that minimizes (3.3).The preceding analysis enables the derivation of an explicit representation of the solution φ b of (3.10) or (3.11).Equation (3.10) implies that φ b satisfies Note that the g j depend onlds only on the data, the eigenfunction χj and the solution φ 0 of the zero-Dirichlet variational problem.Then the eigenfunction equation (4.14) yields that the solution is Theorem 51 Assume that (A1), (B1), (B2) hold and B, g j are defined as above.
Proof.This is a restatement of the preceding results.Note that when div G = 0 on Ω, then φ 0 = 0 so the last sentence holds.
Moreover when the boundary potential ϕ is a good approximation of this φ b in the norm of H 1/2 ( Σ) then the fact that the a-norm and the norm on H 1/2 ( Σ) are isometric implies that such potentials ϕ will provide a good approximation of F on Ω in any equivalent norm on H 1 Σ0 (Ω).
.11) This is the weak form of the equation L A ϕ = 0 on Ω subject to the boundary conditions ϕ(z) = 0 on Σ and A∇ϕ • ν = (G − A∇ϕ 0 ) • ν on Σ. (3.12) trace on Σ of the best approximation is given by the boundary trace of this right hand side.That is imposing Dirichlet boundary data γ(φ b ) on Σ given by (5.18) yields the minimizing potential φ b of E.