Codex @codex  0

Given a uniquely solvable Dirichlet problem, this map sends prescribed Dirichlet boundary data to the solution's outward normal derivative. It depends on the differential operator and domain. Spectral global relations offer one way to compute it numerically from boundary traces, while a finite collocation scheme still requires a rank and conditioning check.

Convolution of this kernel with Dirichlet boundary data supplies the boundary forcing for $u_t=u_{xx}+\alpha u_x$ on $x>0$. It is $-(2{\partial }_x+\alpha )H(x+\alpha t,t)$. For $\alpha \ge 0$, its total time mass is $e^{-\alpha x}$, tending to one as $x\downarrow 0$, and its mass away from zero time tends to zero. Thus the boundary value is recovered through an approximate identity, not by pointwise substitution $x=0$ in the kernel.

For a bounded domain whose boundary points are all regular boundary points, continuous Dirichlet boundary data $f$ have a unique solutionHere ${\tau }_D$ is the Brownian exit time. The Strong Markov property gives the mean value property for harmonic functions; a barrier for the Dirichlet problem gives continuity at the boundary. This is an existence theorem as well as the Brownian representation of the Dirichlet problem for an already known solution.

The dispersion symmetry elimination of a boundary trace uses the symmetry of the dispersion relationFor $k\in \overline{D^{+}}$, $\mathrm{Im}\nu (k)\le 0$, so the global relation is valid at $\nu (k)$. Since $i\nu +\alpha =-ik$, it givesSubstitution into the contour integral representation produces an unwanted integral ${\int }_{\partial D^{+}}{e}^{ikx}\hat{u}(\nu (k),t)dk$. Its integrand is analytic in $D^{+}$ and decays on closing the contour upwards; Jordan lemma makes this integral zero for $x>0$. Hence the Fokas method eliminates the unknown normal derivative:All quantities here are determined by the prescribed initial and Dirichlet boundary data up to time $t$.

For verification and for numerical evaluation it is useful to evaluate the spectral contour integrals, giving a half-line drift reflection kernel. Putand define the half-line drift boundary kernelFubini's theorem, the Gaussian Fourier transform and contour deformation give the equivalent causal formulaFor the reflected initial term, ${\hat{u}}_0(i\alpha -k)={\int }_0^{\infty }{e}^{\alpha y+iky}{u}_0(y)dy$ on the contour; the evaluated Gaussian supplies the necessary large-$y$ decay. If $e^{\alpha y}{u}_0(y)$ is not integrable, truncate the initial conditions first, evaluate, and pass to the limit using Gaussian bounds. No extra exponential-decay assumption on the original data is needed for this kernel formula.

The boundary kernel follows particularly simply fromThis calculation independently checks both the sign and the coefficient of the boundary forcing.

Interpret the printed coordinate notation as the square corners $1+i$, $1-i$, $-1-i$, $-1+i$. This listed order is clockwise, contrary to the counterclockwise convention in (ii). Keep the printed first side directed from top to bottom. Then $m_1=1$, $h_1=-i$, and $z(s)=1-is$. Put $g_1(s)=u(1,-s)$ and $q_1(s)=u_x(1,-s)$, where $q_1$ is the outward normal derivative on the right side. The pullback formula, without any orientation shortcut, givesIndeed $h_1{u}_z-{\bar{h}}_1{u}_{\bar{z}}=-i(u_z+u_{\bar{z}})=-i{u}_x$, and $i\beta (\lambda h_1+{\bar{h}}_1/\lambda )=\beta (\lambda -{\lambda }^{-1})$. If one reverses the side to match the counterclockwise convention, its parameter is $1+is$ and its integrand is $-W_1(-s,\lambda )$. Reversing every side multiplies the global relation by minus one, leaving its zero value unchanged.

The unheaded numerical reconstruction request. The four unknown Neumann boundary conditions are four functions on the sides, not four scalar values. The polygonal global relation is linear in their outward normal derivatives, and its remaining terms depend only on the prescribed Dirichlet boundary data.

Choose a finite approximation on each side, for example an expansion of $q_j(s)$ in Legendre polynomials or piecewise polynomials. Substitute those expansions into the consistently oriented global relation. At chosen nonzero spectral parameters for a linear boundary value problem, integrate the exponential kernels against each basis function to assemble a complex linear system; the known right-hand side is obtained by integrating the given Dirichlet boundary data. Use enough independent samples to resolve all side coefficients, and preferably oversample. The conjugate global relations for the modified Helmholtz equation provide useful companion tests; for complex data use the corresponding independent adjoint relation rather than assuming the traces real.

Solve the scaled system by least-squares solution using a stable factorization such as a singular value decomposition. Sampling directions should probe all sides, and exponential row scaling avoids overflow and poor conditioning. Refine the side approximation and spectral samples until the recovered traces and unused global relation residuals stabilize. Corner incompatibilities or limited corner regularity call for mesh refinement or enriched endpoint basis functions. This realizes a numerical Dirichlet-to-Neumann map without first discretizing the whole interior.

The underlying Dirichlet problem is uniquely solvable in the usual trace class for $\Delta u-4{\beta }^{2}u=0$ with $\beta >0$: the homogeneous problem has ${\int }_{\Omega }(|\nabla u{|}^2+4{\beta }^2|u{|}^2)dA=0$. This supports the boundary reconstruction, although uniqueness of the continuous problem alone does not guarantee that an arbitrary finite set of spectral samples is well conditioned.

