---
title: Discrete Normal Derivative Analysis
url: https://www.emergentmind.com/topics/discrete-normal-derivative
type: topic
---

# Discrete Normal Derivative Analysis

Searching arXiv for papers on discrete normal derivatives and closely related formulations.
A discrete normal derivative is a boundary-oriented discrete operator that approximates, represents, or generalizes the outward normal derivative in a discrete model. In the finite element analysis of elliptic PDEs, the term usually refers either to the classical boundary flux obtained from the discrete gradient on boundary elements or to a variational trace quantity defined by a discrete Green identity; for the Poisson problem on polygonal domains, both notions admit precise \(L^2(\Gamma)\) error estimates and are strongly influenced by corner singularities and mesh grading [1804.05723]. In broader usage, closely related constructions occur as directional derivative stencils on grids, graph- and vertex-based fluxes on fractals, and nonlocal spectral analogues on lattices, so the term denotes a class of boundary-sensitive discrete operators rather than a unique universal formula [1602.04791].

## 1. Continuous antecedent and weak boundary-flux interpretation

For the model problem
\[
-\Delta u=f \quad \text{in }\Omega,\qquad u=0 \quad \text{on }\Gamma=\partial\Omega,
\]
with \(\Omega\subset \mathbb R^2\) polygonal and \(f\in L^2(\Omega)\), the weak formulation is
\[
(\nabla u,\nabla v)_{L^2(\Omega)}=(f,v)_{L^2(\Omega)} \qquad \forall v\in H_0^1(\Omega).
\]
If \(u\in H^2(\Omega)\), the outward normal derivative is the classical boundary flux
\[
\partial_n u=\nabla u\cdot n \quad \text{on }\Gamma.
\]
For weak solutions on polygons, the normal derivative is understood through Green’s identity:
\[
(\partial_n u,w)_{L^2(\Gamma)}
=
(\nabla u,\nabla w)_{L^2(\Omega)}-(f,w)_{L^2(\Omega)}
\qquad \forall w\in H^1(\Omega).
\]
In the weighted regularity regime used for polygonal domains, this identity yields \(\partial_n u\in L^2(\Gamma)\) [1804.05723].

The regularity of \(\partial_n u\) is controlled by corner singularities. If the corners are \(\boldsymbol c_j\) with interior angles \(\omega_j\), the singular exponents are
\[
\lambda_j=\frac{\pi}{\omega_j},\qquad \bar\lambda=\min_j\lambda_j.
\]
These exponents limit the global regularity of both \(u\) and \(\partial_n u\). In weighted Sobolev spaces \(V_{\vec\beta}^{2,p}(\Omega)\), one obtains precise corner-dependent regularity, and pointwise normal derivatives may vanish at convex corners and may blow up at non-convex ones [1804.05723].

A complementary boundary-kernel viewpoint is furnished by Green functions. For a smooth bounded domain \(D\subset\mathbb C\), with Green function \(g_w\) and Poisson kernel
\[
P(z,\zeta)=\partial_{n_\zeta} g_z(\zeta),
\]
the normal derivative of a function \(f\in C^1(\overline D)\cap C^2(D)\) satisfying \(f=0\) on \(\partial D\) obeys
\[
\partial_n f(\zeta)=\int_D \Delta f(z)\,P(z,\zeta)\,dA(z).
\]
This identifies the normal derivative as a boundary functional generated by an interior operator and a boundary kernel [1209.5363].

## 2. Finite element realizations

For a conforming triangulation \(\mathcal T_h\), the standard linear finite element spaces are
\[
V_h=\{v_h\in C(\overline\Omega): v_h|_T\in\mathcal P_1(T)\ \forall T\in\mathcal T_h\},
\qquad
V_{0h}=V_h\cap H_0^1(\Omega),
\]
and the discrete state \(u_h\in V_{0h}\) solves
\[
(\nabla u_h,\nabla v_h)_{L^2(\Omega)}=(f,v_h)_{L^2(\Omega)}
\qquad \forall v_h\in V_{0h}.
\]
Within this setting, two distinct discrete normal derivatives are used [1804.05723].

| Variant | Definition | Representation |
|---|---|---|
| Standard discrete normal derivative | \(\partial_n u_h|_E=\nabla u_h|_T\cdot n_E\) on a boundary edge \(E\subset\Gamma\) of element \(T\) | elementwise constant on boundary edges |
| Variational discrete normal derivative | \((\partial_n^h u_h,w_h)_{L^2(\Gamma)}=(\nabla u_h,\nabla w_h)_{L^2(\Omega)}-(f,w_h)_{L^2(\Omega)}\) for all \(w_h\in V_h\) | \(\partial_n^h u_h\in V_h^\partial=\operatorname{Tr}(V_h)\) |

