Radiative, Reflective, and Dissipative Boundaries

A boundary condition can be mathematically stable yet physically reflective, or physically suggestive yet fail to control the boundary energy. These are separate tests. Maximally dissipative conditions control the incoming characteristic fields with a contraction bound. Their reflection coefficient depends on which fields were diagonalized, the location and curvature of the boundary, background scattering, and the asymptotic expansion used to model radiation. No local condition at a generic finite radius is an exact substitute for the missing exterior Einstein solution.

Required background. The Gibbons–Hawking–York Boundary Term supplies timelike-boundary geometry, while PDEs and Green Functions for Gravitation supplies characteristic energy flux, Sommerfeld conditions, and stability estimates.

Helpful background. The Einstein Equation supplies gravitational radiation and constraints; Hypersurface Normals and Tangential Projectors supplies outward orientation; Formulations of General Relativity orients the volume; and Causality and Global Structure supplies the exterior domain omitted by a finite boundary.

Boundary energy and reflection coefficient

Consider the scalar wave equation on $x\le0$ with boundary at $x=0$ and outward normal toward increasing $x$:

$$u_{tt}-u_{xx}=0.$$

Define the outward and inward characteristic derivatives

$$w_{\mathrm{out}}=u_t-u_x, \qquad w_{\mathrm{in}}=u_t+u_x.$$

The energy

$$E=\frac12\int_{-\infty}^0(u_t^2+u_x^2)\,dx$$

satisfies

$$\frac{dE}{dt} =u_tu_x\big|_{0} =\frac14 \left(w_{\mathrm{in}}^2-w_{\mathrm{out}}^2\right).$$

For

$$w_{\mathrm{in}}=R\,w_{\mathrm{out}}+q,$$

homogeneous data with $\lvert R\rvert\le1$ do not inject energy. The cases $R=-1$ and $R=1$ are perfectly reflective Dirichlet and Neumann conditions; $R=0$ is exactly absorbing for this one-dimensional outgoing mode; and $0\lt \lvert R\rvert\lt 1$ is partially dissipative. Thus “maximally dissipative” describes an energy inequality, not necessarily zero reflection.

Kreiss and collaborators derive boundary and interior estimates for wave equations with Sommerfeld and more general tangential-derivative conditions in Kreiss et al. 2007, Theorems 1–3, pp. 3–10, PDF. Their hypotheses include the sign and spectral restrictions on the boundary operators; a visually similar condition outside those restrictions has no automatic estimate.

Radiative conditions at finite radius

For a three-dimensional outgoing scalar monopole,

$$\Phi=\frac{f(t-r)}{r},$$

the Sommerfeld condition

$$(\partial_t+\partial_r)(r\Phi)=0$$

is exact. Higher multipoles contain additional inverse powers of $r$, curvature backscatter creates tails, and nonlinear gravity couples modes. Consequently a homogeneous first-order Sommerfeld condition at $r=R$ reflects part of a generic outgoing solution. To specify both phase and normalization, consider the flat-space dipole

$$\Phi_\omega =\partial_z\left[ \frac{e^{-i\omega(t-r)}}r +\kappa_\omega\frac{e^{-i\omega(t+r-2R)}}r \right],$$

where $z$ is a Cartesian coordinate, $r=\sqrt{x^2+y^2+z^2}$, and $\omega\gt0$. Both terms solve the three-dimensional wave equation away from the origin. Imposing $(\partial_t+\partial_r)(r\Phi_\omega)=0$ at $r=R$ gives

$$\kappa_\omega =\frac{1}{2\omega^2R^2-2i\omega R-1} =O((\omega R)^{-2})$$

as $\omega R\to\infty$. This is the complex conjugate of Winicour 2012, §5.1, equations (62)–(63), pp. 48–49, whose time factor is $e^{+i\omega t}$. The sign of the imaginary term must change with the Fourier convention.

For a direct check, the outgoing term contributes $R^{-2}$ to the boundary residual after removing its common phase and $\cos\theta$, while the incoming term contributes $\kappa_\omega(-2\omega^2+2i\omega/R+R^{-2})$. Their sum vanishes for the displayed coefficient. The reflected-to-incident far-field power ratio is