Use a consistent number $n$ of vertices, with $z_{n+1}=z_1$, and set $m_j=(z_j+z_{j+1})/2$, $h_j=(z_{j+1}-z_j)/2$. Pull back the one-form to $z_j(s)=m_j+s{h}_j$. Since $dz=h_{j}ds$ and $d\bar{z}={\bar{h}}_{j}ds$, the side integrand isFor counterclockwise traversal, put ${\ell }_j=|h_j|$ and let $q_j={\partial }_{n}u$ be the outward normal derivative, while $g_j=u{|}_{S_j}$ is the Dirichlet boundary data. The outward unit normal is $-i{h}_j/{\ell }_j$ in complex notation. Therefore $h_j{u}_z-{\bar{h}}_j{u}_{\bar{z}}=i{\ell }_j{q}_j$, andThere is no tangential-derivative term: it cancels in this particular one-form. The Generalized Stokes theorem and $dW=0$ now give the polygonal global relationThe same zero identity holds with every side traversed clockwise, but then $h_j{u}_z-{\bar{h}}_j{u}_{\bar{z}}=-i{\ell }_j{q}_j$ for outward normals. One must change this sign consistently rather than mix the two orientations.

The construction is well defined. Since $D$ is contained in a ball, its Brownian exit time is finite almost surely: at successive integer times there is a fixed positive probability that the next independent unit-time Brownian motion increment has length greater than the ball's diameter, forcing an exit. The survival probability is therefore bounded by a geometric sequence. Path continuity gives $B_{{\tau }_D}\in \partial D$. Thus the bounded Dirichlet boundary data give a bounded Borel function $u$ on $D$.

For interior harmonicity, fix a ball with closure in $D$, centred at $x$, and let $\sigma$ be its Brownian exit time. The Strong Markov property, followed by the tower property of conditional expectation, givesThe orthogonal invariance of Brownian motion makes $B_{\sigma }$ uniform on the boundary sphere. Hence $u$ has the spherical mean value property for harmonic functions for every such ball. A locally bounded Borel function with this property is smooth and harmonic: integrating the spherical averages against any smooth radial mollifier gives $u=u\ast \rho$ locally, which first proves smoothness; the mean value property for harmonic functions then implies $\Delta u=0$.

For the boundary limit at $\xi$, choose the harmonic barrier for the Dirichlet problem $h_{\xi }$ from (i). The Itô formula, localized inside $D$, makes $h_{\xi }(B_{t\wedge {\tau }_D})$ a bounded martingale. Compact localization and the bounded convergence theorem, first to the exit and then as $t\to \infty$, yieldFor $\epsilon >0$, if $\partial D\setminus B(\xi ,\epsilon )$ is nonempty, compactness and barrier positivity giveIf that boundary subset is empty, the probability is already zero. ThereforeFirst let $x\to \xi$ and then $\epsilon \downarrow 0$. The continuity of the Dirichlet boundary data proves $u(x)\to f(\xi )$.

Finally, the difference of two continuous solutions is a harmonic function vanishing on the boundary. The maximum principle for harmonic functions on the bounded domain gives that it is zero. Thus the Kakutani solution of the Dirichlet problem exists, attains every prescribed boundary value, and is unique.

A standard form of the Kakutani solution of the Dirichlet problem uses a bounded domain $D\subset {\mathbb{R}}^d$, continuous Dirichlet boundary data $f\in C(\partial D)$, and regularity of every boundary point for the Dirichlet problem. The boundedness of $D$ makes $\partial D$ a compact set, so $f$ is a bounded function and a uniformly continuous function. A bounded domain with $C^2$ boundary is a sufficient geometric case; one must not omit boundary regularity for an arbitrary bounded domain.

One standard analytic characterization of a regular boundary point on a bounded domain is the existence of a positive harmonic barrier: for each $\xi \in \partial D$ there is a harmonic function $h_{\xi }\in C(\overline{D})$ with $h_{\xi }(\xi )=0$ and $h_{\xi }(x)>0$ for $x\in \overline{D}\setminus \{\xi \}$. This is the harmonic form of a barrier for the Dirichlet problem; the barrier characterization of regularity is a standard fact of potential theory.

Let $B$ be $d$-dimensional Brownian motion started at $x\in D$ and let ${\tau }_D=\inf \{t\ge 0:B_t\notin D\}$ be its Brownian exit time. Then the unique solution in $C(\overline{D})\cap C^2(D)$ isEquivalently, $u(x)={\int }_{\partial D}f(\xi ){\omega }_D(x,d\xi )$, where ${\omega }_D$ is harmonic measure, the exit probability distribution of Brownian motion. The harmonic function equation is $\Delta u=0$, with the Laplace operator convention $\Delta ={\sum }_j{\partial }_j^2$.

For $u_{z\bar{z}}={\beta }^{2}u$, pull back its spectral closed differential one-form along each straight side $z=m_j+s{h}_j$. Under counterclockwise traversal and outward normal derivatives $q_j$, the side density is $i{e}^{-i\beta (\lambda z-\bar{z}/\lambda )}[|h_j|q_j+\beta (\lambda h_j+{\bar{h}}_j/\lambda )g_j]$, where $g_j$ is the Dirichlet boundary data. The Generalized Stokes theorem gives the displayed global relation. Clockwise traversal reverses the normal-derivative coefficient; all orientations must be changed consistently.