The standard derivative is the natural elementwise flux. Because \(u_h\) is piecewise linear, \(\nabla u_h\) is piecewise constant, so the boundary flux on each boundary edge is immediately available from the adjacent element.

The variational discrete normal derivative is the \(L^2(\Gamma)\)-Riesz representation of the discrete boundary flux functional
\[
V_h\ni w_h\mapsto (\nabla u_h,\nabla w_h)_{L^2(\Omega)}-(f,w_h)_{L^2(\Omega)}.
\]
Its coefficients are obtained from a boundary mass matrix system: with \(\{\varphi_i\}\) the boundary nodal basis,
\[
\sum_j (\varphi_j,\varphi_i)_{L^2(\Gamma)}\,(\partial_n^h u_h)_j
=
(\nabla u_h,\nabla \varphi_i)_{L^2(\Omega)}-(f,\varphi_i)_{L^2(\Omega)}.
\]

The key error identity is
\[
(\partial_n u-\partial_n^h u_h,w_h)_{L^2(\Gamma)}
=
(\nabla(u-u_h),\nabla w_h)_{L^2(\Omega)}
\qquad \forall w_h\in V_h.
\]
This identity is the starting point for the \(L^2(\Gamma)\) analysis of the variational flux and explains why \(\partial_n^h u_h\) is especially natural in variational boundary formulations [1804.05723].

## 3. Mesh dependence and error estimates

On general shape-regular quasi-uniform meshes, the convergence of discrete normal derivatives is limited by corner regularity. If \(\omega=\max_j\omega_j\), then for both \(\partial_n u_h\) and \(\partial_n^h u_h\) the known quasi-uniform rates are \(s=1\) up to logarithms when \(\omega<2\pi/3\), and otherwise \(s<\pi/\omega-1/2\) [1804.05723].

The principal refinement strategy in the finite element analysis is boundary concentration. Writing
\[
\rho_T=\operatorname{dist}(T,\Gamma),
\]
the mesh grading condition is
\[
h_T\sim
\begin{cases}
h^2, & \rho_T=0,\\[0.3em]
h\sqrt{\rho_T}, & \rho_T>0.
\end{cases}
\]
Thus, elements touching the boundary have diameter \(h^2\), while elements at distance \(O(1)\) from the boundary have size \(\sim h\). The total number of elements is \(O(h^{-2}|\ln h|)\), only a logarithmic factor above the \(O(h^{-2})\) complexity of quasi-uniform meshes [1804.05723].

The analysis uses weighted errors involving
\[
\sigma(x)=\rho(x)+d_I,\qquad \rho(x)=\operatorname{dist}(x,\Gamma),\qquad d_I\sim h^2,
\]
together with weighted regularity, dyadic decomposition by boundary distance, localized energy estimates, and trace inequalities. For the standard boundary flux, a local estimate on a boundary edge \(E\) with adjacent element \(T\) is
\[
\|\partial_n e_h\|_{L^2(E)}
\le
c\big(h_T^{-3/2}\|u-u_h\|_{L^2(T)}+h_T^{1/2}|u|_{H^2(T)}\big),
\qquad e_h=u-u_h.
\]

The resulting global rates on boundary-concentrated meshes are substantially sharper. In convex polygons with \(f\in C^{0,\sigma}(\overline\Omega)\), the standard derivative satisfies
\[
\|\partial_n(u-u_h)\|_{L^2(\Gamma)}
\le
c\,h^{\min\{2,\,-1+2\bar\lambda-2\varepsilon\}}|\ln h|^{3/2}\,\|f\|_{C^{0,\sigma}(\overline\Omega)},
\]
while the variational derivative satisfies
\[
\|\partial_n u-\partial_n^h u_h\|_{L^2(\Gamma)}
\le
c\,h^{\min\{2,\,-1+2\bar\lambda-\varepsilon\}}|\ln h|^{3/2}\,\|f\|_{C^{0,\sigma}(\overline\Omega)}.
\]
In non-convex polygons, the best-case order drops to \(1\) in the corresponding estimates. In particular, if \(\omega<2\pi/3\), then boundary concentration yields \(s=2\) up to logarithms, whereas quasi-uniform meshes yield only \(s=1\) up to logarithms. This is the sense in which boundary-concentrated meshes double the order in \(h\) for normal-derivative approximation [1804.05723].