$$|\kappa_\omega|^2=\frac{1}{1+4(\omega R)^4}.$$

It is $1/5$ when $\omega R=1$ and tends to one at low frequency. “Far away” must therefore be measured relative to wavelength, not only in coordinate units. This finite-radius error persists after grid refinement at fixed $R$.

For Einstein evolution, two physical incoming characteristic combinations can be related to Weyl curvature or metric derivatives, but gauge and constraint incoming fields must be treated simultaneously. Setting a nominal “incoming radiation scalar” to zero does not by itself make the full Einstein boundary condition stable or constraint preserving.

Reflective and AdS-like choices

A reflective boundary sends outgoing disturbances back into the domain. It can be conservative and analytically clean, which is appropriate for a physical wall or for common asymptotically AdS boundary conditions. It is inappropriate when a finite artificial boundary is meant to approximate an isolated system over times long enough for reflected radiation to return.

Reflectivity also depends on the variable. Dirichlet data for a metric representative need not mean a simple sign reversal of every gauge-invariant gravitational polarization. One should compute the frozen-coefficient scattering matrix for the physical characteristic fields and verify the energy flux in the formulation’s symmetrizer norm.

Dissipative and open-system choices

A strict contraction $\lVert H\rVert\lt 1$ in

$$U_{\mathrm{in}}=H U_{\mathrm{out}}+q$$

removes energy from homogeneous perturbations in the controlled norm. With $q\ne0$, the estimate bounds injected energy by the prescribed boundary data. This is useful for an artificial outer boundary, but dissipation can absorb physical, gauge, and constraint fields at different rates. It is not a proof of accurate radiation transmission.

At the continuum level, stability, constraint preservation, and absorption must all pass. At the discrete level, a summation-by-parts or equivalent boundary treatment must reproduce the estimate; a stable continuum condition can be destabilized by its implementation.

Decision table

Choice	Energy behavior	Physical interpretation	Main limitation
Reflective	Zero net homogeneous flux	Closed cavity or AdS-type wall	Returning radiation
Homogeneous Sommerfeld	Outward flux for model modes	Approximate isolation	Finite-radius reflection and backscatter
Inhomogeneous radiative	Controlled incoming signal	Exterior forcing	Boundary data must represent missing exterior
Strictly dissipative	Norm decreases for homogeneous data	Absorbing open boundary	May damp the wrong characteristic content

Common pitfalls

Stable is not nonreflecting. Both Dirichlet and Neumann conditions can obey an energy estimate while reflecting everything.

Zero incoming data are not no incoming physics globally. Exterior backscatter can legitimately return radiation through a finite worldtube.

Moving the boundary and refining the grid test different errors. Only the former diagnoses finite-radius modeling error.

Exercise: classify three scalar conditions

Using the half-line model, find $R$ for $u=0$, $u_x=0$, and $u_t+u_x=0$ at the boundary and classify their energy flux.

Solution

For time-preserved Dirichlet data, $u_t=0$, so $w_{\mathrm{in}}=-w_{\mathrm{out}}$ and $R=-1$: the flux vanishes and reflection is complete. Neumann data give $u_x=0$, hence $w_{\mathrm{in}}=w_{\mathrm{out}}$ and $R=1$, again complete reflection. The Sommerfeld condition is $w_{\mathrm{in}}=0$, so $R=0$ and $dE/dt=-w_{\mathrm{out}}^2/4$; it transmits the outgoing one-dimensional mode without reflection.

References

• Kreiss, Heinz-Otto, Oscar Reula, Olivier Sarbach, and Jeffrey Winicour. “Well-Posed Initial-Boundary Value Problem for the Harmonic Einstein Equations Using Energy Estimates.” Classical and Quantum Gravity 24, no. 23 (2007): 5973–5984. DOI; Open PDF.
• Winicour, Jeffrey. “Characteristic Evolution and Matching.” Living Reviews in Relativity 15, 2 (2012). DOI.