A common source of confusion is that this grading is boundary-concentrated, not corner-concentrated. The refinement law is governed by distance to \(\Gamma\), not only by distance to re-entrant corners [1804.05723].

## 4. Boundary kernels, perturbations, and control problems

Normal derivatives enter directly into boundary representation formulas. For the Laplacian, the Poisson kernel is
\[
P(w,\zeta)=\partial_{n_\zeta} g_w(\zeta),
\]
and the solution of the homogeneous Dirichlet problem is represented by
\[
u(w)=\int_{\partial D} f(\zeta)\,\partial_{n_\zeta}g_w(\zeta)\,ds(\zeta).
\]
This identifies the normal derivative of the Green function as a boundary density reproducing interior harmonic values [1209.5363].

The same boundary-kernel viewpoint persists under operator perturbation. For the Schrödinger operator \(\Delta-\varepsilon u\), the outward normal derivative of the Green function satisfies
\[
\partial_{n_\zeta} g_w^*(\zeta)
=
\partial_{n_\zeta} g_w(\zeta)
+
\varepsilon\int_D u(z)\,g(z,w)\,P(z,\zeta)\,dA(z)
+
o(\varepsilon).
\]
For the Laplace–Beltrami operator \(L=\nabla(\lambda\nabla)\) with \(\lambda=1+\varepsilon u\), the first variation contains a bulk term involving \(\Delta u\) and local multiplicative corrections proportional to \(\partial_n g_w(\zeta)\) [1209.5363]. This operator-theoretic structure explains why variational discrete normal derivatives often appear as natural boundary unknowns.

A direct application is Dirichlet boundary control. With desired state \(y_d\in L^2(\Omega)\), the unconstrained problem is
\[
J(y,u)=\frac12\|y-y_d\|_{L^2(\Omega)}^2+\frac\alpha2\|u\|_{L^2(\Gamma)}^2\to\min
\]
subject to
\[
-\Delta y=0 \text{ in }\Omega,\qquad y=u \text{ on }\Gamma.
\]
The continuous optimality system includes
\[
-\Delta p=y-y_d \text{ in }\Omega,\qquad p=0 \text{ on }\Gamma,\qquad \alpha u-\partial_n p=0 \text{ on }\Gamma,
\]
so the control is \(u=\alpha^{-1}\partial_n p\). In the discrete system,
\[
(\alpha u_h-\partial_n^h p_h,w_h)_{L^2(\Gamma)}=0
\qquad \forall w_h\in V_h^\partial,
\]
hence \(\alpha u_h=\partial_n^h p_h\) in \(L^2(\Gamma)\). The control error is then limited by the discrete variational normal derivative, and on boundary-concentrated meshes one obtains
\[
\|u-u_h\|_{L^2(\Gamma)}+\|y-y_h\|_{L^2(\Omega)}
\le
c\,h^{\min\{2,\,-1+2\bar\lambda-2\varepsilon\}}|\ln h|^{3/2}
\]
under the regularity assumptions stated for the theorem [1804.05723].

## 5. Grid-based directional constructions

On structured grids, a discrete normal derivative is often formed from discrete partial derivatives and a normal vector. One high-order construction is the Vandermonde-based discrete differential operator. Using \(N+1\) samples \(x_i=x_0+h_i\), one solves
\[
W_N
\begin{bmatrix}
a_0^e\\
a_1^e\\
\vdots\\
a_N^e
\end{bmatrix}
=
\begin{bmatrix}
f(x_1)\\
f(x_2)\\
\vdots\\
f(x_{N+1})
\end{bmatrix},
\]
where \(W_N\) is the Vandermonde matrix, and sets \((f^{(n)}(x_0))^e=n!\,a_n^e\). The theorem gives
\[
\bigl|{f^{(i)}(x_0)}^e-f^{(i)}(x_0)\bigr|
\le
M\,K^{2N+1-i}\,h^{N+1-i}\,(N-i)!.
\]
In two dimensions, the estimated gradient components \(f_x^e\) and \(f_y^e\) are extracted from the coefficient vector, and the discrete normal derivative is then
\[
\left(\frac{\partial f}{\partial n}(x_0,y_0)\right)^e
=
f_x^e(x_0,y_0)\,n_x+f_y^e(x_0,y_0)\,n_y.
\]
This construction computes all derivatives up to the chosen order simultaneously and makes the normal derivative a post-processing of the discrete gradient [2507.09480].

A scale-space formulation uses Gaussian derivative operators. For a unit normal \(n=(n_x,n_y)\),
\[
\partial_n L=n_xL_x+n_yL_y,
\qquad
\partial_{nn}L=n_x^2L_{xx}+2n_xn_yL_{xy}+n_y^2L_{yy}.
\]
The analysis of hybrid discretizations distinguishes normalized sampled Gaussian plus central differences and integrated Gaussian plus central differences from direct sampled or integrated Gaussian derivatives and from the genuinely discrete Gaussian \(T_{\mathrm{disc}}(n;s)=e^{-s}I_n(s)\). The reported conclusions are that direct sampled or integrated derivative kernels match continuous spatial spread substantially better than the hybrid methods, while genuinely discrete derivatives behave better than hybrids; hybrid methods remain attractive because one smoothing stage can be shared across multiple derivative orders [2405.05095].

A separate directional framework is Tao General Difference. Its first-order finite-window operator is
\[
f_{\mathrm{TGD}}'(x_0;w,W)
=
\frac{\int_0^W f(x_0+t)w(t)\,dt-\int_0^W f(x_0-t)w(t)\,dt}
{2\int_0^W t\,w(t)\,dt},
\]
and in multiple dimensions the directional operator takes the convolution form
\[
\frac{\partial_{\mathrm{TGD}} f}{\partial \mathbf v}(\mathbf x)
=
K_1(\widetilde T_{\mathbf v}*f)(\mathbf x).
\]
Choosing \(\mathbf v\) to be the unit normal yields a discrete normal derivative. The orthogonal construction replaces the derivative kernel in the normal direction by smoothing kernels in tangential directions, for example
\[
\widetilde T_{x\perp}(x,y)=T(x)\,S(y),
\]
which is particularly suited to axis-aligned boundaries on Cartesian grids [2305.08098].

## 6. Fractal and lattice analogues, and conceptual distinctions

On p.c.f. fractals, discrete normal derivatives appear as one component of Strichartz’s derivatives at vertices. For a boundary vertex \(v_j\),
\[
d_{jk}f(v_j)=\lim_{m\to\infty}\lambda_{jk}^{-m}\,\beta_{jk}\bigl(f|_{F_j^mV_0}\bigr),
\qquad 2\le k\le N_0,
\]
and the case \(k=2\) is, up to normalization, the normal derivative. Here \(\beta_{jk}\) and \(\lambda_{jk}\) are left eigenvectors and eigenvalues of the local harmonic extension matrix \(M_j\). At junction vertices, the branchwise normal derivatives satisfy the compatibility condition
\[
\sum_{j\in J(x)} d_{j'2}f(x)=0.
\]
For functions in \(\operatorname{dom}(\Delta_\mu)\), normal derivatives are uniformly bounded on all vertices. If the normal derivative vanishes at a fixed vertex, then the normal derivatives at neighboring vertices decay to zero with explicit rates governed by the comparison among \(r_j\mu_j\) and \(|\lambda_{j3}|\), with the three cases
\[
O(\mu_j^m),\qquad O(m\mu_j^m),\qquad O\big((\lambda_{j3}r_j^{-1})^m\big)
\]
according to whether \(r_j\mu_j\) is greater than, equal to, or less than \(|\lambda_{j3}|\) [1602.04791].

On \(\mathbb Z^N\), a different, nonlocal notion arises from the derivative at \(s=0\) of the fractional discrete Laplacian. In one dimension,
\[
\log(-\Delta_1)f(n)
=
-\sum_{m\neq n}\frac{f(m)}{|n-m|},
\]
and more generally
\[
\log(-\Delta_N)f(n)
=
-\sum_{m\neq n}K(n-m)f(m)+\rho_N f(n),
\qquad K(m)\sim |m|^{-N}.
\]
The source explicitly interprets this as a discrete logarithmic Laplacian and a nonlocal boundary/normal derivative-type operator associated with the discrete Laplacian and its fractional powers [2603.01039]. This usage is not equivalent to the finite element boundary flux, but it extends the same boundary-oriented vocabulary into a spectral and nonlocal setting.

The main conceptual distinction is therefore structural. In finite elements, the discrete normal derivative is a boundary flux tied to Green’s identity and trace spaces. On structured grids, it is typically a directional derivative projected onto a normal vector. On p.c.f. fractals, it is a renormalized limit of graph-level boundary combinations. On lattices with fractional operators, it can become a nonlocal spectral derivative-type operator. These objects are analogous because each one encodes outward flux, directional variation, or boundary response in a discrete medium, but they are not interchangeable.

Source: https://www.emergentmind.com/topics/discrete-normal-derivative